SEARCH · Search NASA

Results for “Analytic derivatives”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Numerical and Analytical Modeling of the Effect of Cracks on the Self-Inductance of a COTS YJ-41003-TC Toroid

COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A Celestial route to AdS bulk locality

We prove a precise form of AdS bulk locality by deriving analytical two-sided bounds on bulk Wilson coefficients. Our bounds are on the Wilson coefficients themselves, rather than their ratios, as is typically found in the literature. Inspired by the Celestial amplitudes program in flat space, we perform a Celestial transform of the CFT Mellin amplitude of four identical scalars. Using the crossing symmetric dispersion relation (CSDR), we express the resulting amplitude in terms of crossing symmetric conformal partial waves. The partial waves satisy remarkable positivity properties which along with the unitarity of the CFT prove sufficient to derive the bounds. We then employ our methods in the limit of large AdS radius and recover known bounds on flat space Wilson coefficients and new bounds on their large radius corrections. We check that the planar Mellin amplitude of four stress-tensor multiplets in $\mathcal{N}$ = 4 SYM satisfies our bounds. Finally, using null constraints, we derive a form of low-spin dominance in AdS EFTs. We find that the low-spin dominance is strongest in flat space and weakens as we move away from flat space towards higher AdS curvature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Site-decorated model for unconventional frustrated magnets: Ultranarrow phase crossover and two-dimensional spin reversal transition

Here, the site-decorated Ising model is introduced to advance the understanding and experimental realization of the recently discovered one-dimensional (1D) finite-temperature ultranarrow phase crossover in an external magnetic field, while mitigating the geometric complexities of traditional bond-decorated models. The unconventional frustration and physics are clarified by exactly mapping the 1D site-decorated Ising model in a magnetic field onto a zero-field bond-decorated 𝐽 1 −𝐽 2 Ising model with conventional geometrical frustration. Furthermore, although higher-dimensional Ising models in an external field remain unsolved exactly, an exact solution for a spin-reversal transition—driven by an exotic, hidden half-ice, half-fire state induced by site decoration—is derived. This transition, triggered by a slight variation in temperature or magnetic field—without changing its direction—even in the weak-field limit, offers a promising route toward energy-efficient applications such as data storage and processing. The results suggest that site decoration offers an avenue for materials and device design, particularly in systems such as mixed 𝑑−𝑓 compounds, optical lattices, and neural networks, calling for further studies with site-decorated Heisenberg models. In addition, the site-decorated model offers a rigorous test ground for artificial intelligence (AI) in science, as the analytic derivation of the present results was not only validated but also improved by a general-purpose large language model, inspiring the use of AI as scientific discoverer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Multiple electron beam generation with different energies and comparable currents from a single cathode potential for high power traveling wave tubes (TWTs)

We present a technique for generating multiple electron beams with different energies and comparable currents from a single cathode stalk at a single potential using nested magnetically insulated coaxial diodes (MICDs). The application is for a multi-stream traveling wave tube. Particle-in-cell simulations are performed for an experimental parameter-based geometry where two thin-walled intense relativistic electron beams immersed in a strong uniform magnetic field propagate through a cylindrical vacuum channel. The analytically derived results are obtained by extending Fedosov’s solution for generating a hollow electron beam from an MICD on a cathode stalk in an infinite magnetic field. Two electron beams are generated and accelerated downstream assuming zero initial kinetic energy of the electrons from the cathodes. Results show both electron beam currents ranging from 66 to 2.8 kA with an energy difference ranging from 6% to 27% depending on voltages applied from 100 to 600 kV and the geometry of the two MICDs. An optimal geometry is a crucial factor in achieving the maximum energy difference between the electron beams for comparable currents. The analytical and numerical simulation results show good agreement. We are currently in the process of planning experiments using the electron beam accelerator (SINUS-6 at UNM) to validate the analytical and simulation results.

Islam, K. N. (ORCID:0000000214397710)↗

DOE OSTI2024

Size Effect of Negative Capacitance State and Subthreshold Swing in Van der Waals Ferrielectric Field–Effect Transistors

Analytical calculations corroborated by the finite element modeling show that thin films of Van der Waals ferrielectrics covered by a 2D-semiconductor are promising candidates for the controllable reduction of the dielectric layer capacitance due to the negative capacitance (NC) effect emerging in the thin films. The NC state is conditioned by energy-degenerated poly-domain states of the ferrielectric polarization induced in the films under incomplete screening conditions in the presence of a dielectric layer. Calculations performed for the FET-type heterostructure “ferrielectric CuInP 2 S 6 film—2D-MoS 2 single-layer—SiO 2 dielectric layer” reveal the pronounced size effect of the multilayer capacitance. Derived analytical expressions for the electric polarization and multilayer capacitance allow to predict the thickness range of the dielectric layer and ferrielectric film for which the NC effect is the most pronounced in various Van der Waals ferrielectrics, and the corresponding subthreshold swing becomes much less than the Boltzmann's limit. Obtained results can be useful for the size and temperature control of the NC effect in the steep-slope ferrielectric FETs.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Nonadiabatic Excited-State Molecular Dynamics for Open-Shell Systems

Nonadiabatic Molecular Dynamics (NAMD) of excited states has been widely used in the simulation of photoinduced phenomena. However, the inability to treat bond breaking and forming processes with single-reference electronic structure methods limits their application in photochemistry for extended molecular systems. In this work, the extension of excited-state NAMD for open-shell systems is developed and implemented in the NEXMD software. We present the spin-unrestricted CIS and TDSCF formalism for the ground and excited states, analytical derivatives, and nonadiabatic derivative couplings for the respective potential energy surfaces. Next, this methodology is employed to study the photochemical reaction of three model molecules. The results demonstrate the advantage of the open-shell approach in modeling photochemical reactions, especially involving bond breaking processes. We find that the open-shell method lowers the reaction barrier at the bond-breaking limits resulting in larger calculated photochemical quantum yields compared to the respective closed-shell results. We also address problems related to spin contamination in the open-shell method, especially when molecular geometries are far from equilibrium.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Equation of state of neutron-rich matter in d -dimensions

Nuclear systems under constraints, with high degrees of symmetries and/or collectivities may be considered as moving effectively in spaces with reduced spatial dimensions. We first derive analytical expressions for the nucleon specific energy E 0 (ρ), pressure P 0 (ρ), incompressibility coefficient K 0 (ρ) and skewness coefficient J 0 (ρ) of symmetric nucleonic matter (SNM), the quadratic symmetry energy E sym (ρ), its slope parameter L(ρ) and curvature coefficient K sym (ρ) as well as the fourth-order symmetry energy E sym,4 (ρ) of neutron-rich matter in general d spatial dimensions (abbreviated as "dD") in terms of the isoscalar and isovector parts of the isospin-dependent single-nucleon potential according to the generalized Hugenholtz-Van Hove (HVH) theorem. The equation of state (EOS) of nuclear matter in dD can be linked to that in the conventional 3-dimensional (3D) space by the ε-expansion which is a perturbative approach successfully used previously in treating second-order phase transitions and related critical phenomena and more recently in studying the EOS of cold atoms. The ε-expansion of nuclear EOS in dD based on a reference dimension d f = d – ε is shown to be effective with –1 ≲ ε ≲ 1 starting from 1 ≲ d f ≲ 3 in comparison with the exact expressions derived using the HVH theorem. Moreover, the EOS of SNM (with/without considering its potential part) is found to be reduced (enhanced) in lower (higher) dimensions, indicating in particular that the many-nucleon system tends to be deeper bounded but saturate at higher densities in spaces with lower dimensions. Here, the links between the EOSs in 3D and dD spaces from the ε-expansion provide new perspectives to the EOS of neutron-rich matter.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

A Quantitative Method for In Situ Pump-Beam Metrology in 4D Ultrafast Electron Microscopy

In this work, we report a method for measuring spot size and focusing conditions of the femtosecond (fs) excitation laser in situ at the specimen location in 4D ultrafast electron microscopy (UEM). The method makes use of threshold laser ablation behaviors of thin amorphous carbon membranes. For Gaussian beam profiles and for ablation threshold fluence values, we analytically derive expressions describing the relationship between ablated hole size and the actual laser spot size. Using these expressions, we developed experimental procedures for characterizing the shape and spot size of the pump beam at the specimen. We demonstrate the viability of the approach for incident excitation wavelengths of 343 nm and 515 nm, thus illustrating the applicability of the method to a range of optical wavelengths without modification. Further, we show that by measuring ablated hole size as a function of focusing condition, a full metrological characterization of the Gaussian beam propagation properties can be performed. Finally, we find good agreement for spot sizes determined with this method and with those determined by extrapolation from measurements taken outside the microscope column. Overall, this method is a simple, cost-effective means for accurate and precise determination of key pump-beam parameters in situ at the specimen location in UEM experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Maximum entropy distributions of dark matter in ΛCDM cosmology

Context. Small-scale challenges to ΛCDM cosmology require a deeper understanding of dark matter physics. Aims. This paper aims to develop the maximum entropy distributions for dark matter particle velocity (denoted by X ), speed (denoted by Z ), and energy (denoted by E ) that are especially relevant on small scales where system approaches full virialization. Methods. For systems involving long-range interactions, a spectrum of halos of different sizes is required to form to maximize system entropy. While the velocity in halos can be Gaussian, the velocity distribution throughout the entire system, involving all halos of different sizes, is non-Gaussian. With the virial theorem for mechanical equilibrium, we applied the maximum entropy principle to the statistical equilibrium of entire system, such that the maximum entropy distribution of velocity (the X distribution) could be analytically derived. The halo mass function was not required in this formulation, but it did indeed result from the maximum entropy. Results. The predicted X distribution involves a shape parameter α and a velocity scale, v 0 . The shape parameter α reflects the nature of force ( α → 0 for long-range force or α → ∞ for short-range force). Therefore, the distribution approaches Laplacian with α → 0 and Gaussian with α → ∞. For an intermediate value of α , the distribution naturally exhibits a Gaussian core for v ≪ v 0 and exponential wings for v ≫ v 0 , as confirmed by N -body simulations. From this distribution, the mean particle energy of all dark matter particles with a given speed, v , follows a parabolic scaling for low speeds (∝ v 2 for v ≪ v 0 in halo core region, i.e., “Newtonian”) and a linear scaling for high speeds (∝ v for v ≫ v 0 in halo outskirt, i.e., exhibiting “non-Newtonian” behavior due to long-range gravity). We compared our results against N -body simulations and found a good agreement.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

A duality-based coupling of Cosserat crystal plasticity and phase field theories for modeling grain refinement

High-rate deformation processes of metals entail intense grain refinement and special attention needs to be paid to capture the evolution of microstructure. In this article, a new formulation for coupling Cosserat crystal plasticity and phase field is developed. A common approach is to penalize kinematic incompatibility between lattice orientation and displacement-based elastic rotation. However, this can lead to significant solution sensitivity to the penalty parameter, resulting in low accuracy and convergence rates. To address these issues, a duality-based formulation is developed which directly imposes the rotational kinematic compatibility. A weak inf-sup-based skew-symmetric stress projection is introduced to suppress instabilities present in the dual formulation. An additional least squares stabilization is introduced to suppress the spurious lattice rotation with a suitable parameter range derived analytically and validated numerically. The required high-order continuity is attained by the reproducing kernel approximation. It is observed that equal order displacement-rotation-phase field approximations are stable, which allows efficient employment of the same set of shape functions for all independent variables. The proposed formulation is shown to yield superior accuracy and convergence with marginal parameter sensitivity compared to the penalty-based approach and successfully captures the dominant rotational recrystallization mechanism including block dislocation structures and grain boundary migration.

42 ENGINEERING↗

DOE OSTI2021

A Poisson–Nernst–Planck Model of Ion Transport and Interface Segregation in Metal–Insulator–Semiconductor Structures and Solar Cells

A numerical model that describes the transport of mobile ionic species in metal–insulator–semiconductor (MIS) and photovoltaic (PV) devices subject to temperature and voltage stress is presented. The finite element method (FEM) is used to solve the Nernst–Planck equation while imposing Poisson's equation self-consistently as a restriction for the electrostatic potential. This allows the contribution of the ionic species to the potential to be taken into account. Using a variational formulation eases the implementation of diverse boundary conditions, including the incorporation of segregation kinetics at the device interfaces. Segregation across the dielectric–semiconductor interface is relevant to modeling the electronic device degradation in systems where contamination reaches the semiconductor. The model in closed systems with no-flux boundary conditions is validated first. In the limiting case of low contamination levels with respect to the gate bias, the FEM solution matches analytically derived approximations. Then, the implementation is broadened to include an open boundary at the dielectric–semiconductor interface to account for leakage of ions. The predicted time dependence of the flatband voltage in Na-contaminated MIS test structures agrees well with measurements. This model successfully captures the role of long-range ion transport at concentrations of relevance to electronic and PV device instability and neuromorphic computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Neutrinoless double beta decay versus other probes of heavy sterile neutrinos

We make a comparative study of the neutrinoless double beta decay constraints on heavy sterile neutrinos versus other direct and indirect constraints from both lepton number conserving and violating processes, as a sensitive probe of the extent of lepton number violation and possible interference effects in the sterile sector. We introduce a phenomenological parametrisation of the simplified one-generation seesaw model with one active and two sterile neutrino states in terms of experimentally measurable quantities, such as active-sterile neutrino mixing angles, CP phases, masses and mass splittings. This simple parametrisation enables us to analytically derive a spectrum of possible scenarios between the canonical seesaw with purely Majorana heavy neutrinos and inverse seesaw with pseudo-Dirac ones. We then go on to constrain the simplified parameters of this model from various experiments at the energy, intensity and cosmic frontiers. We emphasise that the constraints from lepton number violating processes strongly depend on the mass splitting between the two sterile states and the relative CP phase between them. This is particularly relevant for neutrinoless double beta decay, which is weakened for small mass splitting and opposite CP parities between the sterile states. On the other hand, neutrinoless double beta decay is especially sensitive for Majorana sterile neutrinos with masses around 0.1-10 GeV.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Leptonic sum rules from flavour models with modular symmetries

Sum rules in the lepton sector provide an extremely valuable tool to classify flavour models in terms of relations between neutrino masses and mixing parameters testable in a plethora of experiments. In this manuscript we identify new leptonic sum rules arising in models with modular symmetries with residual symmetries. These models simultaneously present neutrino mass sum rules, involving masses and Majorana phases, and mixing sum rules, connecting the mixing angles and the Dirac CP-violating phase. The simultaneous appearance of both types of sum rules leads to some non-trivial interplay, for instance, the allowed absolute neutrino mass scale exhibits a dependence on the Dirac CP-violating phase. We derive analytical expressions for these novel sum rules and present their allowed parameter ranges as well as their predictions at upcoming neutrino experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

BPS Wilson loop in $\mathcal{N}$ = 2 superconformal SU(N) “orientifold” gauge theory and weak-strong coupling interpolation

We consider the expectation value $ \left\langle \mathcal{W}\right\rangle $ of the circular BPS Wilson loop in $ \mathcal{N} $ = 2 superconformal SU( N ) gauge theory containing a vector multiplet coupled to two hypermultiplets in rank-2 symmetric and antisymmetric representations. This theory admits a regular large N expansion, is planar-equivalent to $ \mathcal{N} $ = 4 SYM theory and is expected to be dual to a certain orbifold/orientifold projection of AdS 5 × S 5 superstring theory. On the string theory side $ \left\langle \mathcal{W}\right\rangle $ is represented by the path integral expanded near the same AdS 2 minimal surface as in the maximally supersymmetric case. Following the string theory argument in [5], we suggest that as in the $ \mathcal{N} $ = 4 SYM case and in the $ \mathcal{N} $ = 2 SU( N ) × SU( N ) superconformal quiver theory discussed in [19], the coefficient of the leading non-planar 1/ N 2 correction in $ \left\langle \mathcal{W}\right\rangle $ should have the universal λ 3/2 scaling at large ’t Hooft coupling. We confirm this prediction by starting with the localization matrix model representation for $ \left\langle \mathcal{W}\right\rangle $ . We complement the analytic derivation of the λ 3/2 scaling by a numerical high-precision resummation and extrapolation of the weak-coupling expansion using conformal mapping improved Padé analysis.

1/N Expansion↗

DOE OSTI2022

Universality aspects of quantum corrections to transverse momentum broadening in QCD media

We study non-linear quantum corrections to transverse momentum broadening (TMB) of a fast parton propagating in dense QCD matter in the leading logarithmic approximation. These non-local corrections yield an anomalous super-diffusive behavior characterized by a heavy tailed distribution which is associated with Lévy random walks. Using a formal analogy with the physics of traveling waves, we show that at late times the transverse momentum distribution tends to a universal scaling regime. We derive analytic solutions in terms of an asymptotic expansion around the scaling limit for both fixed and running coupling. We note that our analytic approach yields a good agreement with the exact numerical solutions down to realistic values of medium length. Finally, we discuss the interplay between system size and energy dependence of the diffusion coefficient $\hat{q}$ and its connection with the gluon distribution function that is manifest at large transverse momentum transfer.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Black hole induced false vacuum decay from first principles

We provide a method to calculate the rate of false vacuum decay induced by a black hole. The method uses complex tunneling solutions and consistently takes into account the structure of different quantum vacua in the black hole metric via boundary conditions. The latter are connected to the asymptotic behavior of the time-ordered Green’s function in the corresponding vacua. We illustrate the technique on a two-dimensional toy model of a scalar field with inverted Liouville potential in an external background of a dilaton black hole. We analytically derive the exponential suppression of tunneling from the Boulware, Hartle-Hawking and Unruh vacua and show that they are parametrically different. The Unruh vacuum decay rate is exponentially smaller than the decay rate of the Hartle-Hawking state, though both rates become unsuppressed at high enough black hole temperature. We interpret the vanishing suppression of the Unruh vacuum decay at high temperature as an artifact of the two-dimensional model and discuss why this result can be modified in the realistic case of black holes in four dimensions.

2D Gravity↗

DOE OSTI2024

Analytically differentiable metrics for phase stability

Here, in this work, a long-established but sparsely documented method of obtaining semi-analytic derivatives of thermodynamic properties with respect to equilibrium conditions is briefly reviewed and rigorously derived. This procedure is then leveraged to construct general forms of derivatives of the residual driving force, a metric for measuring phase stability used in CALPHAD model optimization, with respect to overall system and individual phase compositions. Applied examples – calculating heat capacity in the Al-Fe system, thermodynamic factors in the Nb-V-W system, and residual driving force derivatives in the Ni-Ti system – demonstrate the versatility, accuracy, and extensibility of this method. Using the developed method, residual driving force gradients can be applied directly in CALPHAD model optimizers, as well as in materials design frameworks, to identify regions of phase stability with an efficient, gradient-based approach.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Stress due to electric charge density distribution in a dielectric slab

The spatial distribution of electric field due to an imposed electric charge density profile in an infinite slab of dielectric material is derived analytically by integrating Gauss’s law. Various charge density distributions are considered, including exponential and power-law forms. Here, the Maxwell stress tensor is used to compute a notional static stress in the material due to the charge density and its electric field. Characteristics of the electric field and stress distributions are computed for example cases in polyethylene, showing that field magnitudes exceeding the dielectric strength would be required in order to achieve a stress exceeding the ultimate tensile strength.

42 ENGINEERING↗