Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FUNCTIONALS”

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 451 records · Page 25

Sublinear scaling for time-dependent stochastic density functional theory

A stochastic approach to time-dependent density functional theory is developed for computing the absorption cross section and the random phase approximation (RPA) correlation energy. The core idea of the approach involves time-propagation of a small set of stochastic orbitals which are first projected on the occupied space and then propagated in time according to the time-dependent Kohn-Sham equations. The evolving electron density is exactly represented when the number of random orbitals is infinite, but even a small number (≈16) of such orbitals is enough to obtain meaningful results for absorption spectrum and the RPA correlation energy per electron. We implement the approach for silicon nanocrystals using real-space grids and find that the overall scaling of the algorithm is sublinear with computational time and memory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analytic second derivative of the energy for density functional theory based on the three-body fragment molecular orbital method

Analytic second derivatives of the energy with respect to nuclear coordinates have been developed for spin restricted density functional theory (DFT) based on the fragment molecular orbital method (FMO). The derivations were carried out for the three-body expansion (FMO3), and the two-body expressions can be obtained by neglecting the three-body corrections. Also, the restricted Hartree-Fock (RHF) Hessian for FMO3 can be obtained by neglecting the density-functional related terms. In both the FMO-RHF and FMO-DFT Hessians, certain terms with small magnitudes are neglected for computational efficiency. The accuracy of the FMO-DFT Hessian in terms of the Gibbs free energy is evaluated for a set of polypeptides and water clusters and found to be within 1 kcal/mol of the corresponding full (non-fragmented) ab initio calculation. The FMO-DFT method is also applied to transition states in S{sub N}2 reactions and for the computation of the IR and Raman spectra of a small Trp-cage protein (PDB: 1L2Y). Some computational timing analysis is also presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Native defects in Tl{sub 6}SI{sub 4}: Density functional calculations

Tl{sub 6}SI{sub 4} is a promising room-temperature semiconductor radiation detection material. Here, we report density functional calculations of native defects and dielectric properties of Tl{sub 6}SI{sub 4}. Formation energies and defect levels of native point defects and defect complexes are calculated. Donor-acceptor defect complexes are shown to be abundant in Tl{sub 6}SI{sub 4}. High resistivity can be obtained by Fermi level pinning by native donor and acceptor defects. Deep donors that are detrimental to electron transport are identified and methods to mitigate such problem are discussed. Furthermore, we show that mixed ionic-covalent character of Tl{sub 6}SI{sub 4} gives rise to enhanced Born effective charges and large static dielectric constant, which provides effective screening of charged defects and impurities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Direct assignment of molecular vibrations via normal mode analysis of the neutron dynamic pair distribution function technique

For over a century, vibrational spectroscopy has enhanced the study of materials. Yet, assignment of particular molecular motions to vibrational excitations has relied on indirect methods. Here, we demonstrate that applying group theoretical methods to the dynamic pair distribution function analysis of neutron scattering data provides direct access to the individual atomic displacements responsible for these excitations. Applied to the molecule-based frustrated magnet with a potential magnetic valence-bond state, LiZn{sub 2}Mo{sub 3}O{sub 8}, this approach allows direct assignment of the constrained rotational mode of Mo{sub 3}O{sub 13} clusters and internal modes of MoO{sub 6} polyhedra. We anticipate that coupling this well known data analysis technique with dynamic pair distribution function analysis will have broad application in connecting structural dynamics to physical properties in a wide range of molecular and solid state systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Anatase TiO{sub 2} nanowires functionalized by organic sensitizers for solar cells: A screened Coulomb hybrid density functional study

The adsorption of two different organic molecules cyanidin glucoside (C{sub 21}O{sub 11}H{sub 20}) and TA-St-CA on anatase (101) and (001) nanowires has been investigated using the standard and the range separated hybrid density functional theory calculations. The electronic structures and optical spectra of resulting dye–nanowire combined systems show distinct features for these types of photochromophores. The lowest unoccupied molecular orbital of the natural dye cyanidin glucoside is located below the conduction band of the semiconductor while, in the case of TA-St-CA, it resonates with the states inside the conduction band. The wide-bandgap anatase nanowires can be functionalized for solar cells through electron-hole generation and subsequent charge injection by these dye sensitizers. The intermolecular charge transfer character of Donor-π-Acceptor type dye TA-St-CA is substantially modified by its adsorption on TiO{sub 2} surfaces. Cyanidin glucoside exhibits relatively stronger anchoring on the nanowires through its hydroxyl groups. The atomic structures of dye–nanowire systems re-optimized with the inclusion of nonlinear solvation effects showed that the binding strengths of both dyes remain moderate even in ionic solutions.

36 MATERIALS SCIENCE↗

Optical spectroscopy of molecular junctions: Nonequilibrium Green’s functions perspective

We consider optical spectroscopy of molecular junctions from the quantum transport perspective when radiation field is quantized and optical response of the system is simulated as photon flux. Using exact expressions for photon and electronic fluxes derived within the nonequilibrium Green function (NEGF) methodology and utilizing fourth order diagrammatic perturbation theory (PT) in molecular coupling to radiation field, we perform simulations employing realistic parameters. Results of the simulations are compared to the bare PT which is usually employed in studies on nonlinear optical spectroscopy to classify optical processes. We show that the bare PT violates conservation laws, while flux conserving NEGF formulation mixes optical processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An atomic orbital based real-time time-dependent density functional theory for computing electronic circular dichroism band spectra

One of the challenges of interpreting electronic circular dichroism (ECD) band spectra is that different states may have different rotatory strength signs, determined by their absolute configuration. If the states are closely spaced and opposite in sign, observed transitions may be washed out by nearby states, unlike absorption spectra where transitions are always positive additive. To accurately compute ECD bands, it is necessary to compute a large number of excited states, which may be prohibitively costly if one uses the linear-response time-dependent density functional theory (TDDFT) framework. Here we implement a real-time, atomic-orbital based TDDFT method for computing the entire ECD spectrum simultaneously. The method is advantageous for large systems with a high density of states. In contrast to previous implementations based on real-space grids, the method is variational, independent of nuclear orientation, and does not rely on pseudopotential approximations, making it suitable for computation of chiroptical properties well into the X-ray regime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inequalities for majorizing analytic functions and their applications to rational trigonometric functions and polynomials

New inequalities are established for analytic functions satisfying Meiman's majorization conditions. Estimates for values of and differential inequalities involving rational trigonometric functions with an integer majorant on an interval of length less than the period and with prescribed poles which are symmetrically positioned relative to the real axis, as well as differential inequalities for trigonometric polynomials in some classes, are given as applications. These results improve several theorems due to Meiman, Genchev, Smirnov and Rusak. Bibliography: 27 titles.

97 MATHEMATICS AND COMPUTING↗

AN INITIAL MASS FUNCTION FOR INDIVIDUAL STARS IN GALACTIC DISKS. I. CONSTRAINING THE SHAPE OF THE INITIAL MASS FUNCTION

We derive a semi-empirical galactic initial mass function (IMF) from observational constraints. We assume that the IMF, {psi}(m), is a smooth function of the stellar mass m. The mass dependence of the proposed IMF is determined by five parameters: the low-mass slope {gamma}, the high-mass slope -{Gamma} (taken to be -1.35), the characteristic mass m{sub ch} ({approx} the peak mass of the IMF), and the lower and upper limits on the mass, m{sub l} and m{sub u} (taken to be 0.004 and 120 M{sub sun}, respectively): {psi}(m)dln m {proportional_to} m{sup -}{Gamma}{l_brace}1 - exp [- (m/m{sub ch}){sup {gamma}}+{Gamma}]{r_brace}dln m. The values of {gamma} and m{sub ch} are derived from two integral constraints: (1) the ratio of the number density of stars in the range m = 0.1-0.6 M{sub sun} to that in the range m = 0.6-0.8 M{sub sun} as inferred from the mass distribution of field stars in the local neighborhood and (2) the ratio of the number of stars in the range m = 0.08-1 M{sub sun} to the number of brown dwarfs in the range m = 0.03-0.08 M{sub sun} in young clusters. The IMF satisfying the above constraints is characterized by the parameters {gamma} = 0.51 and m{sub ch} = 0.35 M{sub sun} (which corresponds to a peak mass of 0.27 M{sub sun} ). This IMF agrees quite well with the Chabrier IMF for the entire mass range over which we have compared with data, but predicts significantly more stars with masses <0.03 M{sub sun}; we also compare with other IMFs in current use and give a number of important parameters implied by the IMFs.

79 ASTRONOMY AND ASTROPHYSICS↗

SIMILARITY BETWEEN THE C{sup 18}O (J = 1-0) CORE MASS FUNCTION AND THE INITIAL MASS FUNCTION (IMF) IN THE S140 REGION

We present the results of C{sup 18}O(J = 1-0) mapping observations of a 20' x 18' area in the Lynds 1204 molecular cloud associated with the Sharpless 2-140 (S140) H II region. The C{sup 18}O cube ({alpha}-{delta}-v{sub LSR}) data show that there are three clumps of sizes {approx}1 pc in the region. Two of these have peculiar redshifted velocity components at their edges, which can be interpreted as the results of the interaction between the cloud and the Cepheus Bubble. From the C{sup 18}O cube data, clumpfind identified 123 C{sup 18}O cores, which have mean radius, velocity width in FWHM, and LTE mass of 0.36 {+-} 0.07 pc, 0.37 {+-} 0.09 km s{sup -1}, and 41 {+-} 29 M{sub sun}, respectively. Considering the uncertainty in the C{sup 18}O abundance, all the cores in S140 are most likely to be gravitationally bound. We derived a C{sup 18}O core mass function (CMF), which shows a power-law-like behavior above a turnover at 30 M{sub sun}. The best-fit power-law index of -2.1 {+-} 0.2 is quite consistent with that of the initial mass function (IMF) and the C{sup 18}O CMF in the OMC-1 region as found by Ikeda and Kitamura. Kramer et al. estimated the power-law index of -1.65 in S140 from the C{sup 18}O(J = 2-1) data, which is inconsistent with this study. The C{sup 18}O(J = 2-1) data are spatially limited to the central part of the cloud and are likely to be biased toward high-mass cores, leading to the flatter CMF. Consequently, this study and our previous study strongly support that the power-law form of the IMF has already been determined at the density of {approx}< 10{sup 3}-10{sup 4} cm{sup -3}, traced by the C{sup 18}O(J = 1-0) line.

79 ASTRONOMY AND ASTROPHYSICS↗

THE FERMI LARGE AREA TELESCOPE ON ORBIT: EVENT CLASSIFICATION, INSTRUMENT RESPONSE FUNCTIONS, AND CALIBRATION

The Fermi Large Area Telescope (Fermi-LAT, hereafter LAT), the primary instrument on the Fermi Gamma-ray Space Telescope (Fermi) mission, is an imaging, wide field-of-view, high-energy {gamma}-ray telescope, covering the energy range from 20 MeV to more than 300 GeV. During the first years of the mission, the LAT team has gained considerable insight into the in-flight performance of the instrument. Accordingly, we have updated the analysis used to reduce LAT data for public release as well as the instrument response functions (IRFs), the description of the instrument performance provided for data analysis. In this paper, we describe the effects that motivated these updates. Furthermore, we discuss how we originally derived IRFs from Monte Carlo simulations and later corrected those IRFs for discrepancies observed between flight and simulated data. We also give details of the validations performed using flight data and quantify the residual uncertainties in the IRFs. Finally, we describe techniques the LAT team has developed to propagate those uncertainties into estimates of the systematic errors on common measurements such as fluxes and spectra of astrophysical sources.

79 ASTRONOMY AND ASTROPHYSICS↗

A marked correlation function for constraining modified gravity models

Future large scale structure surveys will provide increasingly tight constraints on our cosmological model. These surveys will report results on the distance scale and growth rate of perturbations through measurements of Baryon Acoustic Oscillations and Redshift-Space Distortions. It is interesting to ask: what further analyses should become routine, so as to test as-yet-unknown models of cosmic acceleration? Models which aim to explain the accelerated expansion rate of the Universe by modifications to General Relativity often invoke screening mechanisms which can imprint a non-standard density dependence on their predictions. This suggests density-dependent clustering as a `generic' constraint. This paper argues that a density-marked correlation function provides a density-dependent statistic which is easy to compute and report and requires minimal additional infrastructure beyond what is routinely available to such survey analyses. We give one realization of this idea and study it using low order perturbation theory. We encourage groups developing modified gravity theories to see whether such statistics provide discriminatory power for their models.

79 ASTRONOMY AND ASTROPHYSICS↗

Density and pair-density scaling for deriving the Euler equation in density-functional and pair-density-functional theory

A link between density and pair density functional theories is presented. Density and pair density scaling are used to derive the Euler equation in both theories. Density scaling provides a constructive way of obtaining approximations for the Pauli potential. The Pauli potential (energy) of the density functional theory is expressed as the difference of the scaled and original exchange-correlation potentials (energies).

74 ATOMIC AND MOLECULAR PHYSICS↗

Generalization of internal density-functional theory and Kohn-Sham scheme to multicomponent self-bound systems, and link with traditional density-functional theory

We generalize the recently developed ''internal'' density-functional theory (DFT) and Kohn-Sham scheme to multicomponent systems. We obtain a general formalism, applicable for the description of multicomponent self-bound systems (such as molecular systems where the nuclei are treated explicitly, atomic nuclei and mixtures of {sup 3}He and {sup 4}He droplets), where the fundamental translational symmetry has been treated correctly. The main difference with traditional DFT is the explicit inclusion of center-of-mass correlations in the functional. A large part of the paper is dedicated to the application to molecular systems, which permits us to clarify the approximations that underly traditional DFT.

74 ATOMIC AND MOLECULAR PHYSICS↗

Antisymmetric wave functions for mixed fermion states and energy convexity

We show how ensembles or mixed states can be described in terms of pure states that for Fermions lead to wave functions that are antisymmetric with respect to interchange of particle coordinates (and spin). The pure states are constructed in an augmented Hilbert space spanned by products of ensemble states projected onto mutually nonoverlapping coordinates that prevents the appearance of interference terms under antisymmetrization. As a demonstration of this new formalism, and under the assumptions of a positive interparticle interaction and a corresponding energy that is extensive in the number of particle pairs (pair extensive), we prove the convexity relation, E v [ N - 1 ] + E v [ N + 1 ] ≤ 2 E v [ N ] , where E v [ N ] denotes the total ground-state energy of N Fermions (electrons) under an external potential v ( r ) .

36 MATERIALS SCIENCE↗

Bypassing the Energy Functional in Density Functional Theory: Direct Calculation of Electronic Energies from Conditional Probability Densities

Density functional calculations can fail for want of an accurate exchange-correlation approximation. The energy can instead be extracted from a sequence of density functional calculations of conditional probabilities (CP DFT). Simple CP approximations yield usefully accurate results for two-electron ions, the hydrogen dimer, and the uniform gas at all temperatures. CP DFT has no self-interaction error for one electron, and correctly dissociates H2, both major challenges. For warm dense matter, classical CP DFT calculations can overcome the convergence problems of Kohn-Sham DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Symmetry-mode analysis for local structure investigations using pair distribution function data

Symmetry-adapted distortion modes provide a natural way of describing distorted structures derived from higher-symmetry parent phases. Structural refinements using symmetry-mode amplitudes as fit variables have been used for at least ten years in Rietveld refinements of the average crystal structure from diffraction data; more recently, this approach has also been used for investigations of the local structure using real-space pair distribution function (PDF) data. Here, the value of performing symmetry-mode fits to PDF data is further demonstrated through the successful application of this method to two topical materials: TiSe2, where a subtle but long-range structural distortion driven by the formation of a charge-density wave is detected, and MnTe, where a large but highly localized structural distortion is characterized in terms of symmetry-lowering displacements of the Te atoms. Here, the analysis is performed using fully open-source code within the DiffPy framework via two packages developed for this work: isopydistort, which provides a scriptable interface to the ISODISTORT web application for group theoretical calculations, and isopytools, which converts the ISODISTORT output into a DiffPy-compatible format for subsequent fitting and analysis. These developments expand the potential impact of symmetry-adapted PDF analysis by enabling high-throughput analysis and removing the need for any commercial software.

36 MATERIALS SCIENCE↗