SEARCH · Engineering Papers

Results for “weak scaling”

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 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Identifying synoptic controls on boundary layer thermodynamic and cloud properties in a regional forecast model

Although most of our understanding of boundary layer cloudiness is based on idealized, subtropical, barotropic marine environments, boundary layer clouds exist across a range of conditions. In this study, we use the Naval Research Laboratory's Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS) and an automated cold-front-relative analysis framework to explore the boundary layer structure associated with low clouds across a transect through the cold front of a midlatitude synoptic cyclone. The model credibly captures boundary layer structure in line with conceptual models of midlatitude cyclones and ground-based observations at Graciosa Island in the Azores. The warm sector is conditionally unstable, with clouds that are too shallow and with too little liquid water, compared to cloud property retrievals from satellite and surface-based instruments. The cold-frontal region exhibits convection associated with weak stability and ascent. Northwest of the cold front, the boundary layer is well mixed, deeper, and capped by a strong inversion maintained by large-scale subsidence. Simulated clouds in frontal and post-frontal regions are mostly too thick, with too much liquid water and too little cloud base drizzle. The post-frontal clouds are associated with grid-scale updrafts, which appear to be the model's attempt to represent mesoscale organization of cellular convection typically observed in the cold sector of midlatitude cyclones. The deep, well-mixed post-frontal boundary layers and cloudiness are maintained by strong surface fluxes, as in cold air outbreaks. Our analysis framework serves as a unique approach to model verification, and our results offer insights into differences in boundary layer cloud behavior between subtropical and synoptic cold-sector regimes.

Eissner, Jordan M. [Univ. of Kansas, Lawrence, KS ↗

DOE OSTI2026

Constraining the Stellar-to-Halo Mass Relation with Galaxy Clustering and Weak Lensing from DES Year 3 Data

We develop a framework to study the relation between the stellar mass of a galaxy and the total mass of its host dark matter halo using galaxy clustering and galaxy-galaxy lensing measurements. We model a wide range of scales, roughly from $\sim 100 \; {\rm kpc}$ to $\sim 100 \; {\rm Mpc}$, using a theoretical framework based on the Halo Occupation Distribution and data from Year 3 of the Dark Energy Survey (DES) dataset. The new advances of this work include: 1) the generation and validation of a new stellar mass-selected galaxy sample in the range of $\log M_\star/M_\odot \sim 9.6$ to $\sim 11.5$; 2) the joint-modeling framework of galaxy clustering and galaxy-galaxy lensing that is able to describe our stellar mass-selected sample deep into the 1-halo regime; and 3) stellar-to-halo mass relation (SHMR) constraints from this dataset. In general, our SHMR constraints agree well with existing literature with various weak lensing measurements. We constrain the free parameters in the SHMR functional form $\log M_\star (M_h) = \log(εM_1) + f\left[ \log\left( M_h / M_1 \right) \right] - f(0)$, with $f(x) \equiv -\log(10^{αx}+1) + δ[\log(1+\exp(x))]^γ/ [1+\exp(10^{-x})]$, to be $\log M_1 = 11.506^{+0.325}_{-0.404}$, $\log ε= -1.632^{+0.306}_{-0.181}$, $α= -1.638^{+0.108}_{-0.099}$, $γ= 0.596^{+0.251}_{-0.210}$ and $δ= 3.810^{+2.045}_{-1.811}$. The inferred average satellite fraction is within $\sim 5-35\%$ for our fiducial results and we do not see any clear trends with redshift or stellar mass. Furthermore, we find that the inferred average galaxy bias values follow the generally expected trends with stellar mass and redshift. Our study is the first SHMR in DES in this mass range, and we expect the stellar mass sample to be of general interest for other science cases.

Zacharegkas, G. [Argonne; Chicago U., KICP] (ORCID↗

DOE OSTI2025

QUOKKA-based understanding of outflows (QED) – III. Outflow loading and phase structure as a function of galactic environment

ABSTRACT We present results from a suite of 3D high-resolution hydrodynamic simulations of supernova-driven outflows from galactic disc regions with a range of gas surface density, metallicity, and supernova scale height. We use this suite to quantify how outflow properties – particularly the loading factors for mass, metallicity, and energy – vary with these parameters. We find that the winds fall into three broad categories: steady and hot, multiphase and moderately bursty, and cool and highly bursty. The first of these is characterized by efficient metal and energy loading but weak mass loading, the second by moderate loading of mass, metals, and energy, and the third by negligible metal and energy loading but substantial mass loading. The most important factor in determining the kind of wind a galaxy will produce is the ratio of supernova to gas scale heights, with the latter set by a combination of supernova rate, metallicity-dependent cooling rate, and the gravitational potential. These often combine in counterintuitive ways – for example increased cooling causes cold clouds to sink into the galactic midplane more rapidly, lowering the volume-filling factor of dense gas and making the environment more favourable for strong winds. Our findings suggest that the nature of galactic winds is likely highly sensitive to phenomena such as runaway stars occurring at a large height and dense gas and are poorly captured in most simulations, and that metal loading factors for type Ia supernovae may be substantially larger than those for type II, with important implications for galactic chemical evolution.

Vijayan, Aditi (ORCID:0000000277142379)↗

DOE OSTI2022

Effect of interfacial vibrational coupling on surface wettability and water transport

Here, we report that the atomic-scale vibrational coupling at the solid-fluid interface can substantially alter the interfacial properties such as wettability and fluid slip. The wettability of water droplets on substrates subjected to various vibrational frequencies is studied using molecular dynamics simulation. The contact angle increases (i.e., becomes more hydrophobic) when the oscillation frequency of the substrate matches the intermolecular bending frequency of liquid water. We investigate the underlying mechanism by examining the dynamics of water molecules at the interface and find that the temporal contact between the solid and fluid is shorter when the frequencies match, resulting in weak solid-fluid adsorption. We further report that the vibrational match at the interface reduces wall-fluid friction and enhances water transport through the nanopore. Our findings demonstrate the importance of the atomic-scale vibrational coupling at the solid-fluid interface on the physicochemical behavior of nanodevices and biological nanochannels.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Exploring the Landscape of Distributed Graph Clustering on Leadership Supercomputers

The rapid growth of large-scale datasets in fields like biology and social networks has driven the need for advanced graph analytics techniques. Community detection, a fundamental task in graph analytics, identifies closely connected groups of nodes within a network, providing valuable insights across various disciplines. This study focuses on two classic community detection methods, the Louvain algorithm and Markov Clustering (MCL), and evaluates the performance of two prominent distributed community detection algorithms: HiPDPL-GPU, our prior implementation, and HipMCL. We conduct experiments on GPU-accelerated heterogeneous HPC systems, Summit and Frontier, to assess their performance under varying conditions. Our objective is to identify the strengths and weaknesses of these algorithms in terms of scalability, and quality of solutions. We evaluate these algorithms on a diverse set of 70+ networks spanning 13 domains, with sizes ranging up to 4.2 billion edges. Our results demonstrate that HiPDPL-GPU consistently outperforms HipMCL, especially for large-scale networks. HiPDPL-GPU achieves significantly faster runtimes (47x to 1439x), higher modularity scores, and improved scalability. These findings highlight HiPDPL-GPU as a promising solution for efficient and effective large-scale graph analytics in diverse application domains, and provide insights into the feasibility of using MCL-based approaches for certain application domains.

Community detection, graph algorithms↗

DOE OSTI2026

Atacama Cosmology Telescope: B-mode delensing with DR6 data and external tracers of large-scale structure

Large-scale -mode polarization of the cosmic microwave background (CMB) is a prime target for current and future experiments in search of primordial gravitational waves. With increasingly sensitive instruments being deployed, secondary -modes induced by the weak gravitational lensing of CMB photons are becoming an important limitation and need to be removed, a process known as delensing. In this work, we combine internally reconstructed CMB lensing maps from the Atacama Cosmology Telescope (ACT) data release 6 (DR6) with galaxy samples from unWISE and a map of the cosmic infrared background (CIB) fluctuations from Planck to produce a well-correlated tracer of the CMB lensing field. Our coadded tracer, shown to be 55%–85% correlated with the true lensing convergence at multipoles , is then convolved with ACT DR6 -mode polarization to yield a template of the lensing -modes. We assess its performance on a wide range of scales by using it to delens ACT DR6 and Planck -modes over 23% of the sky, removing around 39% of the lensing power at and 47% at , respectively. Our template achieves the highest delensing efficiency to date and will be useful for the analysis of early polarization maps from the Simons Observatory. We finally outline prospects for further improvements by including additional large-scale structure tracers from upcoming cosmological surveys.

Hertig, Emilie↗

DOE OSTI2022

THE DESIGN OF LAYER TIME OPTIMIZATION IN LARGE SCALE ADDITIVE MANUFACTURING WITH FIBER REINFORCED POLYMER COMPOSITES

In this study, we have developed a method to optimize a layer deposition time (a.k.a. layer time) for large-scale additive manufacturing (AM) via physics-based simulations. A long layer time leads to an over-cooled surface on which a new layer is deposited, and therefore, it may result in a weak bonding or debonding between layers, cracking, or warping. A short layer time leads to a high temperature of the structure due to insufficient cooling, and therefore, the structure may not be stiff enough and may collapse during manufacturing. Therefore, it is important to estimate the optimal layer time in additive manufacturing for a high-quality product. The temperature of a top layer right before deposition is recommended to be slightly higher than the glass temperature of the material. A temperature cooling was approximated to an exponential function of time, and the optimized layer time was obtained based on a target temperature while maintaining a minimal printing time. The material used is carbon fiber-reinforced polycarbonate (CF/PC), and the large-scale deposition system used is LSAM TM from Thermwood Corporation. Three different layer time cases were used for experiments, and a series of thermal images were obtained via an infra-red (IR) camera during the entire AM processes. AM process simulations were performed using a finite element method and the temperature profiles from the simulation were in good agreement with those from experiments. The layer time optimization was performed based on the temperature profiles from the simulations. A layer temperature with the optimal layer time was confirmed as the target temperature through simulation.

Jo, Eonyeon↗

DOE OSTI2021

Probing large-scale UV background inhomogeneity associated with quasars using metal absorption

ABSTRACT We study large-scale UV background inhomogeneity in three-dimensions associated with the observed quasar populations at high redshift. We do this by measuring metal absorption through quasar absorption spectrum stacking as a function distance to closest quasar in SDSS-IV/eBOSS on 10 s of comoving megaparsec scales. We study both intergalactic medium absorbers and mixed circumgalactic medium absorbers and probe absorption in O vi, C iv, and Si iv, and Si iii. Overall, stronger high ionization species absorption is seen closer to quasars at 2.4 < z < 3.1. O vi absorption shows a particularly strong change, with effects in C iv evident in some cases, and more marginal effects in Si iii and Si iv. We further study 2.05 < z < 2.4 (with weak signs of increasing homogeneity with time) and explore the study of metal absorption as a function of integrated SDSS-r band flux quasar flux (yielding consistent but less significant results). While the metal absorption does show sensitivity to large-scale 3D quasar proximity, the current incomplete quasar samples limit detailed interpretation. This work does, however, demonstrate that UV background inhomogeneities exist on scales of several 10 s of comoving megaparsecs associated with quasars and that they can be measured with precision by examining metal absorption in the intergalactic medium.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Multimodel validation of single wakes in neutral and stratified atmospheric conditions

Previous research has revealed the need for a validation study that considers several wake quantities and code types so that decisions on the trade-off between accuracy and computational cost can be well informed and appropriate to the intended application. In addition to guiding code choice and setup, rigorous model validation exercises are needed to identify weaknesses and strengths of specific models and guide future improvements. Here, we consider 13 approaches to simulating wakes observed with a nacelle-mounted lidar at the Scaled Wind Technology Facility (SWiFT) under varying atmospheric conditions. We find that some of the main challenges in wind turbine wake modeling are related to simulating the inflow. In the neutral benchmark, model performance tracked as expected with model fidelity, with large-eddy simulations performing the best. In the more challenging stable case, steady-state Reynolds-averaged Navier–Stokes simulations were found to outperform other model alternatives because they provide the ability to more easily prescribe noncanonical inflows and their low cost allows for simulations to be repeated as needed. Dynamic measurements were only available for the unstable benchmark at a single downstream distance. These dynamic analyses revealed that differences in the performance of time-stepping models come largely from differences in wake meandering. This highlights the need for more validation exercises that take into account wake dynamics and are able to identify where these differences come from: mesh setup, inflow, turbulence models, or wake-meandering parameterizations. In addition to model validation findings, we summarize lessons learned and provide recommendations for future benchmark exercises.

17 WIND ENERGY↗

DOE OSTI2022

Transition from sub-Rayleigh anticrack to supershear crack propagation in snow avalanches

Snow slab avalanches, characterized by a distinct, broad fracture line, are released following anticrack propagation in highly porous weak snow layers buried below cohesive slabs. The anticrack mechanism is driven by the volumetric collapse of the weak layer, which leads to the closure of crack faces and to the onset of frictional contact. Here, on the basis of snow fracture experiments, full-scale avalanche measurements and numerical simulations, we report the existence of a transition from sub-Rayleigh anticrack to supershear crack propagation. This transition follows the Burridge–Andrews mechanism, in which a supershear daughter crack nucleates ahead of the main fracture front and eventually propagates faster than the shear wave speed. Furthermore, we show that the supershear propagation regime can exist even if the shear-to-normal stress ratio is lower than the static friction coefficient as a result of the loss of frictional resistance during collapse. This finding shows that snow slab avalanches have fundamental similarities with strike-slip earthquakes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Planck 2018 results XII. Galactic astrophysics using polarized dust emission

Observations of the submillimetre emission from Galactic dust, in both total intensity I and polarization, have received tremendous interest thanks to the Planck full-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release of Planck data. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fraction p, polarization angle ψ, and dispersion function of polarization angles S . The joint distribution (one-point statistics) of p and NH confirms that the mean and maximum polarization fractions decrease with increasing NH. The uncertainty on the maximum observed polarization fraction, pmax = 22.0 -1.4 +3.5% at 353 GHz and 80' resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour between p and S found earlier is seen to be present at high latitudes. This follows the S ∝ p -1 relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties of p, ψ, and S for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map S × p, looking for residual trends. While the polarization fraction p decreases by a factor of 3-4 between NH = 10 20 cm -2 and NH = 2 x 10 22 cm -2 , out of the Galactic plane, this product S × p only decreases by about 25%. Because S is independent of the grain alignment efficiency, this demonstrates that the systematic decrease in p with N H is determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of S × p with the dust temperature T d , whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fraction p and low polarization angle dispersion S tend, on the contrary, to have colder dust than lines of sight with low p and high S . We also compare the Planck thermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of the Planck beam. The two polarization emission-to-extinction ratios, R P/p and R S/V , which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction, p V /E(B - V), of 13% at high Galactic latitude, compatible with the polarization fraction p ≈ 20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Spin orbit torque switching of synthetic Co/Ir/Co trilayers with perpendicular anisotropy and tunable interlayer coupling

Spin orbit torque (SOT) has attracted much attention as an energy efficient electrical method to switch the magnetization in single magnetic layer with both in-plane and perpendicular anisotropy. Here, we report SOT switching of synthetic antiferromagnetic and ferrimagnetic Co/Ir/Co trilayers, where the net magnetization and the interlayer coupling strength are highly controllable. Here, a weak external field dependence of the SOT switching was observed due to the robust domain wall structure in the trilayer. The switching current density was observed to scale inversely with the net magnetization. For trilayers with the same net magnetization, switching current is independent of the layer order because switching in the thicker Co layer dominates. In both ferromagnetic and antiferromagetic regimes, the switching current reflects the interlayer coupling strength.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Characterization and controllability of radiated power via extrinsic impurity seeding in strongly negative triangularity plasmas in DIII-D

Experiments with extrinsic impurity seeding in strongly negative triangularity shapes in DIII-D achieved radiated power fractions (relative to input power) of up to ≈85% total radiation and ≈55% core radiation in steady-operating conditions. The relationship between core and total radiation was sensitive to impurity species and input power. Attempts to reach higher radiation levels via higher impurity flows resulted in radiative collapse disruptions. Nitrogen, neon, argon, and krypton were tested. Injection was by gas puffing, usually controlled by feeding back real-time estimates for core or total radiating fractions ($P_\textrm{rad} / P_\textrm{input}$). Argon and krypton were controlled by feeding back total radiating fraction, whereas neon was controlled by feeding back the core radiating fraction, due to the lower efficiency of neon as a divertor radiator. Nitrogen flows were pre-programmed. Reasonable total $P_\textrm{rad}$ control target following was achieved with argon or krypton. Control was more challenging with the neon/core radiation configuration, which was more prone to slow response and overshooting of the control target. Poor particle removal contributed to the control challenge: neon particle inventory within the last closed flux surface was roughly constant for up to 1 s (the longest duration tested) after neon injection was halted. However, separate experiments with laser-blow off of non-recycling impurities measured a short impurity confinement time, on the scale of the energy confinement time of ∼100 ms. Modeling with the Aurora impurity transport simulation matched experimental neon density profiles with full recycling (R = 1.0) and weak pumping but predicted rapid decreases in neon inventory if pumping were increased or recycling decreased. This indicates that changes outside the confined plasma (adding a well-placed pump) would improve controllability for all highly recycling species.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Disentanglement of density and rotation dependences of the field penetration threshold on the J-TEXT tokamak

The field penetration threshold of magnetic perturbations has been observed to vary non-monotonically with an increase of density in ohmic plasmas on the J-TEXT tokamak. This observation appears contradicting the previous empirical density scaling law. Disentanglement of plasma density and rotation dependences of the field penetration threshold has been carried out. It shows that the field penetration threshold depends only weakly on the density but linearly on the plasma rotation. This result is not only important for the prediction of error field tolerance in fusion devices, but also opens a question on the role of density in the forced magnetic reconnection process in magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Responses of summer mesoscale convective systems to irrigation over the North China Plain based on convection-permitting model simulations

Extensive irrigation activities in the North China Plain (NCP) significantly influence regional weather and climate. However, previous studies focusing on the NCP were primarily based on coarse-resolution models, which are unable to explicitly resolve convection systems, causing large uncertainty in precipitation simulations. In this study, a convection-permitting model coupled with a dynamic irrigation scheme is utilized to investigate the impacts of irrigation on summertime mesoscale convective systems (MCSs) over the NCP. Sensitivity experiments with irrigation off and on are conducted for 5 summers and an MCS identification and tracking algorithm is applied to both satellite observations and model simulations. We find that incorporating irrigation in the model increases MCS precipitation, which agrees more with observations. The probability distributions of MCS lifetime, area, propagation speed, and intensity are all better simulated with irrigation. Irrigation increases the occurrence frequency of MCSs throughout the entire day. The nighttime increase is partly because of more frequent local initiation of MCS developed from isolated deep convection, while the daytime increase is mainly attributed to the changes in MCSs initiating elsewhere and then propagating to the NCP. On average, irrigation induces additional moisture that is more thermodynamically favorable for precipitation, but this effect is partially offset by the weakened ascending air motion primarily caused by irrigation surface cooling. Compared to weak MCS precipitation events, strong MCS precipitation events experience greater enhancement in precipitation intensity when including irrigation because the offset effect from the change in large-scale ascending air motion is insignificant. In addition, irrigation makes the variation of MCS precipitation intensity more correlated with the variation in ascending motion but less correlated with that in atmospheric moisture. Our results suggest the pronounced impacts of irrigation on MCSs over the NCP which should be included in numerical models to improve regional precipitation simulation and prediction.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Quasiparticle Energy Relaxation in a Gas of One-Dimensional Fermions with Coulomb Interaction

We consider a system of charged one-dimensional spin-1/2 fermions at low temperature. We study how the energy of a highly excited quasiparticle (or hole) relaxes toward the chemical potential in the regime of weak interactions. The dominant relaxation processes involve collisions with two other fermions. We find a dramatic enhancement of the relaxation rate at low energies, with the rate scaling as the inverse sixth power of the excitation energy. Furthermore, this behavior is caused by the long-range nature of the Coulomb interaction.

1-dimensional systems↗

DOE OSTI2025

Magnitude and timescale of liquid water path adjustments to cloud droplet number concentration perturbations for nocturnal non-precipitating marine stratocumulus

Cloud liquid water path (L) adjusts to perturbations in cloud droplet number concentration (N) over time. We explore the magnitude and timescale of this adjustment in nocturnal non-precipitating marine stratocumuli using large eddy simulations of baseline conditions and aerosol seeding experiments for 22 meteorological conditions. The results confirm that the L adjustment (δL) slope (k) is more negative for simulation pairs with relatively low N and less negative for high N. Overall, k is unlikely to be lower than −0.4 within 24 h of seeding start, meaning the L adjustment is unlikely to fully offset the brightening due to the Twomey effect. After seeding, the δL becomes increasingly negative, which can be characterized by an exponential convergence. This evolution is governed by a short timescale around 5 h and lasts for around 8–12 h. It is driven by the feedback between entrainment, L, and boundary layer (BL) turbulence. Other processes, including radiation, surface fluxes, and subsidence, respond to the seeding weakly. This short timescale is insensitive to the amount of seeding, making the evolution of δL and some other deviations similar for different seeding amounts after appropriate scaling. The timescale of k evolution is closely related to the δL timescale and hence also short, while it could also be affected by the δN evolution. The results are most relevant to conditions where seeding is applied to a large area of marine stratocumulus in well-mixed and overcast BL where shear is not a primary source of turbulence.

54 ENVIRONMENTAL SCIENCES↗

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↗