Electronic Currents and Magnetic Fields in H 2 + Induced by Coherent Resonant Bichromatic Circularly Polarized Laser Pulses: Effects of Orientation, Phase, and Helicity

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.


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][4][5][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][8][9][10][11][12][13]. By creating unidirectional constant valence-type electronic currents in molecules with circularly polarized UV laser pulses, static magnetic fields [7][8][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 2p ± orbital in the hydrogen-like atom, which can be prepared via resonant 1s → 2p ± 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][17][18][19]. Linearly polarized laser pulses can also induce excited ring currents by controlling the rotation direction of π electrons in planar/nonplanar aromatic molecules [20][21][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][28][29][30][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][33][34][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][37][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][42][43] and molecular systems [44][45][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][50][51][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][54][55]. Numerical simulations are performed on the aligned molecular ion H + 2 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 H + 2 to simulate electron currents and magnetic field generation in Section 2. The results of ultrafast magnetic FIGURE 1 | Illustration of ultrafast magnetic field generation B under resonant excitation of H + 2 aligned along the x/z axis by bichromatic co-and counter-rotating circularly polarized XUV pulses E(t) with their field vectors polarized in the (x, y) plane. Two excitation processes for the molecule axis parallel and perpendicular to the laser (x, y) polarization plane, that is, (A) in plane R E and (B) around axis R ⊥ E, are compared. (C) Molecular σ g − π ± u resonant excitation and ionization are shown for different molecular alignments, and corresponding evolutions of coherent electron wave packet density distributions at different times. Protons are at ±R/2 ±1 a.u.
Frontiers in Physics | www.frontiersin.org June 2021 | Volume 9 | Article 675375 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 Z e 2 m e 1 are used unless otherwise stated. ± 1 denotes the helicity of combined fields, that is, corotating (ϵ +1) or counter-rotating (ϵ −1) components. ϕ 1 and ϕ 2 are CEPs of the pulses E 1 (t) and E 2 (t), respectively. A smooth sin 2 (πt/T lp ) pulse envelope f (t) for maximum amplitude E, intensity I I x I y cε 0 E 2 /2, and duration T lp nτ is adopted, where one optical cycle period τ 1,2 2π/ω 1,2 . The 3D TDSE in Eq. 1 is propagated by a second-order split operator method which conserves unitarity in each time step Δt combined with a fifth-order finite difference method and Fourier transform technique in the spatial steps Δρ, Δz, and Δθ [59,60]. The initial electron wave function ψ(r, t 0) is prepared in the ground 1sσ g 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 Δt 0.01 au 0.24 as. The spatial discretization is Δρ Δz 0.25 au for a radial grid range 0 ≤ ρ ≤ 128 au (6.77 nm) and |z| ≤ 32 au (1.69 nm), and the angle grid size Δθ 0.025 radian. To prevent unphysical effects due to the reflection of the wave packet from the boundary, we multiply ψ(ρ, θ, z, t) by a "mask function" or absorber in the radial coordinate ρ with the form cos 1/8 [π(ρ − ρ a )/2ρ abs ]. 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 α d E/ω 2 1/2 of the electron [25]. The time-dependent probability current density is defined by the quantum expression in the length gauge, where ψ(r, t) is the exact Born-Oppenheimer (static nuclei) electron wave function obtained from the TDSE in Eq. 1, and ∇ r e ρ ∇ ρ + e θ ∇ θ /ρ + e z ∇ z is the momentum operator in cylindrical coordinates. Then the corresponding timedependent magnetic field is calculated using the following classical Jefimenko equation [61]: where t r t − r/c is the retarded time and μ 0 4π × 10 − 7 NA −2 (6.692 × 10 −4 au) is the permeability of free space. Units of B(r, t) are teslas (1 T 10 4 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, B(r, t) −e(μ 0 /4π)(v × r)/r 3 [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 j electric −ej. This explains the sign (−) in Eq. 5, e 1 in atomic units.

RESULTS AND DISCUSSIONS
The ground and excited states of H + 2 ion and its electron potentials, 1sσ g , 1sσ u , 2sσ g , 2sσ g , 1π u , etc., are well documented in [65]. Numerical solutions of the TDSE for H + 2 aligned with the laser polarization (x-y) plane, that is, R to the x-axis with the electric field vector E(t) rotating in the (x,y) plane ( Figure 1A), are used for obtaining time-dependent probabilities P(t) of the 1sσ g and 2pπ u states at equilibrium R 2 au. excited by a 70-nm pulse at two different intensities, I 2× 10 14 W/cm 2 and I 1 × 10 15 I 1×10 15 W/cm 2 , by a five-cycle pulse (one cycle 0.234 fs 234 as). Thus, at I 2 × 10 14 , Figure 2A, one sees no significant ionization, whereas at I 1 × 10 15 W/cm 2 , Figure 2B, the 2pπ u 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 2pπ u state at 70 nm and some Rydberg states above the 2pπ u state shown in 1.
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 H + 2 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ω 1 ΔE and the energy difference ΔE E π u − E σ g 0.65 au, a resonant excitation with the ground 1sσ g state and the excited π u electronic states occurs. Moreover, with ω 2 2ω 1 , such bichromatic laser pulses can also produce CEWPs with the same kinetic energies, E e 2ω 1 − I p ω 2 − I p , 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 (x, y) 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 1sσ g state and the excited π u state by the λ 1 70 nm circularly polarized laser pulse is dependent on the molecular alignment. In the case of R E, the in-plane degenerates perpendicular excitation, π ± u ←σ g for different angular momenta Δm ± 1 dominates, and the electron density distribution evolves nearly perpendicular to the molecular axis [63], whereas in the case of R⊥E, the produced electron wave packets move around the molecular axis due to individual π +/− u ←σ g transitions Δm 1 or −1 [64]. Meanwhile, two ω 1 and single ω 2 transitions to the continuum (2ω 1 ω 2 > I p ) 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 (ϵ +1) and counter-rotating (ϵ −1) pulse combinations, are compared as well.

Dependence of Generated Magnetic Fields on Circular Error Probability and Molecular Alignments
We first present the counter-rotating (ϵ −1) dynamics. Figure 3 shows results of the generated maximum magnetic field B at the molecular center r 0 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 H + 2 aligned along a) the x axis, R E, and b) the z axis, R ⊥ E, 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 × 10 14 W/cm 2 (E 0.0534 au) and durations T lp 48.3 au 1.16 fs, corresponding to 580 as FWHM (full width at half maximum), that is, 5τ 1 for ω 1 pulse and 10τ 2 for ω 2 pulse. We also set the pulse CEP ϕ 1 0 and vary the CEP ϕ 2 from 0 to 2π, that is, the relative CEP is ϕ ϕ 2 − ϕ 1 ϕ 2 . Figure 3 shows the CEP ϕ dependence of the generated magnetic fields with different molecular alignments.

Resonant Excitation in the Case of In-Plane R E
For the case of in-plane excitation with the molecule axis parallel to the laser (x, y) 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 2π, B oscillates periodically as cos(ϕ). The maximum and minimum values of the magnetic field are B 0.225 T and 0.165 T at ϕ ≈ 0 (or 2π) 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, B ω 1 0.194 T and B ω 2 0.022 T. Due to the σ g − π + u resonant excitation, the magnetic field induced by the resonant 70 pulse is much stronger than that induced by the nonresonant 35-nm pulse, B ω1 ≈ 8.8B ω2 . The sum value of the generated magnetic fields B s B ω1 + B ω2 is slightly smaller than the 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 D 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, ψ D √ exp(ip · r/Z), where p is the electron momentum tangent to the ring, one gets after performing the ∇ derivation a simple relation j D p m e Dv as shown in [58]. Next, for the case of the electron current mainly localized in a plane one gets from Eq. 5, that is, B ∼ j/r 2 ∼ Dv/r 2 , we thus derived Eq. 6.
Thus, in Eq. 6, the magnetic field is proportional to the electron probability density D and the electron velocity which originates from the electron current density j Dv [58]. Density D 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 D 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 σ g state to the degenerate π ± u state with ω 1 E e − E g via a perpendicular transition. Meanwhile, the absorptions of two ω 1 photons and one ω 2 2ω 1 photon give rise to photoelectron wave packets with the same kinetic energies E e 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 D (1) and D (2) are insensitive to the pulse phases, whereas the interference term D (1,2) 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 D D (1) + D (2) + D (1,2) is also a cosine function of the relative pulse phase ϕ with the form D ∼ cos(ϕ). 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 B ∼ cos(ϕ), giving rise to a CEP dependence of the generated magnetic field, as illustrated in Figure 3A (blue diamond), B (r 0).

Resonant Excitation in the Case of Around-Axis R ⊥ E
For around R axis excitation, Figure 3B shows phase-dependent magnetic field B generation in the process of molecule axis perpendicular to the laser (x, y) polarization plane, leading to around-axis currents. Comparing to the in-plane R E 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 B 0.38 T. The magnetic fields generated by the individual bichromatic pulses are less than the sum of the two 0.019 T. The bichromatic magnetic fields mainly arise from the single-photon 70 nm process, B B ⊥ ω1 , whereas the contribution from the absorption of single 35 nm (ω 2 ) or two 70 nm (ω 1 ) 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 ω 1 /ω 2 pulse case, the generated magnetic fields are also dependent on the molecular alignments and the pulse wavelength. At λ 70 nm, B ⊥ ω 1 ≈ 2B ω 1 , the perpendicular case is stronger, whereas at 35 nm, the generated magnetic fields are nearly equivalent, B ⊥ ω2 ≈ B ω2 . The difference indicates essentially the importance of the around R axis ring electron currents in the resonant intermediate π u electronic state.
For the optical responses in Figure 3B of the molecule H + 2 perpendicular to the laser (x, y) polarization plane, R ⊥ E, Figure 1B, by a single 70-nm circularly polarized pulse, resonant excitation between the ground σ g state and the excited π + u state with magnetic quantum number Δm 1 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 j ⊥ θ (r, t) in the following form Eq. 8 combined with Eqs. 21 and 26 shows that for the superposition electron state ψ c (r, t), the time-dependent electronic density D(r, t) and the induced angular current j ⊥ θ (r, t) evolve in time with the electron coherence period of Δτ (0) 2π/ΔE. 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 ψ c (r, t) in Eq. 21 by the λ 70 nm pulse. The contributions from the continuum electron wave packets with energy E e by the 35-nm pulse are clearly very weak, that is, B ⊥ ω 1 ≈ 20B ⊥ ω 2 , and therefore negligible. As a result, altering the relative CEP ϕ does not lead to a variance of the generated magnetic field B.

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 j θ (r, t) with time t is determined by the coherent resonant excitation and the molecular alignments. Figure 4 illustrates snapshots of angular electron probability currents j ⊥ θ (x, y, z 0, t) at the center of H + 2 obtained from Eq. 4 at different times t (unit of τ 2π/ω 1 τ 1 ) for the molecule H + 2 perpendicular to the laser (x, y) polarization plane, R ⊥ E, Figure 1, leading to around R axis currents, j ⊥ θ (x, y, z 0, t) j ⊥ θ (ρ cos θ, ρ sin θ, z 0, t). The bichromatic counter-rotating (ϵ −1) circularly polarized XUV pulse at CEPs ϕ 1 ϕ 2 0 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 (x, y) 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 σ g state and the excited π + u state by the 70-nm circularly polarized pulse [64]. From such coherent electron currents, the generated magnetic fields are unidirectional, along the z/R axis. At the molecular center z 0, the magnetic field is the sum of those at the two nuclei, ± R/2. 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.
For comparison, in Figure 5, we also plot the in-plane electron probability current j θ (x, y, z 0, t) at the center of the molecule H + 2 aligned along the x axis, parallel to the laser polarization (x, y) plane, R E 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 σ g − π ± u 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 (p y ←s) transitions dominate during the excitation processes. In this case, the coherent exited electronic state is given by where π ± u (r) [π x u (r) ± iπ y u (r)]/ 2 √ and |π ± u (r)| 2 [ π x u (r)| 2 + π y u (r)| 2 ] / 2. Then the corresponding interference term in the time-dependent electron density becomes (in cylindrical coordinates (ρ, θ, z) as shown in [9,64] and in Supplementary Appendix A2) D ( e,g ) (r, t) 2c g c e ψ g (r)ψ ± e ρ, z cos(ΔEt)sin(θ), (10) and the time-dependent current j θ (r, t) −2c g c e ψ g (r)ψ ± e ρ, z ΔE sin(ΔEt)cos(θ), whereψ ± e (ρ, z) is the absolute value of ψ ± e (see Eq. 23). It is found that the coherent density and current superposition terms in Eqs. 10 and 11 follow the forms ∼ sin(θ) and cos(θ). As a result, the coherent wave packets due to the σ g − π ± u 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 ± R/2 has opposite phases [70]. Their overlap leads to a weak magnetic field at the molecule center z 0. The generated magnetic field in Figure 3A therefore mainly results from the coherent excitation by bichromatic circularly polarized laser pulses.

Influence of the Pulse Helicity
We next study the process with a bichromatic corotating (ϵ +1) 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 (ϵ −1) case in Figure 3. For the case of R E, the magnetic field B varies as a sine function of the phase ϕ, as predicted in Supplementary Appendix A1. A phase π/2 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 (ϵ +1) case the generated magnetic field depends on the molecular alignment. In Figure 6B, we see that for the case of R⊥E, 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 B 0.38 T is obtained. The σ g − π + u resonant excitation with Δm 1 gives rise to unidirection ring electron currents, as illustrated in Figure 3. The contribution from the coherent electron wave packets ψ 0 (r, t) in Eq. 21 dominates. The effects of the multiple pathway (2ω 1 and ω 2 ) 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.
It should be noted that in the case of R E in the general pulse E(t) with its field polarization vectors in the (x, y) plane, the strength of the generated magnetic field is slightly sensitive to the pulse helicity ϵ. As shown in Figure 6A, at ϵ +1 for corotating cases, the maximum magnetic field is B 0.206 T at ϕ 0.4π, which is slightly weaker than that at ϵ −1 for counter-rotating cases in Figure 3A, where B 0.225 T at ϕ 0. Eq. 6 predicts that the induced magnetic field, in general, is proportional to the ratio of the electron velocity v and inversely proportional to r 2 , 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 laserinduced dynamics, which depends on the helicity of driving pulses.
FIGURE 6 | Dependence of the maximum magnetic field B (blue diamond) at the molecular center, r 0 on the relative pulse CEP ϕ ϕ 2 − ϕ 1 for the aligned molecule H + 2 by bichromatic corotating (ϵ +1) circularly polarized laser pulses with their field vectors polarized in the (x, y) plane. Two cases are compared for the molecule aligned parallel to (A) the x axis and (B) the z axis, that is, parallel R E in Figure 1A and perpendicular R ⊥ E in Figure 1B Frontiers in Physics | www.frontiersin.org June 2021 | Volume 9 | Article 675375 8 where t 0 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, x(t 0 ) y(t 0 ) 0. 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 c > 1 [72]. The initial electron velocities are nearly equivalent with 2(ω 2 − I p ) . As shown in Figure 7, for the corotating ϵ 1 case, the radius of the continuum electron is much larger than that of the counter-rotating ϵ −1 case. The electron moves way quickly in the corotating field, whereas the counter-rotating field restricts the ionized electron around the molecular center. Because B ∼ j/r 2 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.

CONCLUSION
We present the ultrafast magnetic field generation in molecules from one electron molecular TDSE simulation under effects of coherent resonant excitation in bichromatic (ω 2 2ω 1 ) co-and counter-rotating circularly polarized laser fields. Numerical results are obtained for the aligned molecular ion H + 2 , which can be fully and exactly studied. We evaluate the generated magnetic field B(r 0) 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, B(r 0), which is shown to be dependent on the molecule alignment, and is maximum for in-plane bichromatic excitation.
In a bichromatic (frequency ω 1 ≠ ω 2 ) circularly polarized field, a σ g − π ± u resonant excitation of individual magnetic components m ± 1 is triggered by the resonant ω 1 (70 nm) pulse. As a result, coherent electron currents between the ground σ g state and the degenerate excited π ± u electronic state are induced in molecules. We compare two nonlinear responses for different molecular alignments: • For the around-axis case, R ⊥ E, the molecule aligned along the z axis, perpendicular to the laser (x, y) polarization plane, Figure 1B, and the magnetic field mainly results from the coherent electron currents which are induced by one resonant ω 1 photon around R. The contribution from two ω 1 photons or one ω 2 photon is negligible. Varying the relative CEP ϕ does not influence the generated magnetic field. Similar results are obtained for both counter-rotating ε −1, Figure 3B, and corotating ε +1, Figure 6B, excitations, confirming the main role of the around-axis circular coherent electron current. • For the in-plane case, R E, the molecule aligned parallel to the laser (x, y) 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 ω 1 photon absorption, that is, the σ g − π ± u resonant excitation, and two ω 1 or one ω 2 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 ω 1 and ω 2 photoionization pathways, altering their relative phase ϕ gives rise to a modulation of the generated magnetic field B(r 0) at the molecular center with forms ∼ cos(ϕ) for counter-rotating and counter-rotating (ϵ −1) helicity schemes. The initial position is x(t 0 ) y(t 0 ) 0, where t * 0 indicates the ionization time at the maxima of the combined fields. Green line corresponds to the minimum laser-induced radius r 0 of an ionized electron at frequency ω 1 and circular field amplitude E, r 0 (2E/ω 2 1 ) 1 + (π/2) 2 as shown in [23].
sin(ϕ) 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 H + 2 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.

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.