SEARCH · Search NASA

Results for “fields”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field

Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. Here, in this study, we motivate and explore the role of the volume-averaged root-mean-squared (rms) magnetic coupling term between the turbulent, δB , and large-scale, B 0 , fields, $\langle$(δB·B 0 ) 2 $\rangle$$_{V}^{1/2}$. By considering the second moments of the energy balance equations we show that the rms coupling term is in energy equipartition with the volume-averaged turbulent kinetic energy for turbulence with a sub-Alfvénic large-scale field. Under the assumption of exact energy equipartition between these terms, we derive relations for the magnetic and coupling term fluctuations, which provide excellent, parameter-free agreement with time-averaged data from 280 numerical simulations of compressible magnetohydrodynamic (MHD) turbulence. Furthermore, we explore the relation between the turbulent mean field and total Alfvén Mach numbers, and demonstrate that sub-Alfvénic turbulence can only be developed through a strong, large-scale magnetic field, which supports an extremely super-Alfvénic turbulent magnetic field. This means that the magnetic field fluctuations are significantly subdominant to the velocity fluctuations in the sub-Alfvénic large-scale field regime. Throughout our study, we broadly discuss the implications for observations of magnetic fields and understanding the dynamics in the magnetized ISM.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Ionization of HCCI Neutral and Cations by Strong Laser Fields Simulated With Time Dependent Configuration Interaction

Strong field ionization of neutral iodoacetylene (HCCI) can produce a coherent superposition of the X and A cations. This superposition results in charge migration between the CC π orbital and the iodine π-type lone pair which can be monitored by strong field ionization with short, intense probe pulses. Strong field ionization of the X and A states of HCCI cation was simulated with time-dependent configuration interaction using singly ionized configurations and singly excited, singly ionized configurations (TD-CISD-IP) and an absorbing boundary. Studies with static fields were used to obtain the 3-dimensional angular dependence of instantaneous ionization rates by strong fields and the orbitals involved in producing the cations and dications. The frequency of charge oscillation is determined by the energy separation of the X and A states; this separation can change depending on the direction and strength of the field. Furthermore, fields along the molecular axis can cause extensive mixing between the field-free X and A configurations. For coherent superpositions of the X and A states, the charge oscillations are characterized by two frequencies–the driving frequency of the laser field of the probe pulse and the intrinsic frequency due to the energy separation between the X and A states. For linear and circularly polarized pulses, the ionization rates show marked differences that depend on the polarization direction of the pulse, the carrier envelope phase and initial phase of the superposition. Varying the initial phase of the superposition at the beginning of the probe pulse is analogous to changing the delay between the pump and probe pulses. The charge oscillation in the coherent superposition of the X and A states results in maxima and minima in the ionization yield as a function of the superposition phase

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

DeepM&Mnet: Inferring the electroconvection multiphysics fields based on operator approximation by neural networks

Electroconvection is a multiphysics problem involving coupling of the flow field with the electric field as well as the cation and anion concentration fields. Here, we use electroconvection as a benchmark problem to put forward a new data assimilation framework, the DeepM&Mnet, for simulating multiphysics and multiscale problems at speeds much faster than standard numerical methods using pre-trained neural networks. We first pre-train DeepONets that can predict independently each field, given general inputs from the rest of the fields of the coupled system. DeepONets can approximate nonlinear operators and are composed of two sub-networks, a branch net for the input fields and a trunk net for the locations of the output field. DeepONets, which are extremely fast, are used as building blocks in the DeepM&Mnet and form constraints for the multiphysics solution along with some sparse available measurements of any of the fields. We demonstrate the new methodology and document the accuracy of each individual DeepONet, and subsequently we present two different DeepM&Mnet architectures that infer accurately and efficiently 2D electroconvection fields for unseen electric potentials. Furthermore, the DeepM&Mnet framework is general and can be applied for building any complex multiphysics and multiscale models based on very few measurements using pre-trained DeepONets in a “plug-and-play” mode.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

The effects of field history on magnetic skyrmion formation in [Pt/Co/Ir] 3 multilayers

Here, we report on the formation of Néel-type magnetic bubble skyrmions at room temperature in [Pt/Co/Ir] 3 multilayered thin films after an in-plane magnetic field treatment. Polar magneto-optical Kerr Effect (p-MOKE) microscopy images show that the dendritic magnetic configurations observed after AC demagnetization evolve into magnetic bubble skyrmions after the application and subsequent removal of an in-plane magnetic field. Micromagnetic simulations were used to systematically investigate the role of the in-plane magnetic field magnitude, misalignment of the sample, and the Dzyaloshinskii-Moriya interactions (DMI) in generating bubble skyrmions during the field treatment. The simulations show that in-plane fields slightly below the saturation field are the most effective at producing skyrmions, and, furthermore, a small field angle away from the sample plane not only leads to improved skyrmion formation but also provides a means to select the skyrmion polarity where the direction of the out-of-plane component of the field is opposite to the direction of the skyrmion cores. This field treatment scheme leads to a simple and reliable way to create magnetic bubble skyrmions in multilayered thin films with DMI.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Field-resilient supercurrent diode in a multiferroic Josephson junction

The research on supercurrent diodes has surged rapidly due to their potential applications in electronic circuits at cryogenic temperatures. To unlock this functionality, it is essential to find supercurrent diodes that can work consistently at zero magnetic field and under ubiquitous stray fields generated in electronic circuits. However, a supercurrent diode with robust field tolerance is currently lacking. Here, we demonstrate a field-resilient supercurrent diode by incorporating a 2D multiferroic material into a Josephson junction, and observed a pronounced supercurrent diode effect at zero magnetic field. More importantly, the supercurrent rectification persists over a wide and bipolar magnetic field range beyond industrial standards for field tolerance. By theoretically modeling a multiferroic Josephson junction, we unveil that the interplay between spin-orbit coupling and multiferroicity underlies the unusual field resilience of the observed diode effect. This work introduces multiferroic Josephson junctions as a new field-resilient superconducting device for cryogenic electronics.

Yang, Hung-Yu [Univ. of California, Los Angeles, C↗

DOE OSTI2021

Accelerating phase-field-based microstructure evolution predictions via surrogate models trained by machine learning methods

Abstract The phase-field method is a powerful and versatile computational approach for modeling the evolution of microstructures and associated properties for a wide variety of physical, chemical, and biological systems. However, existing high-fidelity phase-field models are inherently computationally expensive, requiring high-performance computing resources and sophisticated numerical integration schemes to achieve a useful degree of accuracy. In this paper, we present a computationally inexpensive, accurate, data-driven surrogate model that directly learns the microstructural evolution of targeted systems by combining phase-field and history-dependent machine-learning techniques. We integrate a statistically representative, low-dimensional description of the microstructure, obtained directly from phase-field simulations, with either a time-series multivariate adaptive regression splines autoregressive algorithm or a long short-term memory neural network. The neural-network-trained surrogate model shows the best performance and accurately predicts the nonlinear microstructure evolution of a two-phase mixture during spinodal decomposition in seconds, without the need for “on-the-fly” solutions of the phase-field equations of motion. We also show that the predictions from our machine-learned surrogate model can be fed directly as an input into a classical high-fidelity phase-field model in order to accelerate the high-fidelity phase-field simulations by leaping in time. Such machine-learned phase-field framework opens a promising path forward to use accelerated phase-field simulations for discovering, understanding, and predicting processing–microstructure–performance relationships.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Spatial patterns of historical crop yields reveal soil health attributes in US Midwest fields

Abstract Attaining high crop yields and increasing carbon storage in agricultural soils, while avoiding negative environmental impacts on water quality, soil erosion, and biodiversity, requires accurate and precise management of crop inputs and management practices. The long-term analysis of spatial and temporal patterns of crop yields provides insights on how yields vary in a field, with parts of field constantly producing either high yields or low yields and other parts that fluctuate from one year to the next. The concept of yield stability has shown to be informative on how plants translate the effects of environmental conditions (e.g., soil, climate, topography) across the field and over the years in the final yield, and as a valuable layer in developing prescription maps of variable fertilizer rate inputs. Using known relationships between soil health and crop yields, we hypothesize that areas with measured constantly low yield will return low carbon to the soil affecting its heath. On this premises, yield stability zones (YSZ) provide an effective and practical integrative measure of the small-scale variability of soil health on a field relative basis. We tested this hypothesis by measuring various metrics of soil health from commercial farmers’ fields in the north central Midwest of the USA in samples replicated across YSZ, using a soil test suite commonly used by producers and stakeholders active in agricultural carbon credits markets. We found that the use of YSZ allowed us to successfully partition field-relative soil organic carbon (SOC) and soil health metrics into statistically distinct regions. Low and stable (LS) yield zones were statistically lower in normalized SOC when compared to high and stable (HS) and unstable (US) yield zones. The drivers of the yield differences within a field are a series of factors ranging from climate, topography and soil. LS zones occur in areas of compacted soil layers or shallow soils (edge of the field) on steeper slopes. The US zones occurring with high water flow accumulation, were more dependent on topography and rainfall. The differences in the components of the overall soil health score (SHS) between these YSZ increased with sample depth suggesting a deeper topsoil in the US and HS zones, driven by the accumulation of water, nutrients, and carbon downslope. Comparison of the field management provided initial evidence that zero tillage reduces the magnitude of the variance in SOC and soil health metrics between the YSZ.

Fowler, Ames↗

DOE OSTI2021

Guide field effects on the distribution of plasmoids in multiple scale reconnection

The effects of a finite guide field on the distribution of plasmoids in high-Lundquist-number current sheets undergoing magnetic reconnection in large plasmas are investigated with statistical models. Merging of plasmoids is taken into account either assuming that guide field flux is conserved resulting in nonforce-free profiles in general, or that magnetic helicity is conserved and Taylor relaxation occurs to convert part of the summed guide field flux into reconnecting field flux toward minimum energy states resulting in force-free profiles. It is found that the plasmoid distribution in terms of reconnecting field flux follows a power law with index 7/4 or 1 depending on whether merger frequencies are independent of or dependent on their relative velocity to the outflow speed, respectively. This result is approximately the same for the force-free and nonforce-free models, with nonforce-free models exhibiting indices of 2 and 1 for the same velocity dependencies. Distributions in terms of guide field flux yield indices of 3/2 for the nonforce-free model regardless of velocity dependence. This is notably distinct from the indices of 11/8 and 1 for the force-free models independent of and dependent on velocity, respectively. At low guide field fluxes, the force-free models exhibit a second power law index of 1/2 due to nonconstant flux growth rates. The velocity-dependent force-free model predicts the production of slightly more rapidly moving large guide field flux plasmoids which are supported by observational evidence of flux ropes with strong core fields. Implications are discussed on particle acceleration via Fermi processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Near-field spectroscopy of Dirac plasmons in Bi 2 Se 3 ribbon arrays

Plasmons supported in the massless electron surface states of topological insulators (TIs), known as Dirac plasmons, have great potential in next generation optoelectronics. However, their inherent confinement to the surface makes the investigation of Dirac plasmons challenging. Near-field techniques provide the ideal platform to directly probe Dirac plasmons due to the sensitivity to evanescent fields at the surface. Here, we demonstrate the use of aperture near-field spectroscopy for the investigation of localized terahertz (THz) Dirac plasmon resonances in Bi 2 Se 3 ribbon arrays with widths ranging from 10 to 40 µm. Unlike scattering THz near-field techniques, the aperture method is most sensitive to plasmons with the relevant lower-momenta corresponding to plasmon wavelengths on the scale of ~20 µm. The combination of THz time-domain spectroscopy and aperture near-field microscopy enables sampling of localized Dirac plasmons in the near-field zone in the 0.5–2.5 THz range. We map the plasmon dispersion, which reveals a coupled plasmon–phonon polariton interaction. The near-field spectra show a higher contrast of the upper polariton branch in comparison with far-field observations. The information revealed by aperture near-field spectroscopy could deepen our understanding of the behavior of Dirac plasmons, leading to the potential development of real-world TI devices.

42 ENGINEERING↗

DOE OSTI2022

Impact of CH 4 addition on the electron properties and electric field dynamics in a Ar nanosecond-pulsed dielectric barrier discharge

Non-equilibrium plasmas derive their low temperature reactivity from producing and driving energetic electrons and active species under large electric fields. Therefore, the impact of reactants on the plasma properties including electron number density, electric field, and electron temperature is critical for applications such as plasma CH 4 reforming. Due to experimental complexity, electron properties and the electric field are rarely measured together in the same discharge. Here, we combine time-resolved Thomson scattering and electric field induced second harmonic generation (EFISH) to probe electron temperature, electron density, and electric field strength in a 60 Torr CH 4 /Ar nanosecond-pulsed dielectric barrier discharge (ns-DBD) while varying the CH 4 mole fraction from 0 to 8%. These measurements are compared to a 1-D numerical model to benchmark its predictions and identify areas of uncertainty. Nonlinear coupling between CH 4 addition, electron temperature, electron density, and the electric field was directly observed. Contrary to previous measurements in He, the electron temperature increased with CH 4 mole fraction. This rise in electron temperature is identified as electron heating by residual electric fields that increased with larger CH 4 mole fraction. Moreover, the electron number density has been found to decrease rapidly with the increase of methane mole fraction. Comparison of these measurements with the model yielded better agreement at higher CH 4 mole fractions and with the usage of ab initio calculated Ar electron-impact cross-sections from the B-spline R-matrix (BSR) database. Furthermore, the calculated plasma properties are shown to be sensitive to the residual surface charge implanted on the quartz dielectric surfaces. Without considering surface charge in the simulations, the calculated electric field profiles agreed well with the measurements, but the electron properties were underpredicted by more than a factor of three. Therefore, measurements of either the electric field or electron properties measurements alone are insufficient to fully validate modelling predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Optimizing superconducting Nb film cavities by mitigating medium-field Q -slope through annealing

Niobium films are of interest in applications in various superconducting devices, such as superconducting radiofrequency cavities for particle accelerators and superconducting qubits for quantum computing. In this study, we address the persistent medium-field Q-slope issue in Nb film cavities, which, despite their high-quality factor at low RF fields, exhibit a significant Q-slope at medium RF fields compared to bulk Nb cavities. Traditional heat treatments, effective in reducing surface resistance and mitigating the Q-slope in bulk Nb cavities, are challenging for Nb-coated copper cavities. To overcome this challenge, we employed DC bias high-power impulse magnetron sputtering to deposit Nb film onto a 1.3 GHz single-cell elliptical bulk Nb cavity, followed by annealing treatments aimed at modifying the properties of the Nb film. In-situ annealing at 340 °C increased the quench field from 10.0 to 12.5 MV m −1 . Vacuum furnace annealing at 600 °C and 800 °C for 3 h resulted in a quench field increase of 13.5 and 15.3 MV m −1 , respectively. Further annealing at 800 °C for 6 h boosted the quench field to 17.5 MV m −1 . Additionally, the annealing treatments significantly reduced the field dependence of the surface resistance. However, increasing the annealing temperature to 900 °C induced a Q-switch phenomenon in the cavity. The analysis of RF performance and material characterization before and after annealing has provided critical insights into how the microstructure and impurity levels in Nb films influence the evolution of the Q-slope in Nb film cavities. Our findings highlight the significant roles of hydrides, high local misorientation, and lattice and surface defects in driving field-dependent losses. By strategically optimizing film properties and controlling impurity levels, we demonstrate a promising pathway to mitigate the medium-field Q-slope, paving the way for more efficient superconducting RF technologies.

Nb film↗

DOE OSTI2025

Magnetized ICF implosions: non-axial magnetic field topologies

This paper explores four different magnetic field topologies for application to spherical inertial confinement fusion implosions: axial, mirror, cusp and closed field lines. A mirror field is found to enhance the impact of magnetization over an axial field; this is because the mirror field more closely follows the hot-spot surface. A cusp field, while simple to generate, is not found to have any benefits over the tried-and-tested axial field. Closed field lines are found to be of the greatest benefit to hot-spot performance, with the simulated design undergoing a 2× increase in ion temperature before α-heating is considered. The plasma properties of the simulation with closed field lines are radically different from the unmagnetized counterpart, with electron temperatures in excess of 100 keV, suggesting that a fundamental redesign of the capsule implosion is possible if this method is pursued.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗

DOE OSTI2022

Magnetic fields in the formation of the first stars – II. Results

ABSTRACT Beginning with cosmological initial conditions at z = 100, we simulate the effects of magnetic fields on the formation of Population III stars and compare our results with the predictions of Paper I. We use gadget-2 to follow the evolution of the system while the field is weak. We introduce a new method for treating kinematic fields by tracking the evolution of the deformation tensor. The growth rate in this stage of the simulation is lower than expected for diffuse astrophysical plasmas, which have a very low resistivity (high magnetic Prandtl number); we attribute this to the large numerical resistivity in simulations, corresponding to a magnetic Prandtl number of order unity. When the magnetic field begins to be dynamically significant in the core of the minihalo at z = 27, we map it on to a uniform grid and follow the evolution in an adaptive mesh refinement, MHD simulation in orion2. The non-linear evolution of the field in the orion2 simulation violates flux-freezing and is consistent with the theory proposed by Xu & Lazarian. The fields approach equipartition with kinetic energy at densities ∼1010–1012 cm−3. When the same calculation is carried out in orion2 with no magnetic fields, several protostars form, ranging in mass from ∼1 to 30 M⊙; with magnetic fields, only a single ∼30 M⊙ protostar forms by the end of the simulation. Magnetic fields thus suppress the formation of low-mass Pop III stars, yielding a top-heavy Pop III initial mass function and contributing to the absence of observed Pop III stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Nested solitons in two-field fuzzy dark matter

ABSTRACT Dark matter as scalar particles consisting of multiple species is well motivated in string theory where axion fields are ubiquitous. A two-field fuzzy dark matter (FDM) model features two species of ultralight axion particles with different masses, m1 ≠ m2, which is extended from the standard one-field model with $m_a \sim 10^{-22} \, {\rm eV}$. Here we perform numerical simulations to explore the properties of two-field FDM haloes. We find that the central soliton has a nested structure when m2 ≫ m1, which is distinguishable from the generic flat-core soliton in one-field haloes. However, the formation of this nested soliton is subject to many factors, including the density fraction and mass ratio of the two fields. Finally, we study non-linear structure formation in two-field cosmological simulations with self-consistent initial conditions and find that the small-scale structure in two-field cosmology is also distinct from the one-field model in terms of DM halo counts and soliton formation time.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Comptonization by reconnection plasmoids in black hole coronae – III. Dependence on the guide field in pair plasma

ABSTRACT We perform non-radiative two-dimensional particle-in-cell simulations of magnetic reconnection for various strengths of the guide field (perpendicular to the reversing field), in magnetically dominated electron–positron plasmas. Magnetic reconnection under such conditions could operate in accretion disc coronae around black holes. There, it has been suggested that the transrelativistic bulk motions of reconnection plasmoids containing inverse-Compton-cooled electrons could Compton-upscatter soft photons to produce the observed non-thermal hard X-rays. Our simulations are performed for magnetizations 3 ≤ σ ≤ 40 (defined as the ratio of enthalpy density of the reversing field to plasma enthalpy density) and guide field strengths 0 ≤ Bg/B0 ≤ 1 (normalized to the reversing field strength B0). We find that the mean bulk energy of the reconnected plasma depends only weakly on the flow magnetization but strongly on the guide field strength – with Bg/B0 = 1 yielding a mean bulk energy twice smaller than Bg/B0 = 0. Similarly, the dispersion of bulk motions around the mean – a signature of stochasticity in the plasmoid chain’s motions – is weakly dependent on magnetization (for σ ≳ 10) but strongly dependent on the guide field strength – dropping by more than a factor of two from Bg/B0 = 0 to Bg/B0 = 1. In short, reconnection in strong guide fields (Bg/B0 ∼ 1) leads to slower and more ordered plasmoid bulk motions than its weak guide field (Bg/B0 ∼ 0) counterpart.

Astronomy & Astrophysics↗

DOE OSTI2024

Magnetorotational dynamo can generate large-scale vertical magnetic fields in 3D GRMHD simulations of accreting black holes

ABSTRACT Jetted astrophysical phenomena with black hole engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed three-dimensional global general relativistic magnetohydrodynamical simulations of accretion discs to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the non-linear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their non-linear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disc. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disc height. We characterize the complete dynamo mechanism with two time-scales: one for the local magnetic field generation, $t_{\rm gen}$, and one for the large-scale scale advection, $t_{\rm adv}$. Whereas the dynamo we describe is non-linear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $\alpha$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $\alpha$-dynamos.

Astronomy & Astrophysics↗

DOE OSTI2023

Magnetic bulk photovoltaic effect: Strong and weak field

Shift current and ballistic current have been proposed to explain the bulk photovoltaic effect (BPVE), and there have been experiments designed to separate the two mechanisms. These experiments are based on the assumption that under magnetic field, ballistic current can have a Hall effect while the shift current cannot, which is from some energy-scale arguments and has never been proven. A recent work [D. Hornung and R. von Baltz, Phys. Rev. B 103, 195203 (2021)] using quantum transport formalism achieves a conclusion that shift current indeed has a Hall current, seemingly contradicting the previous assumption and making the situation more confusing. Moreover, the behavior of BPVE under strong magnetic field is still unexplored. In this Letter, using a minimal two-dimensional tight-binding model, we carry out a systematic numerical study of the BPVE under weak and strong magnetic field by treating the field in a nonperturbative way. Further, our model clearly shows the appearance of the magnetically induced ballistic current along the transverse direction, which agrees with the previous predictions, and interestingly a sizable longitudinal response of the shift current is also observed, a phenomenon that is not captured by any existing theories where the magnetic field is treated perturbatively. More surprisingly, drastically different shift current is found in the strong-field regime, and the evolution from weak to strong field resembles a phase transition. We hope that our work could resolve the debate over the behavior of BPVE under magnetic field, and the strong-field behavior of shift current is expected to inspire more studies on the relation between nonlinear optics and quantum geometry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗