Search NASA⌕ Search

SEARCH · Search NASA

Results for “polarizability”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Laser pulse induced second- and third-harmonic generation of gold nanorods with real-time time-dependent density functional tight binding (RT-TDDFTB) method

In this paper, we investigate second- and third-harmonic generation processes in Au nanorod systems using the real-time time-dependent density functional tight binding method. Our study focuses on the computation of nonlinear signals based on the time dependent dipole response induced by linearly polarized laser pulses interacting with nanoparticles. We systematically explore the influence of various laser parameters, including pump intensity, duration, frequency, and polarization directions, on harmonic generation. We demonstrate all the results using Au nanorod dimer systems arranged in end-to-end configurations, and disrupting the spatial symmetry of regular single nanorod systems is crucial for second-harmonic generation processes. Furthermore, we study the impact of nanorod lengths, which lead to variable plasmon energies, on harmonic generation, and estimates of polarizabilities and hyper-polarizabilities are provided.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Alignment and dissociation of electronically excited molecular hydrogen with intense laser fields

The dissociation of aligned, electronically excited H 2 (E,F 1 Σ g + ), followed by ionization of the produced H atom, is analyzed via the velocity mapped imaging technique. The dissociation and ionization processes are accomplished, respectively, by a two- and a onephoton absorption from a single 532-nm laser pulse, while the alignment is induced by a separate 1064-nm laser pulse. The velocity of the produced H + photofragments shows a weak perpendicular alignment at low alignment laser field values, evolving to strongly parallel for larger fields. We modeled this alignment behavior with a simple two-state model involving the Stark mixing of the initially-prepared J = 0 with the J = 2 rotational state. This model is able to reproduce all of the observed angular distribution, and permits us to extract from the fit the polarizability anisotropy of H 2 (E,F) electronic state. We determine this value to be (3.7 ± 1.2) x 10 3 a.u. As this value is extremely large in comparison to what one would expect from the pure H 2 (E,F) electronic state, we hypothesize that this value comes from the 1064-nm laser beam mixing nearby electronic states with the initially laser prepared (E,F) state generating a mixed state (EF**) with an extremely large polarizability anisotropy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Fermi-Löwdin-orbital self-interaction correction using the optimized-effective-potential method within the Krieger-Li-Iafrate approximation

Perdew-Zunger self-interaction correction (PZ-SIC) offers a route to remove self interaction errors on an orbital-by-orbital basis. A recent formulation of PZ-SIC by Pederson, Ruzsinszky and Perdew proposes restricting the unitary transformation to localized orbitals called Fermi-Löwdin orbitals. This formulation, called the FLOSIC method, simplifies PZ-SIC calculations and was implemented self-consistently using a Jacobi-like (FLOSIC-Jacobi) iteration scheme. In this work we implement the FLOSIC approach using the Krieger-Li-Iafrate (KLI) approximation to the optimized effective potential (OEP). We compare the results of present FLOSIC-KLI approach with FLOSIC-Jacobi scheme for atomic energies, atomization energies, ionization energies, barrier heights, polarizability of chains of hydrogen molecules etc. to validate the FLOSIC-KLI approach. The FLOSICKLI approach, which is within the realm of Kohn-Sham theory, predicts smaller energy gaps between frontier orbitals due to the lowering of eigenvalues of the lowest unoccupied orbitals. Results show that atomic energies, atomization energies, ionization energy as an absolute of highest occupied orbital eigenvalue, and polarizability of chains of hydrogen molecules between the two methods agree within 2%. Lastly the FLOSIC-KLI approach is used to determine the vertical ionization energies of water clusters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Full versus quasiparticle self-consistency in vertex-corrected GW approaches

Here, using seven semiconductors/insulators with band gaps covering the range from 1 eV to 10 eV we systematically explore the performance of two different variants of self-consistency associated with famous Hedin's system of equations: the full self-consistency and the so called quasi-particle approximation to it. The pros and cons of these two variants of self-consistency are sufficiently well documented in literature for the simplest GW approximation to the Hedin's equations. Our study, therefore, aims primarily at the level of theory beyond GW approximation, i.e. at the level of theory which includes vertex corrections. Whereas quasi-particle self-consistency has certain advantages at GW level (well known fact), the situation becomes quite different when vertex corrections are included. In the variant with full self-consistency, vertex corrections (both for polarizability and for self energy) systematically reduce the calculated band gaps making them closer to the experimental values. In the variant with quasi-particle self-consistency, however, an inclusion of the same diagrams has considerably larger effect and calculated band gaps become severely underestimated. Different effect of vertex corrections in two variants of self-consistency can be related to the Z-factor cancellation which plays positive role in quasi-particle self-consistency at GW level of theory but appears to be destructive for the quasi-particle approximation when higher order diagrams are included. Second result of our study is that we were able to reproduce the results obtained with the Questaal code using our FlapwMBPT code when the same variant of self-consistency (quasi-particle) and the same level of vertex corrections (for polarizability only, static approximation for screened interaction, and Tamm-Dancoff approximation for the Bethe-Salpeter equation) are used.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Third-order Hall effect in the surface states of a topological insulator

Time reversal and inversion symmetric materials fail to yield linear and nonlinear responses since they possess net zero Berry curvature. However, higher-order Hall response can be generated in these systems upon constraining the crystalline symmetries. Motivated by the recently discovered third-order Hall (TOH) response mediated by Berry connection polarizability, namely, the variation of the Berry connection with respect to an applied electric field, here, we investigate the existence of such a Hall effect in the surface states of hexagonal warped topological insulators (e.g., Bi 2 Te 3 ) under the application of only electric field. Using the semiclassical Boltzmann formalism, we investigate the effect of tilt and hexagonal warping on the Berry connection polarizability tensor and, consequently, the TOH effect, provided the Dirac cone remains gapless. We find that the magnitude of the response increases significantly with increasing tilt strength and warping, and therefore, they can provide the tunability of this effect. In addition, we also explore the effect of chemical doping on the TOH response in this system. Interestingly, we show, based on the symmetry analysis, that the TOH can be the leading-order response in this system, which can directly be verified in experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Nuclear multipole responses from chiral effective field theory interactions

We probe nuclear multipole resonances in the framework of the random-phase approximation by using the interaction obtained from the chiral effective field theory. The three-nucleon force is included in a form of the in-medium two-nucleon interaction, which was derived from the chiral three-nucleon force. The isoscalar monopole, isoscalar dipole, isovector dipole, and isoscalar quadrupole resonances of the closed-shell 56,68,78 Ni have been investigated. The calculations reasonably reproduce the experimental multipole resonances of 56,68 Ni, and well describe the pygmy dipole resonance and dipole polarizability measured in 68 Ni. The multipole resonances of 78 Ni, including pygmy dipole resonance and dipole polarizability, are predicted. The detailed effects of the tensor force and three-body force are analyzed by dissecting the chiral interaction. In this work, we find that in general the tensor force effect on electric giant resonances is not as significant as the effect from the three-body force, although the tensor force provides more than half of the binding energy. The effect from three-body force is strong in light nuclei. Particularly, three-body force is crucial for the formation of the pygmy resonance in calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Anisotropic X-Ray Scattering of Transiently Oriented Water

We study the structural dynamics of liquid water by time-resolved anisotropic x-ray scattering under the optical Kerr effect condition. In this way, we can separate the anisotropic scattering decay of 160 fs from the delayed temperature increase of ~ 0.1 K occurring at 1 ps and quantify transient changes in the O-O pair distribution function. Polarizable molecular dynamics simulations reproduce well the experiment, indicating transient alignment of molecules along the electric field, which shortens the nearest-neighbor distances. In addition, analysis of the simulated water local structure provides evidence that two hypothesized fluctuating water configurations exhibit different polarizability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Combined Theoretical Analysis of the Parity-Violating Asymmetry for 48 Ca and 208 Pb

The recent experimental determination of the parity violating asymmetry A PV in 48 Ca and 208 Pb at Jefferson Lab is important for our understanding on how neutrons and protons arrange themselves inside the atomic nucleus. To better understand the impact of these measurements, we present a rigorous theoretical investigation of A PV in 48 Ca and 208 Pb and assess the associated uncertainties. Here, we complement our study by inspecting the static electric dipole polarizability in these nuclei. The analysis is carried out within nuclear energy density functional theory with quantified input. We conclude that the simultaneous accurate description of A PV in 48 Ca and 208 Pb cannot be achieved by our models that accommodate a pool of global nuclear properties, such as masses and charge radii, throughout the nuclear chart, and describe—within one standard deviation—the experimental dipole polarizabilities in these nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Investigation of the effect of local atomic polarization on field-enhanced ion transport in insulating binary oxides: The case of CeO 2

Over the course of more than 70 years, several theoretical models for electric field-enhanced ion diffusion in crystalline insulators have been presented. However, there is no assessment of the validity of these models despite the urge to have a correct description of field-enhanced ion diffusion in several emerging technologies. Herein, an assessment of five models was carried out by computing the field-dependent migration enthalpy ΔH mig (E) for the oxide ion in CeO 2 . The input to these models is the zero-field and zero-temperature migration pathway and/or the activation barrier obtained from an interatomic potential. ΔH mig (E) from these models was compared with reference values obtained from a set of classical molecular-dynamics simulations in the temperature range of 1000 ≤ T ≤ 1600 K and the electric field range of 0 ≤ E ≤ 30 MV/cm. It is revealed that the most successful theoretical models are those that consider local polarization effects induced in the diffusing ion itself and its vicinity. However, all the models did not account for the effect of the electric field and finite temperature on the local polarization, leading to discrepancies between the predicted and the reference ΔH mig (E), particularly at high fields. By comparing results from a rigid-ion potential and a polarizable core-shell potential, it is concluded that an intraionic polarization degree of freedom in the polarizable potential is an important factor in predicting ΔH mig (E), regardless of the field-enhanced diffusion model particularly for low to moderate field. Finally, future more accurate treatments should consider the field- and temperature-dependent interionic and intraionic local polarization effects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure and dynamics of open-shell nuclei from spherical coupled-cluster theory

We extend the spherical coupled-cluster ab initio method for open-shell nuclei where two nucleons are removed from a shell subclosure. Following the recent implementation of the two-particle-attached approach [Phys. Rev. C 110, 044306 (2024)], we focus on the two-particle-removed method. Using the equations-of-motion framework, we address both nuclear structure and dipole response functions by coupling coupled-cluster theory with the Lorentz integral transform technique. We perform calculations using chiral interactions, including three-nucleon forces, and estimate many-body uncertainties by comparing different coupled-cluster truncation schemes. Here, we validate our approach by studying ground-state energies, excited states, and electric dipole polarizabilities in the oxygen and calcium isotopic chains. For binding energies and selected low-lying excited states, we achieve an accuracy comparable to that of the established closed-shell coupled-cluster theory and generally agree with experiment. Finally, we underestimate experimental data for electric dipole polarizabilities, particularly in calcium isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Ultrahigh capacitive energy density in ion-bombarded relaxor ferroelectric films

Dielectric capacitors can store and release electric energy at ultrafast rates and are extensively studied for applications in electronics and electric power systems. Among various candidates, thin films based on relaxor ferroelectrics, a special kind of ferroelectric with nanometer-sized domains, have attracted special attention because of their high energy densities and efficiencies. We show that high-energy ion bombardment improves the energy storage performance of relaxor ferroelectric thin films. Intrinsic point defects created by ion bombardment reduce leakage, delay low-field polarization saturation, enhance high-field polarizability, and improve breakdown strength. We demonstrate energy storage densities as high as ~133 joules per cubic centimeter with efficiencies exceeding 75%. Deterministic control of defects by means of postsynthesis processing methods such as ion bombardment can be used to overcome the trade-off between high polarizability and breakdown strength.

Kim, Jieun↗
DOE OSTI2020

Simplified Ab Initio Molecular Dynamics-Based Raman Spectral Simulations

We describe a simplified approach to simulating Raman spectra from ab initio molecular dynamics (AIMD) calculations. The protocol relies on on-the-fly calculations of approximate molecular polarizabilities using the well-known sum over orbitals (as opposed to states) method. This approach bypasses the more accurate but computationally expensive approach to calculating molecular polarizabilities along AIMD trajectories, i.e., solving the coupled perturbed Hartree–Fock/Kohn–Sham equations. We demonstrate the advantages and limitations of our method through a few case studies targeting molecular systems of interest to surface- and/or tip-enhanced Raman spectroscopy practitioners.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Exploration of a Novel Technique for Waste Heat Recovery Through Molecular Dynamics: Influence of Wettability and Electric Field on Water and Water-Based Nanofluids

Most of the energy produced globally comes by way of a heat engine. The Carnot principle places a limit as to how thermodynamically efficient a heat engine can be. There is no heat engine that can be 100% thermodynamically efficient and as such a substantial proportion of all heat supplied to a heat engine is lost as waste heat. Waste heat therefore is a large energy source ready to be properly utilized. Herein, a novel approach for converting waste heat to electricity is discussed. It involves the use of the liquid to vapor phase change of a material dielectric (water) or electrolyte (nanofluid) in the embodiment of a capacitor for direct thermal to electrostatic energy conversion. While this method of waste heat recovery could potentially be added to the ever expanding portfolio of energy conversion techniques, a number of aspects must be addressed before it can be brought into practice. Water was seen as an ideal dielectric phase change material given its high relative permittivity ratio when in the liquid form as compared to its vapor form. However, given its short voltage holdoff time the phase change of water would need to occur rapidly. This brings up concerns of explosive boiling. Herein, molecular dynamics analysis into the explosive boiling behavior of thin water films gave more insight into how the interaction between the surface and liquid affected explosive boiling onset time. A Lennard-Jones potential with one interaction site and a Morse potential with three interaction sites between water and solid substrate were used. It was found generally that a stronger interaction between water film and substrate led to faster explosive boiling onset times but an increase in the number of interaction sites delayed explosive boiling, even at the same wettability (contact angle). Understanding changes in the density and enthalpy of vaporization of a liquid dielectric such as water in the presence of an electric field is of importance due to the electrostatic nature of the waste heat conversion method under consideration. Specifically, if both density and enthalpy of vaporization are increased, the thermodynamic efficiency of the waste heat conversion method under consideration is decreased. Electric field effects are explored herein via molecular dynamics using two water models, the TIP4P-Ew and SWM4-NDP. The SWM4-NDP model is polarizable while the TIP4P-Ew model is not, which allows for a determination of the importance of model polarizability (i.e. variation in water model dipole moment) on these two properties of water when subjected to an electric field. Herein it was found that both water models respond similarly in terms of density and vaporization enthalpy variance upon the introduction of an electric field. Comparison was also made to the pressure induced by the electric field (electrostriction pressure) by way of a density comparison and it was found that the predicted electrostriction pressure overestimates the pressure experienced by water. Water by itself has a high enthalpy of vaporization, which limits the efficiency of the newly proposed conversion method. Research both experimental and through simulation has shown that the vaporization enthalpy of nanofluids can be engineered via nanoparticle size and material selection. An avenue less explored is manipulating the enthalpy of vaporization by altering the interaction strength between the nanoparticles and the base fluid. In practice this could be achieved through the addition of coatings to the nanoparticles to alter their wettability to the base fluid. This was explored by using a Lennard-Jones potential and Morse potential to model the interaction between base fluid (water) and the nanoparticle. For nanoparticles 2nm in diameter and at weight percentages up to 6%, the change in vaporization enthalpy due to alterations of the interaction strength between nanoparticle and base fluid was not significant (less than a 1% difference) when compared to the effect of altering the weight percentage of nanoparticles in the nanofluid or introducing an electric field. However, the effect of wettability may still become important at other nanoparticle concentrations and sizes. In all, the studies presented here further the understanding of phase change and thermodynamic properties of water and water based nanofluids under an electrostatic field which will help inform the development of a novel approach to waste heat conversion. The reduction of waste heat will improve energy sustainability outlooks.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

Schiff Moment of 129 Xe and 199 Hg [Slides]

Using Chicoma in 2025, we advanced DFT-based calculations needed to connect fundamental T violating interactions to measurable Schiff moments in heavy nuclei. We implemented DFT for polarized odd nuclear systems, a first step toward treating the time-reversal-breaking mean fields required for Schiff moment calculations. As a validation test, we computed polarizabilities using a perturbation with the dipole operator (see figure to the right) for the odd systems shown in Table 1; while no odd-system effect is expected for polarizabilities, this capability is essential for Schiff moments. Initial mean-field solutions for 129 Xe Schiff moment calculations were also generated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Porous Organic Cages CC3 and CC2 as Adsorbents for the Separation of Carbon Dioxide from Nitrogen and Hydrogen

The selective capture of carbon dioxide over nitrogen and hydrogen is of great industrial interest in flue gas and hydrogen purification respectively. Microporous adsorbents are highly suitable materials to preferentially adsorb gases. Here, in particular, Porous organic cages (POCs) with tunable hierarchically ordered micropores, high surface area, and thermodynamic affinity for CO 2 make them appealing candidates for these applications. Herein, we demonstrate that two prototypical POCs denoted as CC3, and CC2 with limiting pore aperture of 3.6 Å can selectively separate CO 2 from N 2 , and H 2 . For CC3 adsorption selectivities as high as ~ 8 and ~20 for CO 2 /N 2 and CO 2 /H 2 respectively were observed. For CC2 adsorption selectivities as high as ~ 9 and ~35 for CO 2 /N 2 and CO 2 /H 2 respectively were observed. Interestingly, the adsorption selectivity of the studied gases correlated linearly with polarizability selectivity.

Carbon Dioxide↗
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 OSTI2022

Multipolar Raman Scattering vs Interfacial Nanochemistry: Case of 4-Mercaptopyridine on Gold

Caution needs to be exercised in associating changes in plasmon-enhanced Raman spectra with chemical transformations. This is demonstrated through a detailed analysis of tip-enhanced Raman (TER) scattering from 4-mercaptopyridine (MPY) on gold. The substrate used consists of gold nanoplates atop a gold surface featuring heterogeneous grooves, all coated with a monolayer of MPY. The brightest spectra across the substrate exhibit features that can only be recovered by considering the generalized polarizability of oriented MPY molecules. Finally, the complex TER spectra we observe do not mark interfacial chemistry, but rather, multipolar TER scattering driven by local field gradients.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗