Abstract
In this work we quantitatively study the reliability of the frozen nuclei approximation for ultrafast dynamics. Specifically we study laser excitation of HCCI+ from its ground state to the first electronically excited state. The population of the first excited state is obtained by both the frozen nuclei approximation and by multidimensional nuclear dynamics. Detailed comparison of the results by the two methods are performed to provide quantitative criteria for the reliability of the frozen nuclei approximation for this system.
1 Introduction
The rapid advances of ultrafast science and technology have made it possible to manipulate electron dynamics in molecular systems with ultrashort laser pulses. In particular, laser induced electron density redistribution such as charge transfer (Marcus, 1956; May and Kühn, 2008) and charge migration (Weinkauf et al., 1996, 1997; Cederbaum and Zobeley, 1999; Calegari et al., 2014; Kraus et al., 2015) have been extensively investigated. In general, charge migration prefers frozen nuclei or small amplitude nuclear motions, while charge transfer is typically accompanied by large amplitude nuclear motions. The research of charge transfer processes has a relatively long history. While ultrafast charge migration emerged as a hot topic during the past two decades (Weinkauf et al., 1996, 1997; Remacle et al., 1998; Cederbaum and Zobeley, 1999; Remacle and Levine, 1999; Barth and Manz, 2006; Kanno et al., 2006; Yudin et al., 2006; Remacle et al., 2007; Kanno et al., 2010; Mineo et al., 2012; Calegari et al., 2014; Kraus et al., 2015; Li et al., 2015; Yamaki et al., 2016; Wörner et al., 2017; Mineo et al., 2021). It should be noted that the first attosecond charge migration phenomenon was already introduced in 1944 (Eyring et al., 1944) and was largely forgotten during the next decades. Surveys of the literature on ultrafast charge migration can be found in Ref. (Jia et al., 2017a; Wörner et al., 2017). Below we summarize some typical features of ultrafast charge migration and its connection to the frozen nuclei approximation (FNA).
Ultrafast charge migration typically represents quantum dynamics of a coherent superposition of more than one electronic state. The typical time scale of ultrafast charge migration ranges from several hundred attoseconds to a few femtoseconds which makes the experimental observation (Kraus et al., 2015) rather difficult. For such a short time, the frozen nuclei approximation has been widely used for theoretical work of ultrafast charge migration. There are also several theoretical investigations which include the effects of nuclear motions (Bandrauk et al., 2009; Kanno et al., 2010; Ulusoy and Nest, 2012; Mendive-Tapia et al., 2013; Mineo et al., 2014; Despré et al., 2015; Mineo et al., 2021). The amplitude of charge migration can be significantly modulated by nuclear motions, in particular for relatively long-time dynamics (Mendive-Tapia et al., 2013; Mineo et al., 2014; Jia et al., 2019a,b). In general, the FNA is widely believed to be only valid for short time pulses, but there are no quantitative criteria for how short the pulses should be. This serves as the motivation for the present work: to seek quantitative criteria for the reliability of the FNA. Specifically, we will investigate short-pulse excitations of HCCI+ by systematically varying the laser parameters in a sufficiently wide region.
The choice of HCCI+ as our model of interest is based on the availability of experimental data (Heilbronner et al., 1971; Kraus et al., 2015) and theoretical techniques (Jia et al., 2019b). The combined experimental and theoretical reconstruction of attosecond charge migration has been reported for ultrafast ionization of HCCI (Kraus et al., 2015). Coherent superposition of the ground and first excited states has been created and analyzed. Subsequent theoretical investigations of ultrafast charge migration in HCCI+ (Jia et al., 2017b; Ding et al., 2017) related to the experimental observation (Kraus et al., 2015) exploit the FNA. In-depth investigations of simulations and manipulations of charge migration in HCCI+ including multidimensional nuclear dynamics have been reported recently (Jia et al., 2019a,b). However, no comparisons between the results of multidimensional nuclear dynamics and the ones of the FNA are available.
In the present work, we will investigate the reliability of the FNA by comparing the FNA and multidimensional nuclear dynamics. The remainder parts of the paper are organized as follows. Section 2 contains the model and methods for numerical calculations. Section 3 presents the results and discussion. The conclusions are drawn in Section 4.
2 Model and Methods
We focus on laser excitations of HCCI+ from its ground state. Full dimensional simulations of the system involve sets of electronic coordinates r = {r1, r2, … } and nuclear coordinates R = {R1, R2, … }. Here ri and Rj are the spatial coordinates of the i-th electron and the j-th nucleus, respectively. It is convenient to use the Dirac notation for the electronic degrees of freedom. The total wavefunction of the system is thusUsing the Born-Huang expansion (Born and Oppenheimer, 1927; Born and Huang, 1954), the total wavefunction can be expressed in terms of the electronic eigenstates |k(R)⟩ which are the solutions of the standard time-independent electronic schrödinger equationThe corresponding electronic eigenenergy Vk(R) is the k-th potential energy surface (PES). Here Vk(R) and ⟨r|k(R)⟩ are the same as the ones used in Refs. (Jia et al., 2019a,b) which are calculated by Molpro (Werner et al., 2012) using the state-averaged CASSCF(15,13) with cc-pVQZ basis set (cc-pVQZ-pp for iodine).
According to (Kraus et al., 2015; Jia et al., 2019a,b), nonadiabatic couplings between different electronic states |k(R)⟩ are negligible. The total Hamiltonian for HCCI+ in an external laser field E(t) can be approximated aswhere T(R) is the nuclear kinetic energy and μkk′(R) = ⟨k(R)|μ|k′(R)⟩ is the transition (or permanent) dipole moment. The laser pulse has a Gaussian shape with maximum amplitude Emax and carrier frequency ωwhere ez is the direction of the electric field. For convenience, the electric field and the molecules are oriented along the z-axis. In the literature, there are different choices of the parameter a in Eq. 4. For the present work we setfor easier characterization of the pulse duration. We define the pulse duration as the full width at half maximum (FWHM) of s(t), which is just T in Eq. 4.
The quantum dynamics of the system can be simulated by the time-dependent schrödinger equation subject to initial condition at t = −∞where χg,v=0(R) is the vibrational ground state wavefunction of the lowest potential energy surface Vg(R). For convenience we use k = g, e to represent the lowest and first excited electronic states, respectively. The wave packet is numerically propagated by means of the split operator method (Leforestier et al., 1991).
Subsequently, we can obtain the population of the electronic state |k(R)⟩ according towhere χk (R,t) = ⟨k(R)|Ψ(R, t)⟩ is the nuclear wave packet on the k-th PES Vk(R). It contains seven vibrational coordinates. According to Ref. (Jia et al., 2019a,b), one-dimensional (1D), three-dimensional (3D), and seven-dimensional (7D) calculations lead to essentially the same results. In the 3D calculations, the H-C, C-C and C-I bond lengths are explicitly taken into account and the four bending degrees of freedom are neglected. This kind of approximation is reasonable for linear molecules, such as HCCI+. In the present work we use the same 3D calculations for the nuclear wave packet χk (R,t) as in Ref. (Jia et al., 2019b). Then we mainly focus on the population of the first electronically excited state Pk(t) for k = e.
To check the reliability of the FNA, we further calculated the population of the first electronically excited state using the FNA. Accordingly, the molecular structure is fixed at the minimum of the lowest PES Vg(R). This structure is called equilibrium structure Req. The corresponding transition dipole moment is μeq ≡μge (R = Req). The electronic wavefunction of the FNA is expanded asthe time dependent coefficient ck(t) can be obtained subject to the initial condition ck(t = −∞) = δkg. The corresponding population is
Throughout this work we fix the carrier frequency of the laser in Eq. 4 as ℏω = Ve (Req) − Vg (Req). We only focus on the final population at t = tf when the laser pulse is off. This leads to the following analytical expression (Jia et al., 2017a)where denotes the z component of the transition dipole at Req. The relative error of the FNA with respect to multi-dimensional nuclear dynamic is defined asFor all the subsequent numerical calculations we set tf = 5T. However, this should not be wrongly interpreted as the FNA is valid even for t = 5T. We choose tf = 5T just to make use of the property that the results presented in this work do not depend on different choices of tf as long as tf ≥ T.
3 Results and Discussion
The equilibrium structure Req of HCCI+ is linear with bond lengths RHC = 1.06 Å, RCC = 1.21 Å and RCI = 1.95 Å. The corresponding vertical excitation energy from ground state |g(Req)⟩ to the first excited state |e(Req)⟩ is ℏω = 2.41 eV. For typical pulse durations, there are sufficient numbers of cycles in E(t) to make the electronic transition resonant. The corresponding transition dipole has only a z-component, which is 3.21 Debye. Subsequently we calculated the population of the first electronically excited state according to Eq. 10 by the FNA. Convergence tests are performed for nuclear dynamics simulations such that the corresponding population of the first electronically excited state Pe(t) does not change subject to further increase of the grid-region or decrease of the spatial or temporal steps. We first analyze the dependence of the results on the pulse durations with the other parameters fixed. Specifically, the maximum amplitude of the electric field Emax is fixed at 2.0 × 109 V/m. The detailed comparison between Pe(t = tf) and is shown in Figure 1 for T ≤ 20 fs.
FIGURE 1
As can be seen from Figure 1A, the deviation between and Pe(tf) gradually increases with the pulse duration T, for the region of short pulses. For T ≤ 20 fs, Pe(t) keeps increasing with T. However, first increases and then decreases for T ≤ 20 fs. This kind of qualitative deviation will be further discussed below. According to Eq. 10 will oscillate periodically with the pulse duration T. In Figure 1A, reaches its maximum for T = 15.45 fs. However, Pe(tf) is still substantially below one even for T = 20 fs.
To quantitatively compare and Pe(tf), the relative error defined in Eq. 11 is shown in Figure 1B. The relative error increases relatively slowly when the pulse duration T is smaller than 5 fs, and increases rapidly when T is larger than 5 fs. For long pulses, say T ≥ 15 fs, Figure 1B shows significant decrease of . However, this is pure coincidence. As can be identified from Figure 1A, the trend of is already qualitatively wrong for T > 15.45 fs. Smaller relative error in this region does not imply better agreement between the frozen nuclei approximation and real physics. Consequently, we focus on short pulses for which the FNA is expected to be reasonable. Accordingly, we add three vertical lines in Figure 1B for relative errors of 40, 20, and 10%, respectively. The corresponding pulse durations with fixed value of Emax = 2.0×109 V/m are T = 6.09, 4.65, and 1.97 fs, respectively.
Figure 2 shows the color-coded contour plots for the dependence of and Pe (tf) on the amplitude of the electric field Emax and the pulse duration T. The full set of the involved parameters span the region 0.5 GV/m ≤ Emax ≤ 4.0 GV/m and T ≤ 120 fs. The region of Emax more or less covers the reported amplitudes of lasers exploited in typical applications in the literature. The region of T reaches the first revival of charge migration in HCCI+ reported in Ref. (Jia et al., 2019b). As can be seen from Figure 2A, oscillates between 1 and 0 periodically with Emax or T. Larger values of Emax or T corresponds to smaller oscillation period of , c.f., Eq. 10. Figure 2B is the same as Figure 2A except that is replaced by Pe(tf) which is obtained by performing 3D nuclear dynamics simulations. We can immediately identify that Figures 2A,B are qualitatively different for relatively long laser lulses (e.g., for T ≥ 20 fs). We therefore consider the results for the FNA are not meaningful for relatively long laser pulses. This is quite natural. Even if is coincidentally very close to Pe(tf) for T ≥ 20 fs, the nuclear wave packet is quite different from the initial one as has been reported in Ref. (Jia et al., 2019b). Accordingly, the basic assumption of the FNA breaks down for T ≥ 20 fs. We further plot Figure 2C which is the same with Figure 2B except that the population of the first electronically excited state is obtained by performing 1D nuclear dynamics simulations which explicitly treats the C-I stretch (Jia et al., 2019a). The results of 1D and 3D simulations agree well with each other, which confirms the findings in Refs. (Jia et al., 2019a,b).
FIGURE 2
To systematically study the reliability of the FNA for short pulses, the difference between in Figure 2A and Pe (tf) in Figure 2B is shown in Figure 3A, for the region T ≤ 20 fs when the FNA may be expected to work. The dependence of on Emax and T is shown by color-code contour plots. Results for which are not resolved, since results with errors larger than 50 percent are in general not helpful. For any fixed value of Emax, the deviation between and Pe (tf) first increases then decreases with T after reaching a maximum. More complicated features can be found for relatively large Emax combined with relatively long T. However, as discussed above for Figure 1, the results of in this complicated region in Figure 3A (specifically, after reaching maximum deviations) can agree better with Pe(tf) by coincidence. In the following we only focus on the left bottom region of Figure 3A before the corresponding deviation reaches its maximum for any fixed value of Emax. As can be seen from Figure 3A, the deviations are rather small for T ≤ 5 fs for all the different values of Emax involved in the present work. This implies that we may roughly use T ≤ 5 fs as the criterion for the reliability of the FNA. The deviation not only increases with T for any given Emax, but also increases with Emax for any given T. According to Eq. 10, increases with the product of Emax and T before reaching its maximum. From Figure 3A we can also find that the deviation essentially increases with the product of Emax and T before reaching its maximum. Better criteria for the reliability of the FNA can be obtained by analyzing the relative error of with respect to Pe (tf) shown in Figure 3B. In the same spirit, we only need to focus on the left region of Figure 3B indicated by the dashed curve in which the relative error never exceeds 50%. For relatively long pulses, say T = 7 fs, the relative error of the FNA is already larger than 50% for any value of Emax shown in Figure 3B. In this case, the FNA is no longer reliable for T ≥ 7 fs.
FIGURE 3
In the short pulse region, the relative error increases with the pulse duration T. To give quantitative criteria for the reliability of the FNA, we define certain characteristic pulse durations as follows:For example, if we want to use the FNA to obtain results with relative errors smaller than 5%, we need to set the pulse durations of the lasers to be smaller than T95. Similarly for the meanings of T90, T80, and T60. According to our model the characteristic pulse durations T95, T90, T80, and T60 only depend on one parameter Emax, which will be investigated subsequently.
The detailed dependence of T95, T90, T80, and T60 on the maximum amplitude of the electric field Emax is shown in Figure 4. A quite good property for the characteristic pulse durations is that T95, T90, and T80 almost do not depend on Emax. The corresponding values are T95 = 1.32 fs, T90 = 1.97 fs, and T80 = 4.65 fs respectively. The value of T60 increases with Emax extremely slightly from 6.09 to 6.14 fs. For relatively high standard criteria, say relative errors below 20%, the corresponding characteristic pulse durations are quite robust with respect to different amplitudes of lasers. This greatly simplifies the criteria for choosing proper lasers for applications of short pulse excitations of HCCI+. Essentially, we only need to care about the durations of the laser pulses with quantitative guidance derived from Figure 4 for the reliability of the FNA.
FIGURE 4
The FNA only considers the electronic degrees of freedom and neglects the nuclear motions. Mathematically this corresponds to a large overlap of the time-dependent and the initial nuclear wave packets. The overlap can be estimated as the product of the corresponding overlap for each normal mode. The overlap for a normal mode may be approximated as for short time dynamics. Here f(t) is the displacement of the normal coordinate with respect to its initial value in units of the standard deviation of the initial wave packet for this normal mode. Typically, there will be only one or a few modes with f2(t) substantially above zero, which are called active modes. The overall overlap is thus mainly determined by the active modes. For the present case, there is only one active mode which is the C-I stretch with period 86 fs (Jia et al., 2019b). Due to the relatively large amplitude of the C-I stretch mode, the function quickly decreases to zero (Jia et al., 2019a,b). In this case, the duration of the pulse must be much shorter than a vibrational period to keep the overlap large enough. According to the results of Figure 4, all the characteristic pulse durations of HCCI+ are smaller than of a vibrational period. However, for some molecules with sufficiently small vibrational amplitudes for all the active modes, the overlap can be relatively large for a rather long time. For such cases, the effects of nuclear dynamics can be neglected for a longer time than just a few femtoseconds (Kanno et al., 2010; Ulusoy and Nest, 2012; Despré et al., 2015).
4 Conclusion
We have systematically investigated the population of the first electronically excited state of HCCI+ excited by different laser pulses. The amplitudes and durations of the laser pulses span a rather large domain for typical applications. The deviations between the results obtained by the frozen nuclei approximation and the ones obtained by multidimensional nuclear dynamics are calculated and analyzed in detail to check the reliability of the FNA. As expected the validity of the FNA can be admitted for sufficiently short laser pulses. Quantitative criteria for the reliability of the FNA are obtained. Specifically if we want to limit the relative errors of the FNA within 5% (or 10, or 20, or 40%), the durations of the laser pulses should be less than T95 = 1.3 fs (or T90 = 2.0 fs, or T80 = 4.7 fs, or T60 = 6.1 fs). For example, ultrafast charge migration in HCCI+ is reconstructed in Ref. (Kraus et al., 2015). for the first period of 1.85 fs. By extrapolation of our results, the error of the reported charge migration in HCCI+ for the first period is less than 10%. For short pulses with durations up to T60, the relative errors of the FNA are found to be almost independent of the amplitudes of the laser pulses. The results of the present work are expected to provide valuable guidance to future investigations of short pulse excitations of HCCI+.
Statements
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
YY proposed the conception and design of the study. DJ carried out all the quantum chemical calculations and the quantum dynamics simulations, and prepared all Figures. YY wrote the zero-order draft. All the authors contributed to the submitted version of the manuscript.
Funding
This work was supported by the National Key Research and Development Program of China (2017YFA0304203), the Program for Changjiang Scholars and Innovative Research Team (IR_17R70), the National Natural Science Foundation of China (11904215), the 111 project (Grant No. D18001), the Fund for “Shanxi 1331 Project,” and the Hundred Talent Program of Shanxi Province.
Acknowledgments
We would like to express our gratitude to Professor Jörn Manz (Berlin) for stimulating discussions and careful reading of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
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.
References
1
BandraukA. D.ChelkowskiS.CorkumP. B.ManzJ.YudinG. L. (2009). Attosecond Photoionization of a Coherent Superposition of Bound and Dissociative Molecular States: Effect of Nuclear Motion. J. Phys. B: Mol. Opt. Phys.42, 134001. 10.1088/0953-4075/42/13/134001
2
BarthI.ManzJ. (2006). Periodic Electron Circulation Induced by Circularly Polarized Laser Pulses: Quantum Model Simulations for Mg Porphyrin. Angew. Chem. Int. Ed.45, 2962–2965. 10.1002/anie.200504147
3
BornM.HuangK. (1954). Dynamical Theory of crystal Lattices. Clarendon Press, Oxford.
4
BornM.OppenheimerR. (1927). Zur Quantentheorie der Molekeln. Ann. Phys.389, 457–484. 10.1002/andp.19273892002
5
CalegariF.AyusoD.TrabattoniA.BelshawL.De CamillisS.AnumulaS.et al (2014). Ultrafast Electron Dynamics in Phenylalanine Initiated by Attosecond Pulses. Science346, 336–339. 10.1126/science.1254061
6
CederbaumL. S.ZobeleyJ. (1999). Ultrafast Charge Migration by Electron Correlation. Chem. Phys. Lett.307, 205–210. 10.1016/S0009-2614(99)00508-4
7
DespréV.MarciniakA.LoriotV.GalbraithM. C. E.RouzéeA.VrakkingM. J. J.et al (2015). Attosecond Hole Migration in Benzene Molecules Surviving Nuclear Motion. J. Phys. Chem. Lett.6, 426–431. 10.1021/jz502493j
8
DingH.JiaD.ManzJ.YangY. (2017). Reconstruction of the Electronic Flux during Adiabatic Attosecond Charge Migration in HCCI+. Mol. Phys.115, 1813–1825. 10.1080/00268976.2017.1287967
9
EyringH.WalterJ.KimballG. E. (1944). Quantum Chemistry. New York: Wiley.
10
HeilbronnerE.MuszkatK. A.SchäublinJ. (1971). An Estimate of the Interatomic Distances in Monohaloacetylene Radical Cations from Photoelectron-Spectroscopic Data. Hca54, 58–76. 10.1002/hlca.19710540107
11
JiaD.ManzJ.PaulusB.PohlV.TremblayJ. C.YangY. (2017a). Quantum Control of Electronic Fluxes during Adiabatic Attosecond Charge Migration in Degenerate Superposition States of Benzene. Chem. Phys.482, 146–159. 10.1016/j.chemphys.2016.09.021
12
JiaD.ManzJ.YangY. (2019a). De- and Recoherence of Charge Migration in Ionized Iodoacetylene. J. Phys. Chem. Lett.10, 4273–4277. 10.1021/acs.jpclett.9b01687
13
JiaD.ManzJ.YangY. (2017b). Generation of Electronic Flux during the Femtosecond Laser Pulse Tailored to Induce Adiabatic Attosecond Charge Migration in. J. Mod. Opt.64, 960–970. 10.1080/09500340.2016.1269216
14
JiaD.ManzJ.YangY. (2019b). Timing the Recoherences of Attosecond Electronic Charge Migration by Quantum Control of Femtosecond Nuclear Dynamics: A Case Study for HCCI+. J. Chem. Phys.151, 244306. 10.1063/1.5134665
15
KannoM.KonoH.FujimuraY. (2006). Control of π-Electron Rotation in Chiral Aromatic Molecules by Nonhelical Laser Pulses. Angew. Chem. Int. Ed.45, 7995–7998. 10.1002/anie.200602479
16
KannoM.KonoH.FujimuraY.LinS. H. (2010). Nonadiabatic Response Model of Laser-Induced Ultrafastπ-Electron Rotations in Chiral Aromatic Molecules. Phys. Rev. Lett.104, 108302. 10.1103/physrevlett.104.108302
17
KrausP. M.MignoletB.BaykushevaD.RupenyanA.HornýL.PenkaE. F.et al (2015). Measurement and Laser Control of Attosecond Charge Migration in Ionized Iodoacetylene. Science350, 790–795. 10.1126/science.aab2160
18
LeforestierC.BisselingR. H.CerjanC.FeitM. D.FriesnerR.GuldbergA.et al (1991). A Comparison of Different Propagation Schemes for the Time Dependent Schrödinger Equation. J. Comput. Phys.94, 59–80. 10.1016/0021-9991(91)90137-a
19
LiH.MignoletB.WachterG.SkruszewiczS.ZherebtsovS.SüßmannF.et al (2015). Coherent Electronic Wave Packet Motion inC60Controlled by the Waveform and Polarization of Few-Cycle Laser Fields. Phys. Rev. Lett.114, 123004. 10.1103/PhysRevLett.114.123004
20
MarcusR. A. (1956). On the Theory of Oxidation‐Reduction Reactions Involving Electron Transfer. I. J. Chem. Phys.24, 966–978. 10.1063/1.1742723
21
MayV.KühnO. (2011). Charge and Energy Transfer Dynamics in Molecular Systems, Third, Revised and Enlarged Edition. Wiley-VCH, Weinheim.
22
Mendive-TapiaD.VacherM.BearparkM. J.RobbM. A. (2013). Coupled Electron-Nuclear Dynamics: Charge Migration and Charge Transfer Initiated Near a Conical Intersection. J. Chem. Phys.139, 044110. 10.1063/1.4815914
23
MineoH.LinS. H.FujimuraY. (2014). Vibrational Effects on UV/Vis Laser-Driven π-electron Ring Currents in Aromatic Ring Molecules. Chem. Phys.442, 103–110. 10.1016/j.chemphys.2014.02.011
24
MineoH.PhanN.-L.LaD.-K.FujimuraY. (2021). Theoretical Study of Dynamic Stark-Induced π-Electron Rotations in Low-Symmetry Aromatic Ring Molecules beyond the Frozen Nuclear Approximation. J. Phys. Chem. A.125, 1476–1489. 10.1021/acs.jpca.0c10216
25
MineoH.YamakiM.TeranishiY.HayashiM.LinS. H.FujimuraY. (2012). Quantum Switching of π-Electron Rotations in a Nonplanar Chiral Molecule by Using Linearly Polarized UV Laser Pulses. J. Am. Chem. Soc.134, 14279–14282. 10.1021/ja3047848
26
RemacleF.LevineR. D. (1999). Charge Migration and Control of Site Selective Reactivity: The Role of Covalent and Ionic States. J. Chem. Phys.110, 5089–5099. 10.1063/1.478406
27
RemacleF.LevineR. D.RatnerM. A. (1998). Charge Directed Reactivity:. Chem. Phys. Lett.285, 25–33. 10.1016/S0009-2614(97)01314-6
28
RemacleF.NestM.LevineR. D. (2007). Laser Steered Ultrafast Quantum Dynamics of Electrons in LiH. Phys. Rev. Lett.99, 183902. 10.1103/PhysRevLett.99.183902
29
UlusoyI. S.NestM. (2012). Remarks on the Validity of the Fixed Nuclei Approximation in Quantum Electron Dynamics. J. Phys. Chem. A.116, 11107–11110. 10.1021/jp304140r
30
WeinkaufR.SchanenP.MetsalaA.SchlagE. W.BürgleM.KesslerH. (1996). Highly Efficient Charge Transfer in Peptide Cations in the Gas Phase: Threshold Effects and Mechanism. J. Phys. Chem.100, 18567–18585. 10.1021/jp960926m
31
WeinkaufR.SchlagE. W.MartinezT. J.LevineR. D. (1997). Nonstationary Electronic States and Site-Selective Reactivity. J. Phys. Chem. A.101, 7702–7710. 10.1021/jp9715742
32
WernerH.-J.KnowlesP. J.KniziaG.ManbyF. R.SchützM.CelaniP.et al (2012). Molpro, Version 2012.1, a Package of Ab Initio Programs. Available at http://www.molpro.net (accessed Feb 18, 2016)
33
WörnerH. J.ArrellC. A.BanerjiN.CannizzoA.CherguiM.DasA. K.et al (2017). Charge Migration and Charge Transfer in Molecular Systems. Struct. Dyn.4, 061508. 10.1063/1.4996505
34
YamakiM.MineoH.TeranishiY.LinS. H.FujimuraY. (2016). Quantum Control of Coherent ρρ‐Electron Dynamics in Chiral Aromatic Molecules. Jnl Chin. Chem. Soc63, 87–92. 10.1002/jccs.201500043
35
YudinG. L.BandraukA. D.CorkumP. B. (2006). Chirped Attosecond Photoelectron Spectroscopy. Phys. Rev. Lett.96, 063002. 10.1103/PhysRevLett.96.063002
Summary
Keywords
frozen nuclei approximation, ultrashort laser pulses, nuclear quantum dynamics, electronic excitation, population transfer
Citation
Jia D and Yang Y (2022) Systematic Investigation of the Reliability of the Frozen Nuclei Approximation for Short-Pulse Excitation: The Example of HCCI+. Front. Chem. 10:857348. doi: 10.3389/fchem.2022.857348
Received
18 January 2022
Accepted
07 February 2022
Published
16 March 2022
Volume
10 - 2022
Edited by
Yuichi Fujimura, Tohoku University, Japan
Reviewed by
Hirobumi Mineo, Ton Duc Thang University, Vietnam
Manabu Kanno, Tohoku University, Japan
Updates
Copyright
© 2022 Jia and Yang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yonggang Yang, ygyang@sxu.edu.cn
This article was submitted to Physical Chemistry and Chemical Physics, a section of the journal Frontiers in Chemistry
Disclaimer
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.