Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Self consistent field methods”

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 289 records · Page 16

Modeling of infrared flux spectra from disk-shaped interstellar dust clouds

Although there is growing observational and theoretical evidence for many disk-shaped objects of astrophysical interest, spherical geometry is assumed in most radiative transfer models. Recently we generalized the quasi-diffusion method developed by Leung (1975, 1976) for spherical geometry to solve the problem of scattering, absorption, and reemission by dust grains in a medium of disk geometry. The method is applicable to a variety of astronomical sources whose dynamics are angular-momentum dominated and hence not accurately treated by spherical geometry, e.g., protoplanetary nebulae, circumstellar disks, bipolar-flow molecular clouds, accretion disks and disk galaxies. Using this technique and realistic grain opacities we construct theoretical models to determine self-consistently the dust temperature distribution and infrared emission from disk-shaped, quiescent dark globules heated externally by the intersteller radiation field. The effects of the following parameters on the temperature structure and the emergent spectrum are studied: grain type (graphite and silicate), optical depth, density inhomogeneity, and degree of disk flattening. The disk models are characterized by an aspect ratio and an optical depth at 0.55 microns. For inhomogeneous models, a gaussian density distribution is assumed such that the ratio of central to surface density is 100. To study the effects of source geometry, we also compare results for models with spherical and disk geometry. In this case both models have the same radius and central optical depth, the disk models having a 1:1 aspect ratio. While the dust temperature distributions in the two cases are very similar, the emergent flux for the disk model depends sensitively on the viewing angle. For clouds which are unresolved, one would expect, since the thermal emission is isotropic in the neighborhood of an emitting grain, and since the emission (in the far infrared) is optically thin, that the emergent flux spectrum should be characteristic only of the dust temperature, and independent of viewing angle. However, because of the lack of complete symmetry and the resulting radiation anisotropy, this is found not to be the case for the disk models. Angle is that it implies a large uncertainty in estimating the radiation dust mass.

Spagna, G. F., Jr.↗

Improving initialization and evolution accuracy of cosmological neutrino simulations

Neutrino mass constraints are a primary focus of current and future large-scale structure (LSS) surveys. Non-linear LSS models rely heavily on cosmological simulations — the impact of massive neutrinos should therefore be included in these simulations in a realistic, computationally tractable, and controlled manner. Here, a recent proposal to reduce the related computational cost employs a symmetric neutrino momentum sampling strategy in the initial conditions. We implement a modified version of this strategy into the Hardware/Hybrid Accelerated Cosmology Code (HACC) and perform convergence tests on its internal parameters. We illustrate that this method can impart Script $\mathcal{O}$(1%) numerical artifacts on the total matter field on small scales, similar to previous findings, and present a method to remove these artifacts using Fourier-space filtering of the neutrino density field. Moreover, we show that the converged neutrino power spectrum does not follow linear theory predictions on relatively large scales at early times at the 15% level, prompting a more careful study of systematics in particle-based neutrino simulations. We also present an improved method for backscaling linear transfer functions for initial conditions in massive neutrino cosmologies that is based on achieving the same relative neutrino growth as computed with Boltzmann solvers. Our self-consistent backscaling method yields sub-percent accuracy in the total matter growth function. Comparisons for the non-linear power spectrum with the Mira-Titan emulator at a neutrino mass of m ν = 0.15 eV are in very good agreement with the expected level of errors in the emulator and in the direct N-body simulation.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamical evolution of twisted magnetic flux tubes. I - Equilibrium and linear stability

The three-dimensional dynamical evolution of twisted magnetic flux tubes is studied using a time-dependent magnetohydrodynamic (MHD) model. The flux tubes are intended to model solar coronal loops, and include the stabilizing effect of photospheric line tying. The model permits the complete evolution of flux tubes to be followed self-consistently, including the formation, equilibrium, linear instability, and nonlinear behavior. Starting from an initial uniform background magnetic field, a twisted flux tube is created by the application of slow, localized photospheric vortex flows. The flux tube evolves quasi-statically through sequences of equilibria with increasing twist, until it becomes linearly unstable to an ideal MHD kink mode. In this paper, the equilibrium properties and the linear stability behavior are discussed. The application of the method to the uniform-twist, Gold-Hoyle field confirms the previous stability threshold for kink instability and provides estimates of the resulting growth rate.

Mikic, Zoran↗

Fokker-Planck simulations of fast ion ICRF and electron EC heating in a mirror plasma using CQL3D-m

The CQL3D-m continuum bounce-average Fokker-Planck code is adapted for magnetic mirror plasmas [1] and is now routinely used in no-free-parameter classical integrated modeling of mirror devices [2, 3]. In the present effort, we report on two RF methods of plasma heating in mirror machine. The fast ions (FI) are heated by Fast waves at 2nd-4th harmonic, where FIs originate from neutral beam injection at 45 degrees to the magnetic field. The scenario shows an efficient ion heating near the FI bouncing point. The electrons are heated by X-mode launched from the high magnetic field side towards the resonance. Different from the tokamak applications, CQL3D-m provides an evolving self-consistent ambipolar parallel electric field, which determines the shape of the loss cone and hence an accurate confinement time of both ions and electrons. Also, it includes a description of ion and electron sources and sinks (related to charge exchange and impact ionization) which are updated at every time step. CQL3D-m utilizes a fully nonlinear Coulomb collision operator that is important for the significantly non-Maxwellian ion distributions typically established in mirror plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards a self-consistent analysis of the anisotropic galaxy two- and three-point correlation functions on large scales: application to mock galaxy catalogues

We establish a practical method for the joint analysis of anisotropic galaxy two- and three-point correlation functions (2PCF and 3PCF, respectively) on the basis of the decomposition formalism of the 3PCF using tripolar spherical harmonics. We perform such an analysis with MultiDark-Patchy mock catalogues to demonstrate and understand the benefit of the anisotropic 3PCF. We focus on scales above $80\, h^{-1}\, {\rm Mpc}$, and use information from the shape and the baryon acoustic oscillation (BAO) signals of the 2PCF and 3PCF. We also apply density field reconstruction to increase the signal-to-noise ratio of BAO in the 2PCF measurement, but not in the 3PCF measurement. In particular, we study in detail the constraints on the angular diameter distance and the Hubble parameter. We build a model of the bispectrum or 3PCF that includes the non-linear damping of the BAO signal in redshift space. We carefully account for various uncertainties in our analysis including theoretical models of the 3PCF, window function corrections, biases in estimated parameters from the fiducial values, the number of mock realizations to estimate the covariance matrix, and bin size. The joint analysis of the 2PCF and 3PCF monopole and quadrupole components shows a $30{{\ \rm per\ cent}}$ and $20{{\ \rm per\ cent}}$ improvement in Hubble parameter constraints before and after reconstruction of the 2PCF measurements, respectively, compared to the 2PCF analysis alone. Furthermore, this study clearly shows that the anisotropic 3PCF increases cosmological information from galaxy surveys and encourages further development of the modelling of the 3PCF on smaller scales than we consider.

79 ASTRONOMY AND ASTROPHYSICS↗

Control simulations of many-body quantum systems by a synergism of discrete real-time learning and optimal control theory

We present a self-consistent algorithm for optimal control simulations of many-body quantum systems. The algorithm features a two-step synergism that combines discrete real-time machine learning (DRTL) with Quantum Optimal Control Theory (QOCT) using the time-dependent Schrödinger equation. Specifically, in step (1), DRTL is employed to identify a compact working space (i.e., the important portion of the Hilbert space) for the time evolution of the many-body quantum system in the presence of a control field (i.e., the initial or previously updated field), and in step (2), QOCT utilizes the DRTL-determined working space to find a newly updated control field for a chosen objective. Steps 1 and 2 are iterated until a self-consistent control objective value is reached such that the resulting optimal control field yields the same targeted objective value when the corresponding working space is systematically enlarged. Furthermore, to demonstrate this two-step self-consistent DRTL-QOCT synergistic algorithm, we perform optimal control simulations of strongly interacting 1D as well as 2D Heisenberg spin systems. In both scenarios, only a single spin (at the left end site for 1D and the upper left corner site for 2D) is driven by the time-dependent control fields to create an excitation at the opposite site as the target. It is found that, starting from all spin-down zero excitation states, the synergistic method is able to identify working spaces and convergence of the desired controlled dynamics with just a few iterations of the overall algorithm. In the cases studied, the dimensionality of the working space scales only quasi-linearly with the number of spins.

Artificial neural networks↗

Magnetotransport properties of the topological nodal-line semimetal CaCdSn

Topological nodal-line semimetals support protected band crossings, which form nodal lines or nodal loops between the valence and conduction bands and exhibit novel transport phenomena. Here we address the topological state of the nodal-line semimetal candidate material CaCdSn, and report magnetotransport properties of its single crystals grown by the self-flux method. In this work, our first-principles calculations show that the electronic structure of CaCdSn harbors a single nodal loop around the Γ point in the absence of spin-orbit coupling effects. The nodal crossings in CaCdSn are found to lie above the Fermi level and yield a Fermi surface that consists of both electron and hole pockets. CaCdSn exhibits high mobility (μ ≈ 3.44 × 10 4 cm 2 V –1 s –1 ) and displays a field-induced metal-semiconductor-like crossover with a plateau in resistivity at low temperature. We observe an extremely large and quasilinear nonsaturating transverse as well as longitudinal magnetoresistance (MR) at low temperatures (≈ 7.44 × 10 3 % and ≈ 1.71 × 10 3 %, respectively, at 4 K). We also briefly discuss possible reasons behind such a large quasilinear magnetoresistance and its connection with the nontrivial band structure of CaCdSn.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Perturbative treatment of nonlocal chiral interactions in auxiliary-field diffusion Monte Carlo calculations

Nuclear many-body systems, ranging from nuclei to neutron stars, are some of the most interesting physical phenomena in our universe, and quantum Monte Carlo (QMC) approaches are among the most accurate many-body methods currently available to study them. In recent decades, interactions derived from chiral effective field theory (EFT) have been widely adopted in the study of nuclear many-body systems. One drawback of the QMC approach is the requirement that the nuclear interactions need to be local, whereas chiral EFT interactions usually contain nonlocalities. In this work, we leverage the capability of computing second-order perturbative corrections to the ground-state energy in order to develop a self-consistent approach to including nonlocal operators in QMC calculations. In conclusion, we investigate both the deuteron and the neutron-matter equation of state in order to show the robustness of our technique and pave the way for future QMC calculations at higher orders in the EFT, where nonlocal operators cannot be avoided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Third General Scientific Assembly of the International Association of Geomagnetism and Aeronomy - Special sessions of auroral processes

Methods of timing magnetic substorms, the rapid fluctuations of aurorae, electromagnetic and electrostatic instabilities observed on the field lines of aurorae, the auroral microstructure, and the relationship of currents, electric field and particle precipitation to auroral form are discussed. Attention is given to such topics as D-perturbations as an indicator of substorm onset, the role of the magnetotail in substorms, spectral information derived from imaging data on aurorae, terrestrial kilometric radiation, and the importance of the mirror force in self-consistent models of particle fluxes, currents and potentials on auroral field lines.

Russell, C. T.↗

Carbon Nanotube-Based Structural Health Monitoring Sensors

Carbon nanotube (CNT)-based sensors for structural health monitoring (SHM) can be embedded in structures of all geometries to monitor conditions both inside and at the surface of the structure to continuously sense changes. These CNTs can be manipulated into specific orientations to create small, powerful, and flexible sensors. One of the sensors is a highly flexible sensor for crack growth detection and strain field mapping that features a very dense and highly ordered array of single-walled CNTs. CNT structural health sensors can be mass-produced, are inexpensive, can be packaged in small sizes (0.5 micron(sup 2)), require less power than electronic or piezoelectric transducers, and produce less waste heat per square centimeter than electronic or piezoelectric transducers. Chemically functionalized lithographic patterns are used to deposit and align the CNTs onto metallic electrodes. This method consistently produces aligned CNTs in the defined locations. Using photo- and electron-beam lithography, simple Cr/Au thin-film circuits are patterned onto oxidized silicon substrates. The samples are then re-patterned with a CNT-attracting, self-assembled monolayer of 3-aminopropyltriethoxysilane (APTES) to delineate the desired CNT locations between electrodes. During the deposition of the solution-suspended single- wall CNTs, the application of an electric field to the metallic contacts causes alignment of the CNTs along the field direction. This innovation is a prime candidate for smart skin technologies with applications ranging from military, to aerospace, to private industry.

Wincheski, Russell↗

Flexible all-organic nanocomposite films interlayered with in situ synthesized covalent organic frameworks for electrostatic energy storage

Poly(vinylidene fluoride)-based terpolymers, known for having the largest dielectric constant among the existing dielectric polymers, are attractive materials for electrostatic film capacitors used in lightweight electrification systems. However, the potential of these terpolymers in film capacitor applications remains constrained by their low electrical insulating and mechanical strengths. To address these limitations, we introduced rigid covalent organic framework (COF) nanospheres into the thin films of soft terpolymers via in situ synthesis and a facile layer-by-layer solution casting method, whereby multilayer films consisting of two polymer outer layers and a COF-containing middle layer were readily obtained. The resultant all-organic thin films exhibit simultaneously high dielectric constant, enhanced breakdown strength, superior energy density (~25 J cm –3 ) at efficiencies over 80%, along with greatly improved mechanical self-supporting capability and excellent mechanical flexibility. Here, this work demonstrates the unprecedented use of COF for electrostatic energy storage, uncovering its potential for flexible electronic applications operating under high electric fields.

25 ENERGY STORAGE↗

Dark Energy Survey Year 3 results: redshift calibration of the weak lensing source galaxies

ABSTRACT Determining the distribution of redshifts of galaxies observed by wide-field photometric experiments like the Dark Energy Survey (DES) is an essential component to mapping the matter density field with gravitational lensing. In this work we describe the methods used to assign individual weak lensing source galaxies from the DES Year 3 Weak Lensing Source Catalogue to four tomographic bins and to estimate the redshift distributions in these bins. As the first application of these methods to data, we validate that the assumptions made apply to the DES Y3 weak lensing source galaxies and develop a full treatment of systematic uncertainties. Our method consists of combining information from three independent likelihood functions: self-organizing map p(z) (sompz), a method for constraining redshifts from galaxy photometry; clustering redshifts (WZ), constraints on redshifts from cross-correlations of galaxy density functions; and shear ratios (SRs), which provide constraints on redshifts from the ratios of the galaxy-shear correlation functions at small scales. Finally, we describe how these independent probes are combined to yield an ensemble of redshift distributions encapsulating our full uncertainty. We calibrate redshifts with combined effective uncertainties of σ〈z〉 ∼ 0.01 on the mean redshift in each tomographic bin.

79 ASTRONOMY AND ASTROPHYSICS↗

Theoretical Study of the First Transition Row Oxides and Sulfides

The first transition row oxides and sulfides are studied using several different levels of theory. The calculations show the bonding mechanism in the sulfides and oxides to be very similar. For the oxides, accurate experimental data allow the theoretical methods to be calibrated. The same level of theory is used to study the sulfides where there is far less experimental information. For ScO through MnO and CuO the coupled cluster singles and doubles technique including a perturbational estimate of the unliked triple excitations [CCSD(T)] yields spectroscopic constants ((tau)e, (omega)e, and D0) in good agreement with experiment. The triple excitations are found to be very important in achieving this accuracy. For FeO to NiO, the single determinant self-consistent-field (SCF) approach yields pi orbitals that are localized on the metal or oxygen. This appears to cause problems for the single reference techniques; this is discussed in detail for NiO. The complete-active-space SCF/internally contracted averaged coupled pair functional approach (CASSCF/ICACPF) works well for FeO to NiO. The calculation of accurate dipole moments is found to be very difficult.

Bauschlicher, Charles W., Jt.↗

Mapping Composition Evolution through Synthesis, Purification, and Depolymerization of Random Heteropolymers

Random heteropolymers (RHPs) consisting of three or more comonomers have been routinely used to synthesize functional materials. While increasing the monomer variety diversifies the side-chain chemistry, this substantially expands the sequence space and leads to ensemble-level sequence heterogeneity. Most studies have relied on monomer composition and simulated sequences to design RHPs, but the questions remain unanswered regarding heterogeneities within each RHP ensemble and how closely these simulated sequences reflect the experimental outcomes. Here, we quantitatively mapped out the evolution of monomer compositions in four-monomer-based RHPs throughout a design-synthesis-purification-depolymerization process. By adopting a Jaacks method, we first determined 12 reactivity ratios directly from quaternary methacrylate RAFT copolymerization experiments to account for the influences of competitive monomer addition and the reversible activation/deactivation equilibria. The reliability of in silico analysis was affirmed by a quantitative agreement (<4% difference) between the simulated RHP compositions and the experimental results. Furthermore, we mapped out the conformation distribution within each ensemble in different solvents as a function of monomer chemistry, composition, and segmental characteristics via high-throughput computation based on self-consistent field theory (SCFT). These comprehensive studies confirmed monomer composition as a viable design parameter to engineer RHP-based functional materials as long as the reactivity ratios are accurately determined and the livingness of RHP synthesis is ensured.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-consistent simulation of resistive kink instabilities with runaway electrons

A new fluid model for runaway electron simulation based on fluid description is introduced and implemented in the magnetohydrodynamics code M3D-C1, which includes self-consistent interactions between plasma and runaway electrons. The model utilizes the method of characteristics to solve the continuity equation for the runaway electron density with large convection speed, and uses a modified Boris algorithm for pseudo particle pushing. The model was employed to simulate magnetohydrodynamics instabilities happening in a runaway electron final loss event in the DIII-D tokamak. Nonlinear simulation reveals that a large fraction of runaway electrons get lost to the wall when kink instabilities are excited and form stochastic field lines in the outer region of the plasma. Plasma current converts from runaway electron current to Ohmic current, and get pinched at the magnetic axis. Here, given the good agreement with experiment, the simulation model provides a reliable tool to study macroscopic plasma instabilities in existence of runaway electron current, and can be used to support future studies of runaway electron mitigation strategies in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Calculation of Accurate Harmonic Frequencies of Large Molecules: The Polycyclic Aromatic Hydrocarbons, a Case Study

The vibrational frequencies and infrared intensities of naphthalene neutral and cation are studied at the self-consistent-field (SCF), second-order Moller-Plesset (MP2), and density functional theory (DFT) levels using a variety of one-particle basis sets. Very accurate frequencies can be obtained at the DFT level in conjunction with large basis sets if they are scaled with two factors, one for the C-H stretches and a second for all other modes. We also find remarkably good agreement at the B3LYP/4-31G level using only one scale factor. Unlike the neutral PAHs where all methods do reasonably well for the intensities, only the DFT results are accurate for the PAH cations. The failure of the SCF and MP2 methods is caused by symmetry breaking and an inability to describe charge delocalization. We present several interesting cases of symmetry breaking in this study. An assessment is made as to whether an ensemble of PAH neutrals or cations could account for the unidentified infrared bands observed in many astronomical sources.

Bauschlicher, Charles W., Jr.↗

Multiscale Ecosystem for solving Maxwell-Schrodinger equations of open quantum systems (OpenMS)

Light-matter interactions play an important role in many branches of physics, chemistry, energy, and materials science. In the strong coupling regime, light-matter interactions are able to tune the materials properties via the formation of new quasiparticles, such as plasmons and polaritons. However, theoretical and numerical modeling of the light-matter interaction-mediated processes remain a big challenge because light-matter interactions are fundamentally multiscale and multiphysics problems involving multiple interactions between electrons, nuclei, and photons at different time/length scales. In the current community, light-matter interactions were treated at different levels of theoretical complexity in quantum chemistry and quantum optics. In quantum optics or quantum photonics, the matter is usually simplified as a few-level system, and the light is treated quantum-mechanically. On the other hand, quantum chemistry explores first-principles methods, including both single-particle and many-body-based techniques, to describe the electronic properties of matter in detail. However, light is usually prescribed as a classical electromagnetic field, and the light-matter interaction is taken into account as an external potential via classical approximations. This software is designed to fill current modeling shortcomings by delivering a first-ever scalable multiscale platform for simulating light-matter interactions in realistic electromagnetic environments. The software solves Maxwell and Schrodinger equations self-consistent on the heterogeneous platforms. It implements HF/DFT, TDDFT, and coupled-cluster counterparts for light-matter interactions and adopts modular programming to offload massively parallel algorithms on a large number of CPU/GPUs.

Zhang, Yu↗

Relaxation time approximations in PAOFLOW 2.0

Abstract Regardless of its success, the constant relaxation time approximation has limited validity. Temperature and energy dependent effects are important to match experimental trends even in simple situations. We present the implementation of relaxation time approximation models in the calculation of Boltzmann transport in 2.0 and apply those to model band-structures. In addition, using a self-consistent fitting of the model parameters to experimental conductivity data, we provide a flexible tool to extract scattering rates with high accuracy. We illustrate the approximations using simple models and then apply the method to GaAs, Si, $$\hbox {Mg}_3\hbox {Sb}_2$$ Mg 3 Sb 2 , and $$\hbox {CoSb}_3$$ CoSb 3 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗