Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME DEPENDENCE”

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.

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

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Solving MaxCut with quantum imaginary time evolution

We introduce a method to solve the MaxCut problem efficiently based on quantum imaginary time evolution (QITE). We employ a linear Ansatz for unitary updates and an initial state involving no entanglement, as well as an imaginary-time-dependent Hamiltonian interpolating between a given graph and a subgraph with two edges excised. We apply the method to thousands of randomly selected graphs with up to fifty vertices. We show that our algorithm exhibits a 93% and above performance converging to the maximum solution of the MaxCut problem for all considered graphs. Our results compare favorably with the performance of classical algorithms, such as the greedy and Goemans–Williamson algorithms. We also discuss the overlap of the final state of the QITE algorithm with the ground state as a performance metric, which is a quantum feature not shared by other classical algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A Hybrid AI/ML and Computational Mechanics Based Approach for Time-Series State and Fatigue Life Estimation of Nuclear Reactor Components

Environmental fatigue modeling is a complex problem due to multiple failure modes and their intermixing. The failure modes are function of various underlying causes in addition to the corrosive effect of reactor coolant environment. Some of the major causes are time-dependence of material associated with cyclic loading, load sequence effect associated with random/variable amplitude loading, effect of strain amplitude and rates, effect of varying temperature (along both temporal and spatial directions) and the effect of mean strain and stress. The nonlinear intermixing of failure modes associated with above mentioned causing parameters makes the environmental fatigue modeling is a challenging task. Because of this challenge, fatigue is traditionally being modeled based on experimental data. However, test based empirical approach often requires hundreds of fatigue tests to model the above-mentioned intermixing failure causes even for a single material system. The problem is further exaggerated for reactor component made from multi-material systems such as made from both carbon and stainless-steel base metals and their similar and dissimilar metal welds. With the difficulty of conducting hundreds of fatigue tests to capture the above-mentioned intermixing failure causes, fatigue modeling approaches often depends on empirical models based on limited available test data such as available through ASME code and NUREG 6909. However, these limited test-data-based models may not be enough to accurately predict the life of reactor components. Accurate prediction of life of reactor component would become a necessity, particularly when the license of the reactors to be extended for long-term-operation (LTO) that is for well beyond its original design life of 40 years. The requirement of extending the license of reactor under LTO requires hundreds of fatigue tests to be conducted to understand the mechanism associated with the above-mentioned interdependent failure causes. However, conducting large number of fatigue tests is not a feasibility due to the cost involved. To address this issues Argonne National Laboratory (ANL) with the sponsorship of DOE Light Water Reactor Sustainability (LWRS) program trying to develop a hybrid predictive modeling approach. This is based on limited experiment-data, Artificial-intelligence (AI) – Machine-Learning (ML) - Deep-Learning (DL) based techniques and Multiphysics-computational-mechanics based modeling tools. The hybrid approach not-only can improve the accuracy of the existing stress analysis and fatigue modeling approach but also can reduce the over-dependency on test-based approach. Towards this goal following are some of the major contributions based on ANL’s FY-20 environmental fatigue modeling activities: 1) A cyclic plasticity material model database for 82/182 dissimilar metal weld, which can be readily shared with US nuclear industry and regulatory agency on request. 2) A well validated analytical modeling methodology to perform cycle-by-cycle stress prediction under both constant amplitude fatigue loading and variable amplitude fatigue loading (with load-sequence effect). 3) An AI/ML/DL based methodology to predict unmeasurable cyclic strain based on other available sensor signals. This type of approach can be used for estimating strain in real reactor components from other sensor readings. 4) An AI/ML based approach to improve the US capability on environmental fatigue testing. This is by improving ANL’s existing environmental fatigue testing capacity to conduct ASME required strain-controlled tests (by controlling strain amplitudes and its rate), while not measuring the strain (due to the difficulty of placing an extensometer in a narrow autoclave in a PWR-water-test system). 5) A simulation and experiment based probabilistic modeling methodology for time-series fatigue state and life estimation of reactor metal such as dissimilar metal weld.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Rehydration of metastudtite in the alteration kinetics of α– and β–U 3 O 8 in dilute aqueous solutions of hydrogen peroxide

The formation of alteration phases on uranium ore concentrates and used nuclear fuels under oxidizing conditions is key to understanding the potential mobility of radionuclides in the environment and designing optimal storage conditions of materials. However, the time-dependent distribution of alteration phases on α– and β–U 3 O 8 under oxidizing conditions has yet to be explored. Here, in this study, crystalline powders of α– and β–U 3 O 8 were submerged in aqueous solutions of hydrogen peroxide (1.6 × 10 −1 to 5.4 × 10 −5 M) with aliquots of solution and solid removed for analysis at 1, 8, 15, 22, 29, 36, 46, 58, 71, and 83 days. Within one day there is significant alteration of U 3 O 8 to the uranyl peroxide metastudtite, [(UO 2 )(O 2 )(H 2 O) 2 ], that is replaced by studtite, [(UO 2 )(O 2 )(H 2 O) 2 ]·2H 2 O, within a week regardless of the polymorph of U 3 O 8 or the initial concentration of H 2 O 2 in solution, as determined by partial least squares regression (PLSR) of Raman spectra collected from the solids. A dissolution/reprecipitation mechanism is proposed for both the alteration of U 3 O 8 to metastudtite and the subsequent alteration of both U 3 O 8 and metastudtite to studtite. The two polymorphs of U 3 O 8 exhibit similar rates and extents of alteration over time. The rehydration of metastudtite to studtite has not been previously reported and highlights the need for future work to determine the mechanism by which metastudtite is converted to studtite and what other conditions facilitate this rehydration.

Alteration of U3O8↗
DOE OSTI2022

Photoinitiator Selection and Concentration in Photopolymer Formulations towards Large-Format Additive Manufacturing

Photopolymers are an attractive option for large-format additive manufacturing (LFAM), because they can be formulated from structural thermosets and cure rapidly in ambient conditions under low-energy ultraviolet light-emitting diode (UV LED) lamps. Photopolymer cure is strongly influenced by the depth penetration of UV light, which can be limited in the 2–4 mm layer thicknesses typical of LFAM. Photoinitiator (PI) systems that exhibit photobleaching have proven useful in thick-section cure applications, because they generate a photoinitiation wavefront, but this effect is time-dependent. This study investigates the light transmission and through-thickness cure behavior in (meth)acrylate photopolymer formulations with the photobleaching initiator bis(2,4,6-trimethylbenzoyl)-phenylphosphine oxide (BAPO). Utilizing an optical model developed by Kenning et al., lower concentrations (0.1 wt% to 0.5 wt%) of BAPO were predicted to yield rapid onset of photoinitiation. In situ cure measurements under continuous UV LED irradiation of 380 mW/cm2 showed that a 0.1 wt% concentration of BAPO achieved peak polymerization rate within 2.5 s at a 3-mm depth. With only 1 s of irradiation at 1.7 W/cm2 intensity, the 0.1 wt% BAPO formulation also achieved the highest level of cure of the formulas tested. For an irradiation dose of 5.5 J/cm2 at a duration of 3.7 s, cured polymer specimens achieved a flexural strength of 108 MPa and a flexural modulus of 3.1 GPa. This study demonstrates the utility of optical modeling as a potential screening tool for new photopolymer formulations, primarily in identifying an upper limit to PI concentration for the desired cure depth. The results also show that photobleaching provides only a limited benefit for LFAM applications with short (1.0 s to 3.7 s) UV irradiation times and indicate that excess PI concentration can inhibit light transmission even under extended irradiation times up to 60 s.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Optimizing multigrid reduction-in-time and Parareal coarse-grid operators for linear advection

Parallel-in-time methods, such as multigrid reduction-in-time (MGRIT) and Parareal, provide an attractive option for increasing concurrency when simulating time-dependent partial differential equations (PDEs) in modern high-performance computing environments. While these techniques have been very successful for parabolic equations, it has often been observed that their performance suffers dramatically when applied to advection-dominated problems or purely hyperbolic PDEs using standard rediscretization approaches on coarse grids. In this paper, we apply MGRIT or Parareal to the constant-coefficient linear advection equation, appealing to existing convergence theory to provide insight into the typically nonscalable or even divergent behavior of these solvers for this problem. To overcome these failings, we replace rediscretization on coarse grids with improved coarse-grid operators that are computed by applying optimization techniques to approximately minimize error estimates from the convergence theory. Therefore, one of our main findings is that, in order to obtain fast convergence as for parabolic problems, coarse-grid operators should take into account the behavior of the hyperbolic problem by tracking the characteristic curves. Our approach is tested for schemes of various orders using explicit or implicit Runge–Kutta methods combined with upwind-finite-difference spatial discretizations. In all cases, we obtain scalable convergence in just a handful of iterations, with parallel tests also showing significant speed-ups over sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Implementation of the D1S Methodology for Shutdown Dose Rate Calculations in the OpenMC Monte Carlo Particle Transport Code

We present an implementation of the direct one-step (D1S) methodology for shutdown dose rate (SDR) calculations in the OpenMC Monte Carlo particle transport code. In addition to being the first fully open-source D1S implementation, it is also the first to require no ad hoc source code or nuclear data library modifications. The code can seamlessly switch between production of prompt and decay photons based on a user input parameter, and the decay data needed for decay photon generation are made available through a depletion chain file, which is already used for OpenMC’s built-in depletion/activation solver. A set of Python functions significantly eases the burden of computing and applying time correction factors needed to properly account for the time dependence of radionuclide activity. To assess the accuracy of the D1S implementation, SDR calculations have been carried out for three problems: a prism of iron irradiated by 14-MeV neutrons, the ITER port plug computational benchmark, and the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For each of these problems, comparisons were made to calculations using the rigorous two-step (R2S) method. The results on the iron prism problem illustrate how the D1S method achieves superior spatial resolution compared to the R2S method without the need for spatial discretization of the activation regions. The D1S and R2S results for the ITER port plug benchmark agree well with previously reported results in the literature. While the D1S results are 10% to 15% lower than the R2S results, this may be due to stochastic uncertainty and/or spatial discretization in the R2S calculations. On the FNG dose rate benchmark problem, the D1S method produces dose rate estimates that are within 4% of the dose rates predicted using a cell-based R2S workflow. The D1S estimates of the SDR are also in reasonable agreement with the experimental measurements and show the same basic trends that have been observed in previous works. A qualitative analysis of the execution time and uncertainty for the R2S and D1S workflows suggests that the D1S method would attain a higher figure of merit.

D1S method↗
DOE OSTI2025

Accessing Transient Isomers in the Photoreaction of Metastable‐State Photoacid

The photoreaction of a metastable‐state photoacid (mPAH) generally involves multiple isomers with various connected pathways of photoinduced structural changes during a single reaction cycle. However, only a limited number of isomers have been identified experimentally so far owing to the inherent complexity in combination with the presence of various competing electronic and vibrational processes, as well as the constantly varying interactions between mPAH isomers and solvent molecules. Here, in this work, an optical spectroscopic study on a benzimidazole‐based mPAH, a novel photoacid using benzimidazole as the structural moiety with the active proton, is reported. Through measurements of linear absorption and steady‐state fluorescence in neat solvents and binary mixtures, a pronounced effect of neat water and its binary mixture with glycerol is discovered on the photoreaction of this benzimidazole‐mPAH, manifested by the remarkably distinct spectral responses to irradiation from that observed for an organic solution under identical conditions. Measurements of time‐ and frequency‐resolved fluorescence emission further enable to access transient isomers and the associated spectral characteristics from other competing electronic excited‐state relaxation processes. Spectral deconvolution analysis and time‐dependent density functional theory calculations are applied to separate distinct spectral components and access their potential origin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Unraveling the hysteretic behavior at double cations-double halides perovskite - electrode interfaces

Despite over a decade of research on metal halide perovskites (MHPs) in the context of photovoltaic applications, understanding the nature of electronic and ionic processes associated with current-voltage (I-V) hysteretic behavior has been limited. Here, we explore the hysteretic behavior in (FAPbI 3 ) 0.85 (MAPbBr 3 ) 0.15 perovskite devices with lateral Cr electrodes by applying first order reversal curve (FORC) bias waveform in I-V, Kelvin probe force microscopy (KPFM) measurements, and in-situ chemical imaging by time-resolved time-of-flight secondary ion mass spectrometry (tr-ToF-SIMS). In dark, we reveal pronounced hysteretic behaviors of charge dynamics in the off-field by probing time-dependent current and contact potential difference (CPD). Under illumination, transient and hysteretic behaviors are significantly reduced. The tr-ToF-SIMS results reveal that the hysteretic behaviors are strongly associated with accumulation of Br- ions at the interfaces. In addition, the low mobility MA + ions result in transient behavior and contribute to the hysteretic phenomena. It was shown that Pb 2+ ions can be reduced at the interfaces due to electrochemical reactions with the electrode in the presence of charge injection and photogenerated charges. Furthermore, these hysteretic behaviors associated with charge dynamics, ion migration, and interfacial electrochemical reaction are critical to further improve the performance and stability of MHPs photovoltaics and optoelectronics.

42 ENGINEERING↗
DOE OSTI2021

Quantification of Mono- and Multivalent Counterion-Mediated Bridging in Polyelectrolyte Brushes

Multivalent counterion-induced bridging interactions have been identified as the key mechanism of drastic collapse of the height of polyelectrolyte (PE) brushes. In this article, we employ all-atom molecular dynamics (MD) simulations to quantify the bridging interactions in PE brushes for counterions of different sizes and valences. We identify that unlike the current notion, bridging interactions are not the sole function of the counterion valence. Rather the bridging interactions depend on the fraction of counterions (of a given type) that get physically condensed on the PE backbone as well as the size of the counterion solvation shell. These mechanisms ensure that certain monovalent counterions demonstrate much stronger bridging interactions than that witnessed for certain divalent and trivalent counterions, while certain counterions of identical valences show drastically different bridging. Here, we argue that these counterion-specific bridging interactions eventually enable not only the significant reduction of the PE brush height in presence of certain multivalent screening counterions, but may also give rise to scenarios where the brush height reduction for certain monovalent counterions is larger than certain divalent and trivalent counterions. This latter observation contradicts the experimental findings where the multivalent counterions invariably led to a larger decrease in the height of the PE brushes: we argue that this discrepancy stems from the fact that in our simulations we only consider densely grafted and short (and hence less flexible) PE brushes that hinder the formation of different laterally inhomogeneous structures (like pinned micelles and cylindrical bundles) that would have led to a larger brush height reduction (in experiments, which invariably consider longer and less densely grafted brushes, the formation of such inhomogeneous structures are primarily responsible for larger brush height reduction in presence of multivalent counterions). Finally, we also probe the dynamic properties of the counterions (i.e., their time-dependent displacements) and their bridging interactions (i.e., lifetime of bridging interactions).

42 ENGINEERING↗
DOE OSTI2024

A generalized Selberg zeta function for flat space cosmologies

Flat space cosmologies (FSCs) are time dependent solutions of three-dimensional (3D) gravity with a vanishing cosmological constant. They can be constructed from a discrete quotient of empty 3D flat spacetime and are also called shifted-boost orbifolds. Using this quotient structure, we build a new and generalized Selberg zeta function for FSCs, and show that it is directly related to the scalar 1-loop partition function. We then propose an extension of this formalism applicable to more general quotient manifolds $\mathcal{M}$/ℤ, based on representation theory of fields propagating on this background. Our prescription constitutes a novel and expedient method for calculating regularized 1-loop determinants, without resorting to the heat kernel. We compute quasinormal modes in the FSC using the zeroes of a Selberg zeta function, and match them to known results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Examination of How Well Long-Range-Corrected Density Functionals Satisfy the Ionization Energy Theorem

For this work, we calculated the vertical ionization energies (VIE) of 99 species in two ways to examine the accuracy of several long-range-corrected (LC) hybrid meta functionals in comparison with a gradient approximation (GA), global hybrids, and doubly hybrids. In the category of LC functionals, we examined both those with meta ingredients (i.e., that depend on the kinetic energy density) and those without them. The LC-hybrid meta functionals examined are M11, revM11, M11plus, and ωB97M-V. The reference data used to assess accuracy consist of 95 molecules and 4 atoms in the GW100 set. The two methods studied are the ΔSCF method (involving the difference of neutral and cation self-consistent field (SCF) energies) and the ionization energy theorem (involving the orbital energy of the highest occupied molecular orbital, HOMO). We calculated linear correlation coefficients (r 2 ) and mean absolute deviations (MADs) between each approach and the reference VIE value from the CCSD(T)/def2-TZVPP level of theory. We compared the new LC-hybrid meta calculations to calculations with the 10 functionals in a previous VIE study by Brémond et al. and to the calculations with LC-BLYP (LC-Becke, Lee–Yang–Parr), CAM-B3LYP (Coulomb-attenuating-method Becke-3-parameter Lee–Yang–Parr), LC-ωHPBE, and ωB97X-D. The results show that Minnesota LC-hybrid meta functionals have the smallest mean absolute deviation of ionization energy theorem VIEs with the reference data; the LC-ωHPBE functional also does quite well in this test. This is very encouraging and indicates that LC-hybrid meta functionals would be the best starting points for the tuning strategy that has been shown to be a very good procedure for improving time-dependent density functional calculations, and it also helps explain the good success of LC-hybrid meta functionals for molecular excitation energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Long-term permeability evolution of shale seal rocks with argon and scCO 2

In this work, we conducted a suite of experiments to evaluate the long-term permeability evolution of shale rocks under constant effective stress, before, during, and after interaction with supercritical carbon dioxide (scCO 2 ). To do so, we measured the time-dependent evolution of argon permeabilities for the pre- and post-reacted samples (i.e., before and after long-term interaction with scCO 2 ). In addition, we obtained permeability evolution during long-term interaction with scCO 2 . The samples showed either relatively constant permeabilities or a moderate decrease during pre reaction long-term argon tests. The permeability evolution during long-term CO 2 tests showed continuous increase, continuous decrease, or cycles of increase/decrease in permeability. The long term response of the samples to CO 2 included phenomena such as 20 (i) salt precipitation, (ii) swelling-induced cracks, and (iii) carbonate dissolution. While it is obvious that salt precipitation and swelling-induced cracks decrease and increase the permeability, respectively, the sample response to carbonate dissolution proved to be more complex and may increase or decrease the permeability. The permeability evolution during post-reaction long-term argon injection is also affected by the contribution of each of these three phenomena, during long-term interaction with CO 2 . We observe increase, decrease, and constant permeability evolution during post-reaction argon tests. Our experiments reveal that the initial permeability of the samples plays a significant role on the long-term permeability response of shales in the presence argon and scCO 2 fluids. This study shows, when shales are hydraulically fractured with CO 2 , their initial permeability has a more significant role than their permeability evolution over time.

58 GEOSCIENCES↗
DOE OSTI2020

Visualizing Energy Transfer at Buried Interfaces in Layered Materials Using Picosecond X-Rays

Understanding the fundamentals of nanoscale heat propagation is crucial for next-generation electronics. For instance, weak van der Waals bonds of layered materials are known to limit their thermal boundary conductance (TBC), presenting a heat dissipation bottleneck. In this paper we present a new non-destructive method to probe heat transport in nanoscale crystalline materials using time-resolved x-ray measurements of photo-induced thermal strain. This technique directly monitors time-dependent temperature changes in the crystal and the subsequent relaxation across buried interfaces by measuring changes in the c-axis lattice spacing after optical excitation. We investigate films of five different layered transition metal dichalcogenides MoX 2 [X= S, Se, Te] and WX 2 [X= S, Se] as well as graphite and a W-doped alloy of MoTe 2 , on c-plane sapphire substrates, finding TBC values in the range 10-30 MWm -2 K -1 at room temperature. In conjunction with molecular dynamics simulations, we show that the high thermal resistances are a consequence of weak interfacial van der Waals bonding and low phonon irradiance. This work paves the way for an improved understanding of thermal bottlenecks in emerging three-dimensional heterogeneously integrated technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Runaway electron deconfinement in SPARC and DIII-D by a passive 3D coil

The operation of a 3D coil—passively driven by the current quench (CQ) loop voltage—for the deconfinement of runaway electrons (REs) is modeled for disruption scenarios in the SPARC and DIII-D tokamaks. Nonlinear magnetohydrodynamic (MHD) modeling is carried out with the NIMROD code including time-dependent magnetic field boundary conditions to simulate the effect of the coil. Further modeling in some cases uses the ASCOT5 code to calculate advection and diffusion coefficients for REs based on the NIMROD-calculated fields, and the DREAM code to compute the runaway evolution in the presence of these transport coefficients. Compared with similar modeling in Tinguely et al (2021 Nucl. Fusion 61 124003), considerably more conservative assumptions are made with the ASCOT5 results, zeroing low levels of transport, particularly in regions in which closed flux surfaces have reformed. Of three coil geometries considered in SPARC, only the n = 1 coil is found to have sufficient resonant components to suppress the runaway current growth. Without the new conservative transport assumptions, full suppression of the RE current is maintained when the thermal quench MHD is included in the simulation or when the RE current is limited to 250kA, but when transport in closed flux regions is fully suppressed, these scenarios allow RE beams on the order of 1–2 MA to appear. Additional modeling is performed to consider the effects of the close ideal wall. In DIII-D, the CQ is modeled for both limited and diverted equilibrium shapes. In the limited shape, the onset of stochasticity is found to be insensitive to the coil current amplitude and governed largely by the evolution of the safety-factor profile. In both devices, prediction of the q-profile evolution is seen to be critical to predicting the later time effects of the coil.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

On the theory of nonhomogeneous nonequilibrium superconductivity in 2D systems with massless fermions

Here we analyze static and nonequilibrium superconducting properties of a 2D relativistic-like model system with local electron-electron interaction, Rashba spin-orbit interaction αR in presence of time-dependent in-plane magnetic field H(t). It is shown that similar to the 2D case with ordinary massive quasiparticle dispersion ε(k)~|k| 2 at large fields, such a system demonstrates a nonhomogeneous superconducting stripe phase with the order parameter Δ(r)=Δ(0)cos(2[μ B B×r]n/$\hbarυ_F$) (B is the magnetic induction, υF is the Fermi velocity, n is the normal to the plane, μ B is the Bohr magneton, and αR$\ll$υF) where the stripes are oriented along the B direction. In the considered system, the inter-stripe period L and the magnitude of the magnetic field B are related by a universal relation BL=$\hbarυ_F$/μ B ≃0.714∙10 -4 Tm. Contrary to the case of massive quasiparticles, where the condition α R ~υ F can be, in principle, satisfied by increasing α R or by charge doping (Fermi velocity decreasing), in a relativistic-like system, where υF is doping-independent and one-two orders of magnitude larger than typical Fermi velocity in the “standard” 2D systems, the stripe phase can be the ground state at a rather low doping level. We also analyzed the nonequilibrium properties of the system with a focus on the melting of the stripe order (when the magnetic field is quenched to a lower value) and stripe dynamics (when the field is rotated by 90° degrees) and found several notable results. In particular, it was shown that the stripe domains melt according to law R~1$\sqrt{t}$ at initial times, while at longer times they shrink exponentially. In the case of the flipped magnetic field, the stripe orientation gradually turns from x- to y-direction, and the intermediate “crossed-stripe” phase takes place during times of order of picoseconds. Such a crossed phase is built of periodic superconducting bubbles that potentially may have applications in modern ultrafast superconducting technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Complex and Unusual Excited-State Relaxation Dynamics of 9,9′-Bifluorenylidene Revealed by Comprehensive Time-Resolved Spectroscopy and MRSF-TDDFT Calculations

Design of novel photochemical molecular motors often requires molecular building blocks that exhibit rather unusual photoactivity, for which conventional analyses of spectroscopic data can lead to conflicting interpretations. We here systematically investigated the excited-state relaxation dynamics of one such molecule, 9,9′-bifluorenylidene (BF), through comprehensive and complementary integration of ultrafast transient absorption (TA) and femtosecond stimulated Raman (FSRS) spectroscopies and first-principles mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT). TA and FSRS identified two sequentially formed transients following photoexcitation. The decay kinetics of the two intermediates differ in response to excitation wavelengths and the viscosity/polarity of solvents. MRSF-TDDFT calculations reveal a direct, barrierless internal-conversion pathway from the bright Franck−Condon state to a dark S1 minimum, where the excited-state population is transiently trapped, accounting for the first transient species observed in spectroscopic experiments. Further tracking down along the PES with MRSF-TDDFT mapped out two nonradiative relaxation pathways via conical intersections that connect the dark S 1 state to three configurations in the ground-state manifolds, within which a ring structure with a C8−C8′ bond and the vibrationally excited ground-state BF were identified from spectroscopic and kinetic data. The complexity of relaxation kinetics was attributed to the flexible torsional and twisting motions about the C9−C9′ bridge bond enabled by the diradical character of the S 1 state. These findings clarify unusual photoactive relaxation dynamics stemming from a novel correlation between structural flexibility and shifting electronic characteristics, and they demonstrate the importance of integrating spectroscopic and advanced electronic structure calculation studies for judicious clarification of complex, competing relaxation pathways of excited states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗