Mach–Zehnder interferometers incorporating electrochromic molecules for controlled single-photon detection

One particularly fruitful research in the fields of integrated photonics, carried out by a good number of physicists and engineers, concerns the study of different types of materials to be used to control the detection of photons, if not just single photons, in interferometers. In a Mach–Zehnder interferometer, which consists of two beam-splitters, two mirrors, and two detectors, a material that can cause a controlled change in the phase of light in one of the two arms of the interferometer consequently allows control of the probability of detection at the two detectors. In this work, we use an electrochromic molecule, N,N′-bis(cysteine)pyromellitic diimide (BCPD), that has a refractive index dependent on the applied electric field. We simulate the single-photon detection probability in a Mach–Zehnder interferometer with direct light transmission and a waveguide-based Mach–Zehnder interferometer, consisting of two 3-dB couplers connected by two optical channel waveguides. With the employment of the non-equilibrium Green’s function formalism, we have simulated the conductance of BCPD. The results could be of interest in quantum communication.

Introduction

Electrochromic molecules are materials that can reversibly change their optical properties, such as color or transparency, through redox reactions triggered by an external electrical potential. This phenomenon, known as electrochromism, has gained great interest due to its potential applications in energy-efficient smart windows, displays, wearable electronics, and multifunctional energy-storage devices (Pradhan et al., 2025; Chen et al., 2025; Cho et al., 2025). Among the electrochromic molecules, viologens and viologen-based materials are the most intensely studied species (Bird and Kuhn, 1981; Stolar, 2020; Ambrose et al., 2024; Otaegui et al., 2025). Other noteworthy electrochromic molecules are graphene molecules, also called polycyclic aromatic hydrocarbons (Ji et al., 2015; Stec et al., 2017; Chapkin et al., 2020). Moreover, the small molecule N,N′-bis(cysteine)pyromellitic diimide (BCPD) is a good example of a short molecule that exhibits a distinct and reversible color change (Stolar, 2020; Abdinejad et al., 2017; Abdinejad et al., 2018; Guo et al., 2014).

The color change, which is related to a change in the complex refractive index of the material, can also be exploited as a phase shift component in a Mach–Zehnder interferometer (Rioux, 2019). Mach–Zehnder interferometers with photochromic molecules or polymers—and with vanadium dioxide exploiting its insulating-to-metallic phase change—have been reported (Scotognella, 2024a; Scotognella, 2024b). Mach–Zehnder interferometers that are based on direct transmission of light are composed of a light source, two beam-splitters, two mirrors, and two detectors (Rioux, 2019). Mach–Zehnder interferometers that are based on channel waveguides, e.g., in a glass slide, are composed of a light source, two 3-dB couplers that have the same function as the two beam-splitters, and two detectors (De Ridder et al., 2006; Ng et al., 2012).

In this study, we use an electrochromic molecule, namely, BCPD (Stolar, 2020; Abdinejad et al., 2017; Abdinejad et al., 2018), as the material in both arms of a Mach–Zehnder interferometer for direct transmission of light and in a Mach–Zehnder interferometer composed of channel waveguides, e.g., buried in a glass slab. We simulate the application of the electric field on BCPD via the addition of an electron charge to the molecule, with a subsequent change of its complex refractive index. The electric field-induced change in the BCPD refractive index in one of the arms of the interferometers allows control of the photon detection in the two detectors of the apparatus. Moreover, by using the non-equilibrium Green’s function formalism, we found the conductance of BCPD.

Methods

To study the optical response of BCPD, we employed the Hückel theory. Specifically, with this model, we simulated the absorption bands of the BCPD molecule in its neutral its charged states. The Hückel theory is an established molecular orbital theory in which the molecular orbitals are written as the linear combination of atomic orbitals, in which the nuclei positions are fixed (this is a strong approximation, and more accurate calculations can go beyond this approximation (Janoš et al., 2025), but such approaches are not considered in this work) and the molecular orbitals are linear combinations of carbon p_z orbitals (other elements can also be considered), neglecting electron–electron interactions (Guy and Troy, 2022). Equation 1 presents the time-independent Schrödinger equation:

$H\psi =E\psi .(1)$

In this study, the HuLiS package (Huckel Theory and HuLiS, 2013) was employed to determine the Hamiltonian for BCPD (neutral and charged). Then, eigenvalues and eigenvectors were calculated. To find the transition probabilities, we employed the Fermi golden rule, as reported in Equation 2:

${\mathrm{\Gamma }}_{i\to j}=\frac{2\pi }{\hslash }{\left|〈{\psi }^j\left|H\right|{\psi }^i〉\right|}^2.(2)$

Using Equation 2, the transition probabilities were determined between the highest occupied molecular orbital (HOMO) and the unoccupied orbitals. The simulated absorption spectra were built using a set of Gaussian peaks, with a peak linewidth c of 0.1 eV. Considering the transitions from the ground ith state to the different j^th excited states, the absorption coefficient can be written as shown in Equation 3:

$\alpha \left(E\right)={\sum }_j{\mathrm{\Gamma }}_{i\to j}\mathit{exp}\left(\frac{{\left(E-\left(E_j-E_i\right)\right)}^2}{{2c}^2}\right).(3)$

For clarity, we emphasize that Hückel’s theory finds six unoccupied states for BCPD. Therefore, the simulated absorption spectrum is given by the six transitions from HOMO to the six unoccupied states. For the calculations based on the density functional theory (DFT), the following procedure was followed. The molecules were designed using the Avogadro package (Hanwell et al., 2012). We optimized the geometries and calculated the electronic transitions via density functional theory with the ORCA package developed by Neese (2012). In these calculations, we used the B3LYP functional (Lee et al., 1988). Moreover, we employed the Ahlrichs split valence basis set (Schäfer et al., 1992) together with the all-electron nonrelativistic basis set SVPalls1 (Schäfer et al., 1994; Eichkorn et al., 1997). Finally, for these calculations, we used the Libint library (Valeev, 2014) and the Libxc library (Lehtola et al., 2018; Marques et al., 2012).

Using Equation 4, modified from the equation reported by Kohandani and Saini (2022), the imaginary part of the complex refractive index can be extracted from the absorption coefficient:

$k\left(E\right)=\frac{A\alpha \left(E\right)}{4\pi }.(4)$

In Equation 4, A is a parameter (its value is $1\times {10}^{15}$). Another formalism for determining the imaginary part of the refractive index was reported by Dresselhaus (2017). Notably, the parameter A has been chosen to reproduce the experimental data. The real part of the complex refractive index can be derived via the Kramers–Kronig relations (Pankove, 1975) (Equation 5):

$n\left(\overline{\nu }\right)-1=\frac{2}{\pi }\mathcal{P}{\int }_0^{\infty }\frac{\omega k\left(\omega \right)}{{\omega }^2-{\overline{\nu }}^2}d\omega .(5)$

The Kramers–Kronig relations are actually Hilbert transforms, as in the formalism reported by Ogilvie and Fee (2013). For the real part of the refractive index, a constant offset of 1.5 has been used, as detailed by Wiebeler et al. (2014) and Scotognella (2020). Such an offset is used to consider the experimental factors, such as light scattering. Thus, Equation 6 can be used to describe the complex refractive index of BCPD in the neutral (n_n) and charged (n_c) forms:

$n_{n,c}^{complex}\left(E\right)=n_{n,c}\left(E\right)+i{k}_{n,c}\left(E\right).(6)$

In the model, we simulated the contribution of the electric field applied to the molecule by studying the optical properties of the charged molecule, i.e., the molecule to which an electron has been added. The change in the electronic configuration of the molecule, i.e., the different filling of the molecular orbitals in the charged molecule compared to that in the neutral molecule, gives rise to different absorption bands, with a consequent variation in the complex refractive index of the molecule.

The non-equilibrium Green’s function formalism was employed for the calculation of the elastic transmission, following Solomon et al. (2011). Assuming that only a single site of the molecule couples to each electrode, for BCPD, the vector for the left electrode is provided in Equation 7:

$V_L=\left[\begin{array}{cccccccccccccccc}\gamma & 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0\end{array}\right].(7)$

On the other hand, the vector for the right electrode is provided in Equation 8:

$V_R=\left[\begin{array}{cccccccccccccccc}0& 0& 0& 0& 0& 0& 0& 0& 0& 0& \gamma & 0& 0& 0& 0& 0\end{array}\right].(8)$

Following the sketch of BCPD in Figure 1, the non-zero element of V_L is the first one, representing the connection between the nitrogen atom of BCPD indexed “1” and the electrode (which will replace the hydrogen atom). Similarly, the non-zero element of VR is the 11th one, representing the connection between the nitrogen atom of BCPD indexed “11” and the electrode (which also replaces the hydrogen atom next to the nitrogen atom). The value of $\gamma$ is −1 eV, as reported by Solomon et al. (2011) for different studied molecule/electrode junctions. A similar value was used for another molecular/electrode system by Chang et al. (2014). The broadening function is provided in Equation 9:

${\mathrm{\Gamma }}^{L\left(R\right)}=2\pi \rho V_{L\left(R\right)}{V}_{L\left(R\right)}^{†}.(9)$

Figure 1

Figure 1. Chemical sketch of N,N′-bis(cysteine)pyromellitic diimide. The atoms are numbered as in the Hückel Hamiltonian (the four hydrogen atoms are not numbered).

Here, $\rho$ is the density of state of the electrode, which is set to 1/2π (eV)⁻¹ (following the formalism presented by Solomon et al. (2011)). With the determination of the broadening function from Equation 9, we are able to find the tunneling self-energy, which is purely of imaginary value, using Equation 10:

${\mathrm{\Sigma }}_T^{L\left(R\right)}=-\frac{i}{2}{\mathrm{\Gamma }}^{L\left(R\right)}.(10)$

In Equation 11, the energy-dependent retarded Green’s function is obtained as follows:

$G^r\left(E\right)={\left(E-H_{mol}-{\mathrm{\Sigma }}_T^L-{\mathrm{\Sigma }}_T^R\right)}^{-1}.(11)$

It is noteworthy that the advanced Green’s function $G^a\left(E\right)$ is the conjugated transpose of the retarded Green’s function. Using the broadening functions and Green’s functions, it is possible to determine the energy-dependent elastic transmission, as provided in Equation 12:

$T\left(E\right)=Tr\left[{\mathrm{\Gamma }}^L{G}^r\left(E\right){\mathrm{\Gamma }}^R{G}^a\left(E\right)\right].(12)$

The transmission is related, through Equation 13, to the conductance by the relation (Chang et al., 2014).

$G=G_{0}T\left(E_F\right).(13)$

Here, $G_0$ is the quantum of conductance (with value $7.748\times {10}^{-5}S$).

For the Mach–Zehnder interferometer with direct transmission of light, the state vectors and operators are implemented in agreement with Rioux (2019). Equation 14 describes the photon moving along the x-axis as the state vector:

$x=\left[\begin{array}{c}1\\ 0\end{array}\right].(14)$

Meanwhile, in Equation 15, the photon moving on the y-axis is represented by the state vector:

$y=\left[\begin{array}{c}0\\ 1\end{array}\right].(15)$

The two mirrors are represented by the operator M (Equation 16):

$M=\left[\begin{array}{cc}0& 1\\ 1& 0\end{array}\right].(16)$

Equation 17 describes the two beam-splitters that are represented by the operator BS:

$BS=\frac{1}{\sqrt{2}}\left[\begin{array}{cc}1& i\\ i& 1\end{array}\right].(17)$

A 90° phase shift is assigned to the reflection. The BCPD-based electrochromic (BCPD) layer is represented by the operator (Equation 18):

$BCPD=\left[\begin{array}{cc}\mathit{exp}\left(i\frac{2\pi d{n}_{n,c}^{complez}}{\lambda }\right)& 0\\ 0& \mathit{exp}\left(i\frac{2\pi d{n}_{n,c}^{complez}}{\lambda }\right)\end{array}\right].(18)$

Here, d is the thickness of the layer. The thickness of the BCPD layers is 2 µm. This value for the thickness was chosen to obtain a significant phase difference between the BCPD film in the neutral form and BCPD in the charged form.

With the BCPD-based electrochromic layers in the Mach–Zehnder interferometer, the probability of measuring the photon using detector D1 is provided in Equation 19:

${D1}_{BCPD}={\left|x^{T}BSMBCPDBSx\right|}^2.(19)$

Meanwhile, the probability of measuring the photon using detector D2 is obtained as shown in Equation 20:

${D2}_{BCPD}={\left|y^{T}BSMBCPDBSx\right|}^2.(20)$

Here, $x^T$ is the transpose of vector $x$, while $y^T$ is the transpose of vector $y$.

In the waveguide-based Mach–Zehnder interferometer, the phase shift for two branches of equal length L is provided in Equation 21 (De Ridder et al., 2006):

$\Delta \phi ={\beta }_{2}L-{\beta }_{1}L.(21)$

The parameters ${\beta }_1$ and ${\beta }_2$ in Equation 21 can be defined in the following way (Equation 22):

${\beta }_{1,2}=\frac{2\pi n_{1,2}}{\lambda }.(22)$

In Equation 22, $n_{1,2}$ is the complex refractive index for the materials in the two branches of the Mach–Zehnder interferometer. For this study, $n_1=n_n^{complex}$ and $n_2=n_c^{complex}$. The output powers at ports 1 and 2 (Figure 2b) are obtained using Equations 23 and 24, respectively:

$P_{out,1}={\mathit{sin}}^2\left(\frac{\Delta \phi }{2}\right)P_{in,1}+{\mathit{cos}}^2\left(\frac{\Delta \phi }{2}\right)P_{in,2},(23)$

$P_{out,2}={\mathit{cos}}^2\left(\frac{\Delta \phi }{2}\right)P_{in,1}+{\mathit{sin}}^2\left(\frac{\Delta \phi }{2}\right)P_{in,2}.(24)$

Figure 2

Figure 2. (a) Sketch of the Mach–Zehnder interferometer with direct transmission of light, showing the photon source, the two mirrors M, the two 50/50 beam-splitters BS1 and BS2, the two electrochromic layers BCPDx (along the x-direction transmitted by BS1) and BCPDy (along the y-direction reflected by BS1), and the two detectors D1 and D2. (b) Waveguide-based Mach–Zehnder interferometer with inputs In1 and In2, two 3-dB couplers, and two outputs Out1 and Out2. In one of the branches of the interferometer, the electric field induces a change in the refractive index of BCPD.

Results and discussion

Figure 1 shows the chemical sketch of the molecule BCPD. The carbon (black), oxygen (red), and nitrogen (blue) atoms are indexed as the Hückel Hamiltonian.

In the Hamiltonian of BCPD reported in Equation 25, the diagonal elements correspond to the different atoms of the molecules, while the non-null off-diagonal elements correspond to nearest neighbor elements:

$H=\left[\begin{array}{c}\begin{array}{c}\begin{array}{cccccccccccccccc}\alpha N& \beta N& 0& 0& \beta N& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0\end{array}\\ \begin{array}{cccccccccccccccc}\beta N& \alpha & \beta & 0& 0& 0& 0& 0& 0& 0& 0& 0& \beta O& 0& 0& 0\end{array}\\ \begin{array}{cccccccccccccccc}0& \beta & \alpha & \beta & 0& \beta & 0& 0& 0& 0& 0& 0& 0& 0& 0& 0\end{array}\\ \begin{array}{cccccccccccccccc}0& 0& \beta & \alpha & \beta & 0& 0& 0& \beta & 0& 0& 0& 0& 0& 0& 0\end{array}\\ \begin{array}{c}\begin{array}{cccccccccccccccc}\beta N& 0& 0& \beta & \alpha & 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& \beta O\end{array}\\ \begin{array}{cccccccccccccccc}0& 0& \beta & 0& 0& \alpha & \beta & 0& 0& 0& 0& 0& 0& 0& 0& 0\end{array}\\ \begin{array}{cccccccccccccccc}0& 0& 0& 0& 0& \beta & \alpha & \beta & 0& 0& 0& \beta & 0& 0& 0& 0\end{array}\\ \begin{array}{cccccccccccccccc}0& 0& 0& 0& 0& 0& \beta & \alpha & \beta & \beta & 0& 0& 0& 0& 0& 0\end{array}\end{array}\end{array}\\ \begin{array}{c}\begin{array}{cccccccccccccccc}0& 0& 0& \beta & 0& 0& 0& \beta & \alpha & 0& 0& 0& 0& 0& 0& 0\end{array}\\ \begin{array}{cccccccccccccccc}0& 0& 0& 0& 0& 0& 0& \beta & 0& \alpha & \beta N& 0& 0& 0& \beta O& 0\end{array}\\ \begin{array}{cccccccccccccccc}0& 0& 0& 0& 0& 0& 0& 0& 0& \beta N& \alpha N& \beta N& 0& 0& 0& 0\end{array}\\ \begin{array}{cccccccccccccccc}0& 0& 0& 0& 0& 0& \beta & 0& 0& 0& \beta N& \alpha & 0& \beta O& 0& 0\end{array}\\ \begin{array}{c}\begin{array}{cccccccccccccccc}0& \beta O& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& \alpha O& 0& 0& 0\end{array}\\ \begin{array}{cccccccccccccccc}0& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& \beta O& 0& \alpha O& 0& 0\end{array}\\ \begin{array}{cccccccccccccccc}0& 0& 0& 0& 0& 0& 0& 0& 0& \beta O& 0& 0& 0& 0& \alpha O& 0\end{array}\\ \begin{array}{cccccccccccccccc}0& 0& 0& 0& \beta O& 0& 0& 0& 0& 0& 0& 0& 0& 0& 0& \alpha O\end{array}\end{array}\end{array}\end{array}\right].(25)$

Here, the on-site energy α is set to 0 eV, and the nearest neighbor element β is set to −3 eV, following the publication of Solomon et al. (2011). Equation 26 reports the other matrix elements:

$\left\{\begin{array}{c}\begin{array}{c}\alpha N=0.51\beta \\ \alpha O=0.97\beta \end{array}\\ \begin{array}{c}\beta N=1.02\beta \\ \beta O=1.02\beta \end{array}\end{array}\right..(26)$

In the literature, the absorption spectrum of BCPD has been studied (Abdinejad et al., 2018). Neutral BCPD shows two absorption peaks at approximately 320 and 395 nm, while charged BCPD shows an additional peak at 605 nm. It is worth noting that the spectrum of charged BCPD is related to photocoloration, i.e., exposure to light at 365 nm. However, the authors show that the effects of photocoloration are very similar to the effects of electrocoloration in terms of the optical response of BCPD (Abdinejad et al., 2018). In Supplementary Figure S1 of Supplementary Material, the calculated light absorption spectra via the Hückel theory and via DFT have been reported for a comparison. Both the calculations are in fair agreement with the experimental data, considering that both predict the formation of the peak in red upon electrocoloration. Supplementary Material contains the transitions with the corresponding oscillator strengths and the geometry of BCPD related to the DFT calculations. Finally, in Supplementary Material, we list the weight of the individual excitations for each transition with which we have built the absorption spectrum employing density functional theory. In this work, although all single-particle transitions contribute in principle to the absorption spectrum, only a selected subset has been considered (i.e., the transitions between the HOMO and the six lowest unoccupied orbitals).

We studied the complex refractive index of BCPD in the neutral form and with the addition of a charge. We focused on a wavelength of 800 nm since commercially available single-photon sources are present in this wavelength.

The extinction coefficient k for the neutral BCPD is $k_n=1.1788\times {10}^{-23}$, while the extinction coefficient k for the charged BCPD is $k_c=1.1078\times {10}^{-16}$. The real part of the refractive index for the neutral BCPD is n_n = 1.5773, while the real part of the refractive index for the charged BCPD is n_c = 1.4901.

We have simulated the Mach–Zehnder interferometer with direct transmission of light, taking into account the 2-μm-thick slabs of BCPD. Examples of electrochromic organic films are reported in the literature (Rendón-Enríquez et al., 2023). Given that the possible formation of different types of molecular aggregates in thin films is an important aspect (Zhang et al., 2020), we would like to emphasize that in this work, we have not studied in depth the possibility of having BCPD aggregates in this slab, but this could be an interesting topic for future research. With BCPD in the neutral form along the x-axis and BCPD in the charged form along the y-axis (Figure 2a), the probability of detecting a photon at detector D1 is 0.8056, while the probability of detecting a photon at detector D2 is 0.1944.

We then simulated the waveguide-based Mach–Zehnder interferometer in a glass slab, with 2-μm-long arms. Considering the input power at port 2 equal to 0 ($P_{in,2}=0$, Figure 2b) and considering branches of the interferometer with length L = 2 μm, the probability of detecting a photon at output port 1 is 0.4001, while the probability of detecting a photon at output port 2 is 0.5999.

Finally, with the non-equilibrium Green’s function formalism, we determined the conductance of BCPD. In Figure 3, we show the elastic transmission spectrum (in the logarithm scale) of BCPD. A node close to −3 eV is evident due to destructive interference (Solomon et al., 2011). Other nodes occurring in the spectrum are higher-order nodes, i.e., supernodes, at −2 eV and 2 eV, which can be ascribed to cross-conjugation in BCPD—the delocalization of π electrons branches into two or more non-continuous paths (Solomon et al., 2011; Zhao et al., 2025). When setting the Fermi energy at 0 eV in Equation 13, the conductance of BCPD is 0.127 µS.

Figure 3

Figure 3. Simulated elastic transmission spectrum as a function of the energy of N,N′-bis(cysteine)pyromellitic diimide.

Conclusion

In this work, we studied an electrochromic molecule, BCPD. Employing a simple method such as the Hückel theory, we were able to investigate the optical properties of BCPD in the neutral and charged forms, with good agreement with the experimental findings. The electrochromism of BCPD is related to the refractive index of the molecule, which depends on the applied electric field. Taking advantage of this refractive index control, we simulated the probability of single-photon detection on a Mach–Zehnder interferometer with direct light transmission and a waveguide-based Mach–Zehnder interferometer—for example, buried in a glass slab. The waveguide-based Mach–Zehnder interferometer can be useful in integrated photonic platforms. We also studied the conductance of BCPD with the non-equilibrium Green’s function formalism. The results of this study could be of interest for quantum communication.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author contributions

FS: Conceptualization, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review and editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This project received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement No. 816313) and the EIC LEAF (grant agreement No. 101186701).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphot.2025.1685128/full#supplementary-material

References

Abdinejad, T., Zamanloo, M. R., Alizadeh, T., and Mahmoodi, N. O. (2017). Dual photo-electrochromic diimides derived from aliphatic aminothiols and π-electron deficient aromatic dianhydrides. Dyes Pigments 146, 203–209. doi:10.1016/j.dyepig.2017.06.063

CrossRef Full Text | Google Scholar

Abdinejad, T., Zamanloo, M. R., Alizadeh, T., Mahmoodi, N. O., and Pouran, S. R. (2018). Photochromic and electrochromic diimide synthesized simply from inexpensive compounds: a multidisciplinary experiment for undergraduate students. J. Chem. Educ. 95, 1642–1647. doi:10.1021/acs.jchemed.7b00540

CrossRef Full Text | Google Scholar

Ambrose, B., Krishnan, M., Ramamurthy, K., and Kathiresan, M. (2024). Recent advances in molecular engineering for viologen-based electrochromic materials: a mini-review. Adv. Photonics Res. 5, 2400016. doi:10.1002/adpr.202400016

CrossRef Full Text | Google Scholar

Bird, C. L., and Kuhn, A. T. (1981). Electrochemistry of the viologens. Chem. Soc. Rev. 10, 49–82. doi:10.1039/CS9811000049

CrossRef Full Text | Google Scholar

Chang, W. B., Mai, C.-K., Kotiuga, M., Neaton, J. B., Bazan, G. C., and Segalman, R. A. (2014). Controlling the thermoelectric properties of thiophene-derived single-molecule junctions. Chem. Mater 26, 7229–7235. doi:10.1021/cm504254n

CrossRef Full Text | Google Scholar

Chapkin, K. D., Bursi, L., Clark, B. D., Wu, G., Lauchner, A., Tsai, A.-L., et al. (2020). Effects of electronic structure on molecular plasmon dynamics. J. Phys. Chem. C 124, 20450–20457. doi:10.1021/acs.jpcc.0c06072

CrossRef Full Text | Google Scholar

Chen, D., Chua, M. H., He, Q., Zhu, Q., Wang, X., Meng, H., et al. (2025). Multifunctional electrochromic materials and devices recent advances and future potential. Chem. Eng. J. 503, 157820. doi:10.1016/j.cej.2024.157820

CrossRef Full Text | Google Scholar

Cho, H. M., Hwang, Y. J., Oh, H. S., Hwang, H. S., Kim, K., Kim, S.-U., et al. (2025). Recent advances in electrochromic devices: from multicolor to flexible applications. Adv. Photonics Res. 6, 2400103. doi:10.1002/adpr.202400103

CrossRef Full Text | Google Scholar

De Ridder, R. M., and Roeloffzen, C. G. H. (2006). “Interleavers,” in Wavelength filters in fibre optics springer series in optical sciences. Editor H. Venghaus (Berlin, Heidelberg: Springer Berlin Heidelberg), 381–432. doi:10.1007/3-540-31770-8_10

CrossRef Full Text | Google Scholar

Dresselhaus, M. S. (2017). Solid state physics part II. Available online at: https://web.mit.edu/course/6/6.732/www/6.732-pt2.pdf.

Google Scholar

Eichkorn, K., Weigend, F., Treutler, O., and Ahlrichs, R. (1997). Auxiliary basis sets for main row atoms and transition metals and their use to approximate Coulomb potentials. Theor. Chem. Acc. Theory Comput. Model Theor. Chim. Acta 97, 119–124. doi:10.1007/s002140050244

CrossRef Full Text | Google Scholar

Guo, X., Facchetti, A., and Marks, T. J. (2014). Imide- and amide-functionalized polymer semiconductors. Chem. Rev. 114, 8943–9021. doi:10.1021/cr500225d

PubMed Abstract | CrossRef Full Text | Google Scholar

Guy, G. R., and Troy, V. V. (2022). Physical chemistry. MIT OpenCourseWare. Available online at: https://ocw.mit.edu/courses/chemistry/5-61-physical-chemistry-fall-2007/ (Accessed February 5, 2022).

Google Scholar

Hanwell, M. D., Curtis, D. E., Lonie, D. C., Vandermeersch, T., Zurek, E., and Hutchison, G. R. (2012). Avogadro: an advanced semantic chemical editor, visualization, and analysis platform. J. Cheminformatics 4, 17. doi:10.1186/1758-2946-4-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Huckel Theory and HuLiS (2013). Huckel theory and HuLiS: a calculator that also describes mesomerism. Available online at: http://www.hulis.free.fr/ (Accessed June 8, 2020).

Google Scholar

Janoš, J., Slavíček, P., and Curchod, B. F. E. (2025). Selecting initial conditions for trajectory-based nonadiabatic simulations. Acc. Chem. Res. 58, 261–270. doi:10.1021/acs.accounts.4c00687

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, Z., Doorn, S. K., and Sykora, M. (2015). Electrochromic graphene molecules. ACS Nano 9, 4043–4049. doi:10.1021/acsnano.5b00093

PubMed Abstract | CrossRef Full Text | Google Scholar

Kohandani, R., and Saini, S. S. (2022). Extracting optical absorption characteristics from semiconductor nanowire arrays. Nanotechnology 33, 395204. doi:10.1088/1361-6528/ac74cc

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, C., Yang, W., and Parr, R. G. (1988). Development of the colle-salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 37, 785–789. doi:10.1103/PhysRevB.37.785

PubMed Abstract | CrossRef Full Text | Google Scholar

Lehtola, S., Steigemann, C., Oliveira, M. J. T., and Marques, M. A. L. (2018). Recent developments in libxc — a comprehensive library of functionals for density functional theory. SoftwareX 7, 1–5. doi:10.1016/j.softx.2017.11.002

CrossRef Full Text | Google Scholar

Marques, M. A. L., Oliveira, M. J. T., and Burnus, T. (2012). Libxc: a library of exchange and correlation functionals for density functional theory. Comput. Phys. Commun. 183, 2272–2281. doi:10.1016/j.cpc.2012.05.007

CrossRef Full Text | Google Scholar

Neese, F. (2012). The ORCA program system. WIREs Comput. Mol. Sci. 2, 73–78. doi:10.1002/wcms.81

CrossRef Full Text | Google Scholar

Ng, J. C., Li, C., Herman, P. R., and Qian, L. (2012). Femtosecond laser writing of a flat-top interleaver via cascaded Mach-Zehnder interferometers. Opt. Express 20, 17894–17903. doi:10.1364/OE.20.017894

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogilvie, J. F., and Fee, G. J. (2013). Equivalence of Kramers-Kronig and fourier transforms to convert between optical dispersion and optical spectra. MATCH – Commun. Math. Comput. Chem. 69, 249–262.

Google Scholar

Otaegui, J. R., Mena, S., Dorsainvil, J. M., Guirado, G., Ruiz-Molina, D., Hernando, J., et al. (2025). Enhanced electrochromic smart windows based on supramolecular viologen tweezers. Chem. Mater 37, 2220–2229. doi:10.1021/acs.chemmater.4c03174

PubMed Abstract | CrossRef Full Text | Google Scholar

Pankove, J. I. (1975). Optical processes in semiconductors. Unabridged republication, with slight corrections. New York, NY: Dover Publications, Inc.

Google Scholar

Pradhan, L., Kumar Das, T., and Kumar Jena, B. (2025). Advances in next-generation electrochromic supercapacitors: emerging applications, challenges and future perspectives. Sustain Energy Fuels 9, 1917–1956. doi:10.1039/D4SE01551G

CrossRef Full Text | Google Scholar

Rendón-Enríquez, I., Palma-Cando, A., Körber, F., Niebisch, F., Forster, M., Tausch, M. W., et al. (2023). Thin polymer films by oxidative or reductive electropolymerization and their application in electrochromic windows and thin-film sensors. Molecules 28, 883. doi:10.3390/molecules28020883

PubMed Abstract | CrossRef Full Text | Google Scholar

Rioux, F. (2019). 8.86: simulating a quantum computer with a Mach-Zehnder interferometer. Chem. Libr. Available online at: https://chem.libretexts.org/Bookshelves/Physical_and_Theoretical_Chemistry_Textbook_Maps/Quantum_Tutorials_(Rioux)/08%3A_Quantum_Teleportation/8.86%3A_Simulating_a_Quantum_Computer_with_a_Mach-Zehnder_Interferometer (Accessed June 12, 2025).

Google Scholar

Schäfer, A., Horn, H., and Ahlrichs, R. (1992). Fully optimized contracted Gaussian basis sets for atoms Li to Kr. J. Chem. Phys. 97, 2571–2577. doi:10.1063/1.463096

CrossRef Full Text | Google Scholar

Schäfer, A., Huber, C., and Ahlrichs, R. (1994). Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr. J. Chem. Phys. 100, 5829–5835. doi:10.1063/1.467146

CrossRef Full Text | Google Scholar

Scotognella, F. (2020). Multilayer plasmonic photonic structures embedding photochromic molecules or optical gain molecules. Phys. E Low-Dimens Syst. Nanostructures 120, 114081. doi:10.1016/j.physe.2020.114081

CrossRef Full Text | Google Scholar

Scotognella, F. (2024a). Quantum computation with photochromic films in a Mach–Zehnder interferometer. Front. Nanotechnol. 6. doi:10.3389/fnano.2024.1422573

CrossRef Full Text | Google Scholar

Scotognella, F. (2024b). Control of the phase shift in a Mach-Zehnder interferometer with temperature via VO_2 phase transition and with light via pDTE photochromic switching. Available online at: https://zenodo.org/records/12172152 (Accessed May 23, 2025).

Google Scholar

Solomon, G. C., Bergfield, J. P., Stafford, C. A., and Ratner, M. A. (2011). When “small” terms matter: coupled interference features in the transport properties of cross-conjugated molecules. Beilstein J. Nanotechnol. 2, 862–871. doi:10.3762/bjnano.2.95

PubMed Abstract | CrossRef Full Text | Google Scholar

Stec, G. J., Lauchner, A., Cui, Y., Nordlander, P., and Halas, N. J. (2017). Multicolor electrochromic devices based on molecular plasmonics. ACS Nano 11, 3254–3261. doi:10.1021/acsnano.7b00364

PubMed Abstract | CrossRef Full Text | Google Scholar

Stolar, M. (2020). Organic electrochromic molecules: synthesis, properties, applications and impact. Pure Appl. Chem. 92, 717–731. doi:10.1515/pac-2018-1208

CrossRef Full Text | Google Scholar

Valeev, E. F. (2014). A library for the evaluation of molecular integrals of many-body operators over Gaussian functions. Available online at: http://libint.valeyev.net/.

Google Scholar

Wiebeler, C., Bader, C. A., Meier, C., and Schumacher, S. (2014). Optical spectrum, perceived color, refractive index, and non-adiabatic dynamics of the photochromic diarylethene CMTE. Phys. Chem. Chem. Phys. 16, 14531–14538. doi:10.1039/C3CP55490B

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Zhao, Z., Turley, A. T., Wang, L., McGonigal, P. R., Tu, Y., et al. (2020). Aggregate science: from structures to properties. Adv. Mater 32, 2001457. doi:10.1002/adma.202001457

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, N., Jeon, S. J., and Li, Y. (2025). Cross-conjugated polymer semiconductors. Macromol. Rapid Commun. 46, e00281. doi:10.1002/marc.202500281

PubMed Abstract | CrossRef Full Text | Google Scholar