Abstract
We theoretically study pulse phase and helicity effects on ultrafast magnetic field generation in intense bichromatic circularly polarized laser fields. Simulations are performed on the aligned molecular ion H2+ from numerical solutions of corresponding time-dependent Schrödinger equations. We demonstrate how electron coherent resonant excitation influences the phase and helicity of the optically induced magnetic field generation. The dependence of the generated magnetic field on the pulse phase arises from the interference effect between multiple excitation and ionization pathways, and is shown to be sensitive to molecular alignment and laser polarization. Molecular resonant excitation induces coherent ring electron currents, giving enhancement or suppression of the phase dependence. Pulse helicity effects control laser-induced electron dynamics in bichromatic circular polarization excitation. These phenomena are demonstrated by a molecular attosecond photoionization model and coherent electron current theory. The results offer a guiding principle for generating ultrafast magnetic fields and for studying coherent electron dynamics in complex molecular systems.
1 Introduction
Imaging and manipulating molecular electron dynamics is one of the main goals in photophysical processes and photochemical reactions. Advances in synthesizing ultrashort intense laser pulses [1, 2] allow one to visualize and control electrons on their natural attosecond (1 as = 10−18 s) timescale and sub-nanometer dimension [3–6]. One important application of ultrashort circularly polarized attosecond pulses is to produce strong magnetic field pulses from electronic ring currents in atomic and molecular systems [7–13]. By creating unidirectional constant valence-type electronic currents in molecules with circularly polarized UV laser pulses, static magnetic fields [7–9] can be efficiently generated by the excitation of resonant degenerate orbitals. These laser-induced magnetic fields are much larger than those obtained by traditional static field methods [14]. In [8], it has been found that for the hydrogen-like atom, the existence of ring currents is related to the presence of the states having nonzero magnetic orbital momentum magnetic quantum numbers. Surprisingly, the strongest magnetic field originates from the orbital in the hydrogen-like atom, which can be prepared via resonant transition. One also finds that, in general, ring electronic currents are dependent on the symmetry of the molecular orbitals. The helicity of driving circularly polarized pulses [15] can be used to reconstruct attosecond charge migration [16–19]. Linearly polarized laser pulses can also induce excited ring currents by controlling the rotation direction of π electrons in planar/nonplanar aromatic molecules [20–22]. We have proposed methods previously to create “spinning” continuum electrons which can be generated and remain localized on sub-nanometer molecular dimensional scales [23, 24], offering a way to produce high-order harmonic generation (HHG). Time-dependent circular coherent electron wave packets (CEWPs) and currents are created as superposition of bound-continuum states. They thus become the source of intense time-dependent internal magnetic fields generated on attosecond timescale. The induced attosecond magnetic fields have been shown to be a function of various laser pulse parameters, such as the pulse intensity, wavelength, and duration [25, 26], thus providing new tools for control of ultrafast optical magnetism generation [27–31].
Investigating ultrafast electron dynamics by bichromatic circularly polarized attosecond laser pulses with corotating or counter-rotating components has been attracting considerable attention in the field of light–matter interactions. It has already been shown that counter-rotating intense ultrafast circularly polarized pulses can induce re-collision, thus ensuring efficient HHG [32–35], the new source of circularly polarized X-ray attosecond pulses. These counter-rotating laser fields are now being adopted to produce circularly polarized HHG with nonzero initial angular momenta [36–38]. With counter-rotating circularly polarized laser pulses, the technique of double optical gating can be efficiently employed for producing isolated elliptically polarized attosecond pulses [39]. Bichromatic laser fields have also been adopted to probe atomic and molecular structure by photoelectron momentum distributions [40]. By combination of two circularly polarized attosecond ultraviolet (UV) pulses, spiral electron vortices in photoionization momentum distributions have been predicted theoretically in both atomic [41–43] and molecular systems [44–46], which are shown to be sensitive to the helicity of the bichromatic fields. Recent experiments have demonstrated this fact by focusing on multiphoton femtosecond ionization of potassium atoms [47, 48]. Most recently, above-threshold ionization obtained previously by a bicircular field has been reported [49–52].
In this work, we present attosecond magnetic field generation and electron currents under molecular resonant excitation in bichromatic attosecond circular polarization processes. Such ultrafast attosecond pulses have been generated by current laser techniques from circularly polarized HHG [53–55]. Numerical simulations are performed on the aligned molecular ion at equilibrium by numerically solving the corresponding three-dimensional (3D) time-dependent Schrödinger equation (TDSE). Ultrafast magnetic field generation has been studied previously by bichromatic circularly polarized laser pulses [56, 57]. It has shown that the interference effect between multiple ionization pathway influences the magnetic field generation. However, the effect of the coherent electron currents in bicircular magnetic field processes under molecular resonant excitation has not been presented. We focus here on the pulse phase and helicity-dependent magnetic field generation. We demonstrate molecular resonant excitation effects by comparing the dependence of the generated ultrafast magnetic field on the relative carrier-envelope phase (CEP) and the helicity of pulses at different molecular alignments. It is also found that attosecond charge migration arising from coherent resonant excitation induces coherent electron ring currents in molecules, leading to an absence of the CEP dependence. Induced electron ring currents resulting from molecular coherent resonant excitation are shown to be an important factor in bichromatic magnetic field generation. These results allow to control ultrafast magnetic fields, leading to molecular attosecond charge migration dynamics. Since molecular vibrational and rotational effects occur on the femtosecond (1 fs = 10−15 s) and picosecond (1 ps = 10−12 s) timescales, fixed nuclei simulations are valid and used to describe ultrafast magnetic field generation processes on the attosecond timescale.
The article is organized as follow: We briefly describe the computational method for solving TDSEs of the aligned molecular ion to simulate electron currents and magnetic field generation in Section 2. The results of ultrafast magnetic fields by intense bichromatic circularly polarized attosecond XUV laser pulses are presented and discussed in Section 3. In Section 4, we finally summarize our findings. Throughout this article, atomic units (au) which are defined by setting are used unless otherwise stated.
2 Numerical Methods
For the aligned molecule ion within Born–Oppenheimer approximation static nuclear frames, the corresponding 3D TDSE in the cylindrical coordinate reads aswith the Laplacianwhere is the electron-nuclear potential in which we fixed the internuclear separation to the equilibrium separation and |R| = 2 au. The circularly polarized laser pulse propagates along the z axis, perpendicular to the (x,y) plane, with and , as shown in Figure 1. The radiative interaction between the laser field and the electron is described in the length gauge for circularly polarized pulses of frequencies and .where is the laser polarization direction, and the symbol denotes the helicity of combined fields, that is, corotating () or counter-rotating () components. and are CEPs of the pulses and , respectively. A smooth pulse envelope for maximum amplitude E, intensity , and duration is adopted, where one optical cycle period .
FIGURE 1
The 3D TDSE in Eq. 1 is propagated by a second-order split operator method which conserves unitarity in each time step combined with a fifth-order finite difference method and Fourier transform technique in the spatial steps , , and [59, 60]. The initial electron wave function is prepared in the ground state calculated by propagating an initial appropriate wave function in imaginary time using the zero-field TDSE in Eq. 1. The time step is taken to be au = 0.24 as. The spatial discretization is au for a radial grid range au (6.77 nm) and |z| ≤ 32 au (1.69 nm), and the angle grid size radian. To prevent unphysical effects due to the reflection of the wave packet from the boundary, we multiply by a “mask function” or absorber in the radial coordinate ρ with the form . For all results reported here, we set the absorber domain at ρa = ρmax−ρabs = 104 au with ρabs = 24 au, exceeding well the field-induced electron oscillation of the electron [25]. The time-dependent probability current density is defined by the quantum expression in the length gauge,where is the exact Born–Oppenheimer (static nuclei) electron wave function obtained from the TDSE in Eq. 1, and is the momentum operator in cylindrical coordinates. Then the corresponding time-dependent magnetic field is calculated using the following classical Jefimenko equation [61]:where is the retarded time and NA−2 (6.692 × 10−4 au) is the permeability of free space. Units of are teslas (1 T = 104 Gauss) if the elementary charge e is in Coulombs. For the static zero-field time-independent conditions occurring after the pulse duration, Eq. 5 reduces to the classical Biot–Savart law, that is, [61]. Note that the retardation effects due to r/c = 0.35 attoseconds (where for an estimate r = R = 2 au. is used) are negligible.
Of note is that equation Eq. 4 defines the quantum probability current (not the electric current) as defined in any quantum mechanics textbook. The electron electric current used in the Biot–Savart law in Eq. 5 is therefore . This explains the sign (−) in Eq. 5, e = 1 in atomic units.
3 Results and Discussions
The ground and excited states of ion and its electron potentials, , , , , , etc., are well documented in [65]. Numerical solutions of the TDSE for aligned with the laser polarization (x–y) plane, that is, to the x-axis with the electric field vector rotating in the (x,y) plane (Figure 1A), are used for obtaining time-dependent probabilities P(t) of the and states at equilibrium R = 2 au. excited by a 70-nm pulse at two different intensities, I = 2× 1014 W/cm2 and I = 1×1015 W/cm2, by a five-cycle pulse (one cycle = 0.234 fs = 234 as). Thus, at , Figure 2A, one sees no significant ionization, whereas at I = 1 × 1015 W/cm2, Figure 2B, the state is 85% occupied, and still one sees little ionization. In conclusion, short few-cycle intense and resonant pulses contribute little ionization with major excitation of the resonant state, such as the state at 70 nm and some Rydberg states above the state shown in 1.
FIGURE 2
We investigate laser-induced highly nonlinear optical effects using pairs of bichromatic circularly polarized laser pulses. We use λ1 = 70 nm (ω1 = 0.65 au) circularly polarized pulse in combination with λ2 = 35 nm (ω2 = 2ω1 = 1.3 au) circularly polarized pulse. Pairs of circularly polarized harmonics of different frequency and helicity can easily be prepared by a combination of pairs of counter-rotating circularly polarized laser pulses at different frequencies [62]. The molecular ion is aligned along the x/z axis, the two X-ray ultra violet (XUV) pulses with their field polarization vectors in the (x,y) plane propagate along the z axis, as illustrated in Figure 1A, for in-plane excitation, and Figure 1B for around axis excitation. With the pulse frequency and the energy difference au, a resonant excitation with the ground state and the excited electronic states occurs. Moreover, with , such bichromatic laser pulses can also produce CEWPs with the same kinetic energies, , by combination of multiple multiphoton transitions to Rydberg and the continuum state, thus leading to electron interference between the two ionization pathways. Since the induced electron currents are localized in the laser () polarization plane, the generated magnetic field is concentrated along the z axis. We therefore only present the results of the magnetic field B along the z axis.
As illustrated in Figure 1C, the coherent resonant excitation between the ground state and the excited state by the λ1 = 70 nm circularly polarized laser pulse is dependent on the molecular alignment. In the case of , the in-plane degenerates perpendicular excitation, for different angular momenta dominates, and the electron density distribution evolves nearly perpendicular to the molecular axis [63], whereas in the case of , the produced electron wave packets move around the molecular axis due to individual transitions or [64]. Meanwhile, two and single transitions to the continuum () occur, leading to interference effects between these two excitation ionization channels. We show how coherent resonant excitation influences the interference effect on the magnetic field generation. We present the effects of the pulse phase on the attosecond magnetic field generation under the resonant excitation with various molecular alignments. Processes with different helicities ϵ, that is, corotating () and counter-rotating () pulse combinations, are compared as well.
3.1 Dependence of Generated Magnetic Fields on Circular Error Probability and Molecular Alignments
We first present the counter-rotating () dynamics. Figure 3 shows results of the generated maximum magnetic field B at the molecular center by intense single and bichromatic circularly polarized attosecond pulses described in Eq. 3 [57]. Two cases of molecular excitation–ionization processes in the molecular ion aligned along a) the x axis, , and b) the z axis, , are compared for the bichromatic laser pulses with their field polarization vectors in the (x,y) plane, as illustrated in Figure 1. The corresponding generated magnetic fields mainly lie in both cases along the z axis. The pulse wavelengths are, respectively, λ1 = 70 nm (ω1 = 0.65 au) and λ2 = 35 nm (ω2 = 1.3 au). We always fix the pulse intensities I = 1 × 1014 W/cm2 (E = 0.0534 au) and durations Tlp = 48.3 au = 1.16 fs, corresponding to 580 as FWHM (full width at half maximum), that is, for pulse and for pulse. We also set the pulse CEP and vary the CEP from 0 to , that is, the relative CEP is . Figure 3 shows the CEP ϕ dependence of the generated magnetic fields with different molecular alignments.
FIGURE 3
3.1.1 Resonant Excitation in the Case of In-Plane
For the case of in-plane excitation with the molecule axis parallel to the laser () polarization plane as shown in Figure 1A, the generated magnetic field B at the molecular center is critically sensitive to the relative pulse CEP ϕ, as shown in Figure 3A. It is found that as the CEP ϕ varies from 0 to , B oscillates periodically as . The maximum and minimum values of the magnetic field are T and 0.165 T at 0 (or and π, respectively. For comparison, we also simulate the results of the induced magnetic field B with single-color circularly polarized attosecond pulses at separated wavelengths λ = 70 and 35 nm. The other parameters are the same as those used in Figure 3A. For the two single-pulse processes, the corresponding maximum strengths of induced magnetic fields at the molecular center are, respectively, T and T. Due to the resonant excitation, the magnetic field induced by the resonant 70 pulse is much stronger than that induced by the nonresonant 35-nm pulse, . The sum value of the generated magnetic fields is slightly smaller than the maximum field T of the bicircular process. Comparing to the single-color circular processes, one notes that the sensitivity of the induced magnetic field B to the CEP ϕ in Figure 3A results from the interference effects between λ1 = 70 nm and λ2 = 35 nm optical processes.
From Eq. 4, 5, for a simple ring current flowing in the ring having the radius r, one can derive a simple relation:where is the probability density and v is the electron speed. This simple relation can be obtained in the following way: from the quantum definition of the current given in Eq. 4 in which we assume that ψ is the plane wave, that is, , where is the electron momentum tangent to the ring, one gets after performing the derivation a simple relation as shown in [58]. Next, for the case of the electron current mainly localized in a plane one gets from Eq. 5, that is, , we thus derived Eq. 6.
Thus, in Eq. 6, the magnetic field is proportional to the electron probability density and the electron velocity which originates from the electron current density [58]. Density for the case of resonant transitions corresponds to the excited state transition probability, which is determined by the intrinsic transition dipole of molecules and the electric field strength as discussed in detail in Supplementary Appendix A2. Therefore, the dependence of the generated magnetic field in Figure 3 on the pulse CEP ϕ mainly comes from the transition probability which is influenced by the interference effect between color λ1 = 70 nm and λ2 = 35 nm pulse nonlinear optical responses, as presented in Supplementary Appendix A1.
With the bichromatic counter-rotating circularly polarized pulse at λ1 = 70 nm and λ2 = 35 nm, two nonlinear responses can be triggered. By the λ1 = 70 nm pulse, resonance-enhanced excitation ionization occurs where the molecule is resonantly excited from the state to the degenerate state with via a perpendicular transition. Meanwhile, the absorptions of two photons and one photon give rise to photoelectron wave packets with the same kinetic energies in the continuum. The total photoionization probability is the sum of the two ionization excitation processes and their interference. As shown in Supplementary Appendix A1, the two ionization probabilities and are insensitive to the pulse phases, whereas the interference term depends on the relative phase of the laser pulses and on the relative phase of the electron wave packets in the continuum. As a result, the total excitation probability density is also a cosine function of the relative pulse phase ϕ with the form . The interference effect between the two processes modulates the total transition and current probabilities. Combining Eq. 6 with Eqs 19 and 20 in Supplementary Appendix A1, one obtains , giving rise to a CEP dependence of the generated magnetic field, as illustrated in Figure 3A (blue diamond), B (r = 0).
3.1.2 Resonant Excitation in the Case of Around-Axis
For around R axis excitation, Figure 3B shows phase-dependent magnetic field B generation in the process of molecule axis perpendicular to the laser () polarization plane, leading to around-axis currents. Comparing to the in-plane case in Figure 3A, one sees that the oscillation of the magnetic field B with the relative phase ϕ is, however, strongly suppressed. As shown in Figure 3B, the magnetic field is almost insensitive to the phase ϕ, with a constant value of T. The magnetic fields generated by the individual bichromatic pulses are less than the sum of the two single 70 and 35 nm processes, that is, T, where T and T. The bichromatic magnetic fields mainly arise from the single-photon 70 nm process, , whereas the contribution from the absorption of single 35 nm () or two 70 nm () photons is weak and negligible, as shown in Figure 3B. This implies that the interference effect between the two ionization processes does not influence the bichromatic magnetic field generation. Of note is that for the single / pulse case, the generated magnetic fields are also dependent on the molecular alignments and the pulse wavelength. At λ = 70 nm, , the perpendicular case is stronger, whereas at 35 nm, the generated magnetic fields are nearly equivalent, . The difference indicates essentially the importance of the around R axis ring electron currents in the resonant intermediate electronic state.
For the optical responses in Figure 3B of the molecule perpendicular to the laser () polarization plane, , Figure 1B, by a single 70-nm circularly polarized pulse, resonant excitation between the ground state and the excited state with magnetic quantum number leads to ring electron currents in molecules [64, 67]. As shown in Supplementary Appendix A2, using the electronic angular continuity equation [16, 17],One obtains the laser-induced electron current in the following form
Eq. 8 combined with Eqs. 21 and 26 shows that for the superposition electron state , the time-dependent electronic density and the induced angular current evolve in time with the electron coherence period of . Therefore, the total magnetic field B in the bichromatic circularly laser field is mainly generated from the electron currents in the coherent electron wave packet in Eq. 21 by the λ =70 nm pulse. The contributions from the continuum electron wave packets with energy by the 35-nm pulse are clearly very weak, that is, , and therefore negligible. As a result, altering the relative CEP ϕ does not lead to a variance of the generated magnetic field B.
3.1.3 Laser-Induced Electron Currents in Molecules
For the magnetic field B generated by coherent electron currents in molecules, the evolution of the induced electron currents with time t is determined by the coherent resonant excitation and the molecular alignments. Figure 4 illustrates snapshots of angular electron probability currents at the center of obtained from Eq. 4 at different times t (unit of ) for the molecule perpendicular to the laser polarization plane, , Figure 1, leading to around R axis currents, . The bichromatic counter-rotating () circularly polarized XUV pulse at CEPs is used to excite the molecule. It is found that the induced electron currents are asymmetric with respect to the molecular center and rotate with a period of τ around the z or molecular R axis in the polarization plane with an anticlockwise direction, as predicted in Eq. 8. This confirms that the ring electron currents mainly arise from the coherent resonant excitation between the ground state and the excited state by the 70-nm circularly polarized pulse [64]. From such coherent electron currents, the generated magnetic fields are unidirectional, along the axis. At the molecular center , the magnetic field is the sum of those at the two nuclei, . In the bichromatic circular polarized processes, the contributions of the magnetic field generation from the 70 nm pulse are dominant, which do not depend on the pulse phase. Therefore, varying the relative pulse CEP ϕ does not influence the generated magnetic field B, as illustrated in Figure 3B.
FIGURE 4
For comparison, in Figure 5, we also plot the in-plane electron probability current at the center of the molecule aligned along the x axis, parallel to the laser polarization () plane, in Figure 1A, of bichromatic counter-rotating circularly polarized pulses. Combining Figure 4 with Figure 5, one sees that these time-dependent electron currents are sensitive to the molecular alignments. The joint resonant excitations induce coherent electron currents. The induced electron probability currents are mainly localized along the molecular internuclear axis, that is, the resonant perpendicular atomic () transitions dominate during the excitation processes. In this case, the coherent exited electronic state is given bywhere and . Then the corresponding interference term in the time-dependent electron density becomes (in cylindrical coordinates () as shown in [9, 64] and in Supplementary Appendix A2) and the time-dependent currentwhere is the absolute value of (see Eq. 23). It is found that the coherent density and current superposition terms in Eqs. 10 and 11 follow the forms and . As a result, the coherent wave packets due to the superposition oscillate along the y direction, perpendicular to the molecular axis, Figure 1C, whereas the corresponding currents oscillate mainly along the molecular axis, Figure 4. The generated magnetic field at the two molecular nuclear centers has opposite phases [70]. Their overlap leads to a weak magnetic field at the molecule center . The generated magnetic field in Figure 3A therefore mainly results from the coherent excitation by bichromatic circularly polarized laser pulses.
FIGURE 5
3.2 Influence of the Pulse Helicity
We next study the process with a bichromatic corotating () circularly polarized laser pulse. Figure 6 shows the maximum generated magnetic field B at various relative pulse phases ϕ. The other laser parameters are the same as those used in Figure 3. It was found that similar results are obtained in Figure 6, as the counter-rotating () case in Figure 3. For the case of , the magnetic field B varies as a sine function of the phase ϕ, as predicted in Supplementary Appendix A1. A phase shift occurs in the CEP ϕ-dependent magnetic field generation by the corotating bichromatic pulses. This mainly results from the different electron dynamics induced by the two pulses with opposite helicity. It was also found for the corotating () case the generated magnetic field depends on the molecular alignment. In Figure 6B, we see that for the case of , the total generated magnetic field (blue diamond) is shown to be almost insensitive to the relative CEP ϕ, similar as in Figure 3B. At various pulse phases ϕ, the magnetic field T is obtained. The resonant excitation with gives rise to unidirection ring electron currents, as illustrated in Figure 3. The contribution from the coherent electron wave packets in Eq. 21 dominates. The effects of the multiple pathway ( and ) excitation interference can be neglected. Consequently, the generated magnetic field does not depend on the pulse CEP ϕ. The independence of the generated magnetic field on the pulse helicity ϵ also confirms the importance of the charge migration. This offers an approach to explore molecular structure and orbitals.
FIGURE 6
It should be noted that in the case of in the general pulse with its field polarization vectors in the plane, the strength of the generated magnetic field is slightly sensitive to the pulse helicity ϵ. As shown in Figure 6A, at for corotating cases, the maximum magnetic field is T at , which is slightly weaker than that at for counter-rotating cases in Figure 3A, where T at . Eq. 6 predicts that the induced magnetic field, in general, is proportional to the ratio of the electron velocity v and inversely proportional to , where r is the radius of the excitation state of an electron under the influence of a strong laser field. The difference between the corotating and counter-rotating generated magnetic fields in Figures 3A,6A results from the laser-induced dynamics, which depends on the helicity of driving pulses.
According to the classical laser-induced electron motion models [68, 69], the electron velocity and radius are determined by the pulse amplitude E, frequency , and helicity ϵ, which are given byand the corresponding electron displacements arewhere is the ionization time. Figure 7 shows bichromatic circularly polarized laser field vectors with corotating and counter-rotating components and corresponding induced electron displacements, where zero initial position of one ionized electron is assumed, that is, . Such electron displacements can be measured by interfering spirals in photoelectron momentum distributions [71]. With the XUV pulses, a multiphoton ionization process occurs since the Keldysh parameter [72]. The initial electron velocities are nearly equivalent with . As shown in Figure 7, for the corotating case, the radius of the continuum electron is much larger than that of the counter-rotating case. The electron moves way quickly in the corotating field, whereas the counter-rotating field restricts the ionized electron around the molecular center. Because in Eq. 5, 6, weaker maximum magnetic field is generated in the corotating case in Figure 5A. The dependence of generated magnetic fields on the pulse helicity ϵ reflects the laser-induced electron dynamics in bichromatic fields. The present simulation confirms that steering the radius of the induced electron currents allows to control generated magnetic fields with bichromatic circularly polarized pulses.
FIGURE 7
4 Conclusion
We present the ultrafast magnetic field generation in molecules from one electron molecular TDSE simulation under effects of coherent resonant excitation in bichromatic co- and counter-rotating circularly polarized laser fields. Numerical results are obtained for the aligned molecular ion , which can be fully and exactly studied. We evaluate the generated magnetic field at the molecular center at different relative CEP ϕ of the two circularly polarized pulses with both counter-rotating and corotating combinations. It is found that altering the CEP ϕ varies the maximum values of the generated magnetic field, , which is shown to be dependent on the molecule alignment, and is maximum for in-plane bichromatic excitation.
In a bichromatic (frequency
) circularly polarized field, a
resonant excitation of individual magnetic components
is triggered by the resonant
(70 nm) pulse. As a result, coherent electron currents between the ground
state and the degenerate excited
electronic state are induced in molecules. We compare two nonlinear responses for different molecular alignments:
• For the around-axis case, , the molecule aligned along the z axis, perpendicular to the laser () polarization plane, Figure 1B, and the magnetic field mainly results from the coherent electron currents which are induced by one resonant photon around R. The contribution from two photons or one photon is negligible. Varying the relative CEP ϕ does not influence the generated magnetic field. Similar results are obtained for both counter-rotating , Figure 3B, and corotating , Figure 6B, excitations, confirming the main role of the around-axis circular coherent electron current.
• For the in-plane case, , the molecule aligned parallel to the laser () polarization plane, Figure 1A, and the magnetic field generation is shown to be strongly dependent on the pulse phase ϕ, confirming the effects of multiple pathway ionization. The in-plane coherent electron currents generate the magnetic fields in the two nuclei with opposite phases and evolve perpendicular to the molecular R axis. Their superposition suppresses the magnetic field generation at the molecular center, Figure 3. The results we present arise mainly from the total electron currents by one photon absorption, that is, the resonant excitation, and two or one transition to the continuum. The generated magnetic field depends on the photoionization probability which is a function of the relative CEP ϕ Eq. 20. As a consequence of interference effects between the and photoionization pathways, altering their relative phase ϕ gives rise to a modulation of the generated magnetic field at the molecular center with forms for counter-rotating and for corotating combinations. It is also shown that the maximum of the magnetic field for the counter-rotating case is larger than that for the corotating case because of the difference in laser-induced electron displacements, which depends strongly on the helicity of the driving pulse, Figure 6.
The present results in principle provide the importance of coherent electron dynamics and of control magnetic fields by bichromatic circularly polarized laser pulses. The dependence of the generated magnetic field on the relative phase and helicity of driving laser pulses also allows to characterize the property of laser pulses and probe coherent electron currents and to charge migration in molecules. Although a simple single electron molecular ion is used, similar electron dynamics phenomena should be predicted in more complex molecular systems [7, 27, 30, 73], thus offering an approach for controlling ultrafast magnetic field generation.
The above laser-induced molecular magnetic field generation on the electron’s quantum timescale, the asec, was studied in the Born–Oppenheimer Approximation, that is, with static nuclei. Nuclear motion effects, that is, non–Born–Oppenheimer, are now being pursued on the near femtosecond timescale in order to include nuclear motion effects with bound and dissociation molecular states [83], de-and re-coherence in charge migration [84] and isotope effects in HD+ ultrafast ionization [85]. In the case of laser pulses propagating perpendicular to the molecular R-axis with the pulse electric fields in the molecular plane, Figure 1A, re-collision of electron currents with nuclei is an important nonlinear optical effect shown in Figure 6 to be examined in detail for moving nuclei. Finally, the strong magnetic fields generated by intense ultrafast laser pulses are expected to interact with the electron currents themselves. Proton beams have been shown recently to be useful tools to measure intense magnetic field directions generated by current solenoids [86], thus confirming that laser-generated magnetic fields can interact also with nuclei in matter.
Statements
Author contributions
AB as a research leader directed this research formulated research goals and contributed in the scientific interpretation of results. KY wrote codes, executed codes, and prepared graphics but died last March before completing this research. SC participated in initial algorithm preparation for computer codes and did final preparation.
Funding
This work is supported in part by the National Natural Science Foundation of China (Grant Nos. 11974007 and 11574117) and the NSERC-RGPIN2019-05291.
Acknowledgments
The authors also thank Compute Canada for access to massively parallel computer clusters, and the Natural Sciences and Engineering Research Council of Canada and the Fonds de Recherche du Québec-Nature et Technologies for supporting their research work.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2021.675375/full#supplementary-material
References
1.
KrauszF.IvanovM.Attosecond Physics. Rev Mod Phys (2009) 81:163, 234. 10.1103/revmodphys.81.163
2.
ChangZ.CorkumP. B.LeoneS. R.Attosecond optics and technology: progress to date and future prospects [Invited]. J Opt Soc Am B (2016) 33:1081. 10.1364/josab.33.001081
3.
ChelkowskiS.YudinG. L.BandraukA. D.Observing Electron Motion in Molecules. J Phys B At Mol Opt Phys (2006) 39:S409, S417. 10.1088/0953-4075/39/13/s14
4.
NiikuraHVilleneuveDMCorkumPB. Mapping Attosecond Electron Wave Packet Motion. Phys Rev Lett (2005) 94:083003. 10.1103/physrevlett.94.083003
5.
ShaoHCStaraceAF. Detecting Electron Motion in Atoms and Molecules. Phys Rev Lett (2010) 105:263201. 10.1103/physrevlett.105.263201
6.
VrakkingM. J. J.ElsaesserT.X-rays inspire electron movies. Nature Photon (2012) 6:645, 647. 10.1038/nphoton.2012.247
7.
BarthI.ManzJ.ShigetaY.YagiK.Unidirectional Electronic Ring Current Driven by a Few Cycle Circularly Polarized Laser Pulse: Quantum Model Simulations for Mg−Porphyrin. J Am Chem Soc (2006) 128:7043, 7049. 10.1021/ja057197l
8.
BarthIManzJ. Electric Ring Currents in Atomic Orbitals and Magnetic fields Induced by Short Intense Circularly Polarized π Laser Pulses. Phys Rev A (2007) 75:012510. 10.1103/physreva.75.012510
9.
BarthI.Serrano-AndrésL.SeidemanT.Nonadiabatic Orientation, Toroidal Current, and Induced Magnetic Field in BeO Molecules. The Journal of Chemical Physics (2008) 129:164303. 10.1063/1.2994737
10.
BarthIManzJ. Quantum Switching of Magnetic Fields by Circularly Polarized Optimized Laser Pulses. In: YamanouchiKGerberGBandraukAD, editors. Progress in Ultrafast Intense Laser Science VI, 2021. Springer Verlag Berlin (2010) p. 2.
11.
BarthI.ManzJ.Serrano-AndrésL.Quantum Simulations of Toroidal Electric Ring Currents and Magnetic fields in Linear Molecules Induced by Circularly Polarized Laser Pulses. Chemical Physics (2008) 347:263–271. 10.1016/j.chemphys.2007.09.037
12.
KöksalK.KoçF.The Effect of Twisted Light on the Ring-Shaped Molecules: The Manipulation of the Photoinduced Current and the Magnetic Moment by Transferring Spin and Orbital Angular Momentum of High Frequency Light. Computational and Theoretical Chemistry (2017) 1099:203, 208. 10.1016/j.comptc.2016.11.031
13.
ZhangGMurakamiMLefkidisGHübnerWGeorgeTF. Introduction to Ultrafast , F͡rom Femtosecond Magnetism to High-Harmonic Generation. CRC Press, Taylor & Francis Group (2020) 10.1201/9780429194832
14.
NobusadaKYabanaK. Photoinduced Electric Currents in Ring-Shaped Molecules by Circularly Polarized Laser Pulses. Phys Rev A (2007) 75:032518. 10.1103/physreva.75.032518
15.
JiaD.ManzJ.PaulusB.PohlV.TremblayJ. C.YangY.Quantum Control of Electronic Fluxes during Adiabatic Attosecond Charge Migration in Degenerate Superposition States of Benzene. Chemical Physics (2017) 482:146, 159. 10.1016/j.chemphys.2016.09.021
16.
HermannG.LiuC.ManzJ.PaulusB.Pérez-TorresJ. F.PohlV.et alMultidirectional Angular Electronic Flux during Adiabatic Attosecond Charge Migration in Excited Benzene. J Phys Chem A (2016) 120:5360, 5369. 10.1021/acs.jpca.6b01948
17.
HermannG.LiuC.ManzJ.PaulusB.PohlV.TremblayJ. C.Attosecond Angular Flux of Partial Charges on the Carbon Atoms of Benzene in Non-aromatic Excited State. Chemical Physics Letters (2017) 683:553, 558. 10.1016/j.cplett.2017.01.030
18.
DingH.JiaD.ManzJ.YangY.Reconstruction of the electronic flux during adiabatic attosecond charge migration in HCCI+. Molecular Physics (2017) 115:1813, 1825. 10.1080/00268976.2017.1287967
19.
JiaD.ManzJ.YangY.Generation of electronic flux during the femtosecond laser pulse tailored to induce adiabatic attosecond charge migration in. Journal of Modern Optics (2017) 64:960, 970. 10.1080/09500340.2016.1269216
20.
KannoMKonoHFujimuraYLinSH. Nonadiabatic Response Model of Laser-Induced Ultrafast π-electron Rotations in Chiral Aromatic Molecules. Phys Rev Lett (2010) 104:108302. 10.1103/physrevlett.104.108302
21.
MineoH.LinS. H.FujimuraY.Coherent π-electron dynamics of (P)-2,2′-biphenol induced by ultrashort linearly polarized UV pulses: Angular momentum and ring current. The Journal of Chemical Physics (2013) 138:074304. 10.1063/1.4790595
22.
KannoM.KonoH.FujimuraY.Laser-Control of Ultrafast π-Electron Ring Currents in Aromatic Molecules: Roles of Molecular Symmetry and Light Polarization. Applied Sciences (2018) 8:2347. 10.3390/app8122347
23.
YuanK.-J.BandraukA. D.Circularly polarized attosecond pulses from molecular high-order harmonic generation by ultrashort intense bichromatic circularly and linearly polarized laser pulses. J Phys B At Mol Opt Phys (2012) 45:074001. 10.1088/0953-4075/45/7/074001
24.
YuanKJLuHBandraukAD. High-order-harmonic Generation in Molecular Sequential Double Ionization by Intense Circularly Polarized Laser Pulses. Phys Rev A (2015) 92:023415. 10.1103/physreva.92.023415
25.
YuanKJBandraukAD. Attosecond-magnetic-field-pulse Generation by Intense Few-Cycle Circularly Polarized UV Laser Pulses. Phys Rev A (2013) 88:013417. 10.1103/physreva.88.013417
26.
YuanKJBandraukAD. Attosecond-magnetic-field-pulse Generation by Coherent Circular Molecular Electron Wave Packets. Phys Rev A (2015) 91:042509. 10.1103/physreva.91.042509
27.
ZhangG. P.HübnerW.LefkidisG.BaiY.GeorgeT. F.Paradigm of the Time-Resolved Magneto-Optical Kerr Effect for Femtosecond Magnetism. Nature Phys (2009) 5:499, 502. 10.1038/nphys1315
28.
ZhangGPGeorgeTF. Total Angular Momentum Conservation in Laser-Induced Femtosecond Magnetism. Phys Rev B (2008) 78:052407. 10.1103/physrevb.78.052407
29.
KimJ-WVomirMBigotJ-Y. Ultrafast Magnetoacoustics in Nickel Films. Phys Rev Lett (2012) 109:166601. 10.1103/physrevlett.109.166601
30.
BigotJ.-Y.Down to the nanometre scale. Nature Mater (2013) 12:283, 284. 10.1038/nmat3605
31.
ZhangXZhuXWangDLiLLiuXLiaoQet alUltrafast Oscillating-Magnetic-Field Generation Based on Electronic-Current Dynamics. Phys Rev A (2019) 99:013414. 10.1103/physreva.99.013414
32.
ZuoT.BandraukA. D.High-order Harmonic Generation in Intense Laser and Magnetic fields. J Nonlinear Optic Phys Mat (1995) 04:533, 546. 10.1142/s0218863595000227
33.
BandraukADLuHZ. Controlling Harmonic Generation in Molecules with Intense Laser and Static Magnetic fields: Orientation Effects. Phys Rev A (2003) 68:043408. 10.1103/physreva.68.043408
34.
LongS.BeckerW.McIverJ. K.Model Calculations of Polarization-dependent Two-Color High-Harmonic Generation. Phys Rev A (1995) 52:2262, 2278. 10.1103/physreva.52.2262
35.
MiloševićDBBeckerWKopoldR. Generation of Circularly Polarized High-Order Harmonics by Two-Color Coplanar Field Mixing. Phys Rev A (2000) 61:063403.
36.
FleischerA.KfirO.DiskinT.SidorenkoP.CohenO.Spin Angular Momentum and Tunable Polarization in High-Harmonic Generation. Nature Photon (2014) 8:543, 549. 10.1038/nphoton.2014.108
37.
MedišauskasLWraggJvan der HartHIvanovMY. Generating Isolated Elliptically Polarized Attosecond Pulses Using Bichromatic Counterrotating Circularly Polarized Laser fields. Phys Rev Lett (2015) 115:153001. 10.1103/PhysRevLett.115.153001
38.
MiloševićDB. Circularly Polarized High Harmonics Generated by a Bicircular Field from Inert Atomic Gases in the P State: A Tool for Exploring Chirality-Sensitive Processes. Phys Rev A (2015) 92:043827.
39.
ChangZ. Single Attosecond Pulse and Xuv Supercontinuum in the High-Order Harmonic Plateau. Phys Rev A (2004) 70:043802. 10.1103/physreva.70.043802
40.
MancusoCAHicksteinDDGrychtolPKnutRKfirOTongX-Met alStrong-field Ionization with Two-Color Circularly Polarized Laser fields. Phys Rev A (2015) 91:031402(R). 10.1103/physreva.91.031402
41.
Ngoko DjiokapJMHuSXMadsenLBManakovNLMeremianinAVStaraceAF. Electron Vortices in Photoionization by Circularly Polarized Attosecond Pulses. Phys Rev Lett (2015) 115:113004. 10.1103/physrevlett.115.113004
42.
Ngoko DjiokapJMMeremianinAVManakovNLHuSXMadsenLBStaraceAF. Multistart Spiral Electron Vortices in Ionization by Circularly Polarized UV Pulses. Phys Rev A (2016) 91:013408.
43.
LiM.ZhangG.KongX.WangT.DingX.YaoJ.Dynamic Stark Induced Vortex Momentum of Hydrogen in Circular fields. Opt Express (2018) 26:878. 10.1364/oe.26.000878
44.
YuanKJChelkowskiSBandraukAD. Photoelectron Momentum Distributions of Molecules in Bichromatic Circularly Polarized Attosecond UV Laser fields. Phys Rev A (2016) 93:053425. 10.1103/physreva.93.053425
45.
YuanK.-J.LuH.BandraukA. D.Photoionization of triatomic molecular ions H2+3 by intense bichromatic circularly polarized attosecond UV laser pulses. J Phys B At Mol Opt Phys (2017) 50:124004. 10.1088/1361-6455/aa72fa
46.
Ngoko DjiokapJMMeremianinAVManakovNLMadsenLBHuSXStaraceAF. Dynamical Electron Vortices in Attosecond Double Photoionization of H. Phys Rev A (2018) 98:063407. 10.1103/physreva.98.063407
47.
PengelDKerbstadtSJohannmeyerDEnglertLBayerTWollenhauptM. Electron Vortices in Femtosecond Multiphoton Ionization. Phys Rev Lett (2017) 118:053003. 10.1103/physrevlett.118.053003
48.
PengelDKerbstadtSEnglertLBayerTWollenhauptM. Control of Three-Dimensional Electron Vortices from Femtosecond Multiphoton Ionization. Phys Rev A (2017) 96:043426. 10.1103/physreva.96.043426
49.
HasovicEBeckerWMilosevicDB. Electron Rescattering in a Bicircular Laser Field. Opt Express (2016) 24:6413.
50.
MancusoCAHicksteinDDDorneyKMEllisJLHasovicEKnutRet alControlling Electron-Ion Rescattering in Two-Color Circularly Polarized Femtosecond Laser fields. Phys Rev A (2016) 93:053406. 10.1103/physreva.93.053406
51.
LiMJiangW-CXieHLuoSZhouYLuP. Strong-field Photoelectron Holography of Atoms by Bicircular Two-Color Laser Pulses. Phys Rev A (2018) 97:023415. 10.1103/physreva.97.023415
52.
YuanKJBandraukAD. Above-threshold Ionization in Molecules by Intense Multiple-Frequency Circularly Polarized Laser Pulses. Phys Rev A (2018) 98:023413. 10.1103/physreva.98.023413
53.
XieXScrinziAWickenhauserMBaltuškaABarthIKitzlerM. Internal Momentum State Mapping Using High Harmonic Radiation. Phys Rev Lett (2008) 101:033901. 10.1103/physrevlett.101.033901
54.
YuanKJBandraukAD. Single Circularly Polarized Attosecond Pulse Generation by Intense Few Cycle Elliptically Polarized Laser Pulses and Terahertz fields from Molecular media. Phys Rev Lett (2013) 110:023003. 10.1103/physrevlett.110.023003
55.
Hernández-GarcíaCDurfeeCGHicksteinDDPopmintchevTMeierAMurnaneMMet alSchemes for Generation of Isolated Attosecond Pulses of Pure Circular Polarization. Phys Rev A (2016) 93:043855. 10.1103/physreva.93.043855
56.
YuanKJBandraukAD. Attosecond-magnetic-field-pulse Generation by Electronic Currents in Bichromatic Circularly Polarized UV Laser fields. Phys Rev A (2015) 92:063401. 10.1103/physreva.92.063401
57.
GuoJYuanKJLuHBandraukAD. Spatiotemporal Evolution of Ultrafast Magnetic Field Generation in Molecules with Intense Bichromatic Circularly Polarized Uv Laser Pulses. Phys Rev A (2019) 99:053416. 10.1103/physreva.99.053416
58.
SederbergSKongFCorkumPB. Tesla-scale Terahertz Magnetic Impulses. Phys Rev X (2020) 10:011063. 10.1103/physrevx.10.011063
59.
BandraukA. D.ShenH.Exponential split operator methods for solving coupled time‐dependent Schrödinger equations. The Journal of Chemical Physics (1993) 99:1185, 1193. 10.1063/1.465362
60.
BandraukA. D.LuH.Exponential Propagators (Integrators) for the Time-dependent Schrödinger Equation. J Theor Comput Chem (2013) 12:1340001. 10.1142/s0219633613400014
61.
JefimenkoOD. Electricity and Magnetism: An Introduction to the Theory of Electric and Magnetic Fields. 2nd ed.Star City, West Virginia: Electret Scientific Co. (1989)
62.
YuanKJBandraukAD. Symmetry in Circularly Polarized Molecular High-Order Harmonic Generation with Intense Bicircular Laser Pulses. Phys Rev A (2018) 97:023408. 10.1103/physreva.97.023408
63.
YuanK.-J.BandraukA. D.Time-resolved Photoelectron Imaging of Molecular Coherent Excitation and Charge Migration by Ultrashort Laser Pulses. J Phys Chem A (2018) 122:2241, 2249. 10.1021/acs.jpca.7b11669
64.
YuanK.-J.ShuC.-C.DongD.BandraukA. D.Attosecond Dynamics of Molecular Electronic Ring Currents. J Phys Chem Lett (2017) 8:2229, 2235. 10.1021/acs.jpclett.7b00877
65.
SharpT. E.Potential-energy curves for molecular hydrogen and its ions. Atomic Data and Nuclear Data Tables (1970) 2:119, 169. 10.1016/s0092-640x(70)80007-9
66.
JuanKJGuoJBandraukAD. Circularly Polarized Attosecond Pulse Generation and Applications to Ultrafast Magnetism. J Opt (2017) 19:124016.
67.
YuanK.-J.BandraukA. D.Applied Sciences (2019) Probing Attosecond Electron Coherence in Molecular Charge Migration by Ultrafast X-Ray Photoelectron Imaging9:1941. 10.3390/app9091941
68.
CorkumP. B.BurnettN. H.BrunelF.Above-threshold Ionization in the Long-Wavelength Limit. Phys Rev Lett (1989) 62:1259, 1262. 10.1103/physrevlett.62.1259
69.
CorkumP. B.Plasma Perspective on strong Field Multiphoton Ionization. Phys Rev Lett (1993) 71:1994, 1997. 10.1103/physrevlett.71.1994
70.
YuanKJGuoJBandraukAD. Generation of Ultrafast Magnetic fields from Molecular Coherent Electron Currents. Phys Rev A (2018) 98:043410. 10.1103/physreva.98.043410
71.
XiaoX RWangM XLiangHGongQPengL-Y. Proposal for Measuring Electron Displacement Induced by a Short Laser Pulse. Phys Rev Lett (2019) 122:053201. 10.1103/physrevlett.122.053201
72.
KeldyshLV. Ionization in the Field of a strong Electromagnetic Wave. Sov Phys JETP (1965) 20:1307.
73.
MineoH.FujimuraY.Quantum Design of π-Electron Ring Currents in Polycyclic Aromatic Hydrocarbons: Parallel and Antiparallel Ring Currents in Naphthalene. J Phys Chem Lett (2017) 8:2019, 2025. 10.1021/acs.jpclett.7b00704
74.
ProninEAStaraceAFFrolovMVManakovNL. Perturbation Theory Analysis of Attosecond Photoionization. Phys Rev A (2009) 80:063403. 10.1103/physreva.80.063403
75.
ProninEAStaraceAFPengLY. Perturbation-theory Analysis of Ionization by a Chirped Few-Cycle Attosecond Pulse. Phys Rev A (2011) 84:013417. 10.1103/physreva.84.013417
76.
YuanK-JBandraukAD. Molecular Above-Threshold-Ionization Angular Distributions with Attosecond Bichromatic Intense XUV Laser Pulses. Phys Rev A (2012) 85:013413. 10.1103/physreva.85.013413
77.
BakosJ.KissA.SzabóL.TendlerM.Light Intensity Dependence of the Multiphoton Ionization Probability in the Resonance Case. Physics Letters A (1972) 41:163, 164. 10.1016/0375-9601(72)91095-x
78.
AgostiniP.BensoussanP.Resonant three‐photon ionization of potassium with circularly and linearly polarized light. Appl Phys Lett (1974) 24:216, 219. 10.1063/1.1655158
79.
YuanK.-J.BandraukA. D.Exploring Coherent Electron Excitation and Migration Dynamics by Electron Diffraction with Ultrashort X-ray Pulses. Phys Chem Chem Phys (2017) 19:25846, 25852. 10.1039/c7cp05067d
80.
YuanK.-J.BandraukA. D.Ultrafast X-ray Photoelectron Imaging of Attosecond Electron Dynamics in Molecular Coherent Excitation. J Phys Chem A (2019) 123:1328, 1336. 10.1021/acs.jpca.8b12313
81.
LeoneS. R.McCurdyC. W.BurgdörferJ.CederbaumL. S.ChangZ.DudovichN.et alWhat will it take to observe processes in 'real time'?Nature Photon (2014) 8:162, 166. 10.1038/nphoton.2014.48
82.
KrausPMMignoletBBaykushevaDRupenyanAHornýLPenkaEFet alMeasurement and Laser Control of Attosecond Charge Migration in Ionized Iodoacetylene. Science (2015) 350:790, 5. 10.1126/science.aab2160
83.
BandraukADChelkowskiSCorkumPBManzJYudinGL. Attosecond Photoinization of Coherent Superpositions of Bound-Dissociative Molecular States-Effects of Nuclear Motion. J Phys (2009) B42:143001.
84.
JiaD.ManzJ.YangY.De- and Recoherence of Charge Migration in Ionized Iodoacetylene. J Phys Chem Lett (2019) 10:4273, 4277. 10.1021/acs.jpclett.9b01687
85.
ParamonovGKKlamrothTLuHZBandraukAD. Quantum Dynamic Isotope Effects and Power Spectra of H2+ and HD+ in Strong One-Cycle Laser Pulses”. Phys Rev (2018) A98:063431.
86.
ZhuBJet al30 Co-authors), Strong Magnetic fields with Simple Coils Irradiated by High Power Lasers. Appl Phys Lett (2015) 107:2619003. 10.1063/1.4939119
Summary
Keywords
magnetic field generation, intense laser pulses, coherent ring currents, multiple ionization pathways, bichromatic circularly polarized pulse
Citation
Bandrauk AD, Chelkowski S and Yuan K-J (2021) Electronic Currents and Magnetic Fields in Induced by Coherent Resonant Bichromatic Circularly Polarized Laser Pulses: Effects of Orientation, Phase, and Helicity. Front. Phys. 9:675375. doi: 10.3389/fphy.2021.675375
Received
03 March 2021
Accepted
14 May 2021
Published
14 June 2021
Volume
9 - 2021
Edited by
Robert Gordon, University of Illinois at Chicago, United States
Reviewed by
Hirohiko Kono, Tohoku University, Japan
Joern Manz, Freie Universität Berlin, Germany
Updates
Copyright
© 2021 Bandrauk, Chelkowski and Yuan.
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: André D. Bandrauk, Andre.Dieter.Bandrauk@usherbrooke.ca
†Deceased
This article was submitted to Physical Chemistry and Chemical Physics, a section of the journal Frontiers in Physics
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.