Science.gov

Sample records for effective potential calculation

  1. {alpha}-decay calculations of heavy and superheavy nuclei using effective mean-field potentials

    SciTech Connect

    Pei, J. C.; Lin, Z. J.; Xu, F. R.; Zhao, E. G.

    2007-10-15

    Using an effective potential that is based on the Skyrme-Hartree-Fock mean-field model, systematic {alpha}-decay properties of even-even heavy and superheavy nuclei have been investigated. Calculations do not raise any adjustable parameter. The obtained {alpha}-decay half-lives agree reasonably well with experimental data. The characteristics of the effective potential and the deformation effect on the {alpha} decay are discussed.

  2. Use of Relativistic Effective Core Potentials in the Calculation of Electron-Impact Ionization Cross Sections

    NASA Technical Reports Server (NTRS)

    Huo, Winifred M.; Kim, Yong-Ki

    1999-01-01

    Based on the Binary-Encounter-Bethe (BEB) model, the advantage of using relativistic effective core potentials (RECP) in the calculation of total ionization cross sections of heavy atoms or molecules containing heavy atoms is discussed. Numerical examples for Ar, Kr, Xe, and WF6 are presented.

  3. Exchange-only optimized effective potential calculation of excited state spectra for He and Be atoms.

    SciTech Connect

    Desjarlais, Michael Paul; Muller, Richard Partain

    2006-02-01

    The optimized effective potential (OEP) method allows orbital-dependent functionals to be used in density functional theory (DFT), which, in particular, allows exact exchange formulations of the exchange energy to be used in DFT calculations. Because the exact exchange is inherently self-interaction correcting, the resulting OEP calculations have been found to yield superior band-gaps for condensed-phase systems. Here we apply these methods to the isolated atoms He and Be, and compare to high quality experiments and calculations to demonstrate that the orbital energies accurately reproduce the excited state spectrum for these species. These results suggest that coupling the exchange-only OEP calculations with proper (orbital-dependent or other) correlation functions might allow quantitative accuracy from DFT calculations.

  4. Effect of deformation on the calculated half-lives of cluster emission using a proximity potential

    NASA Astrophysics Data System (ADS)

    Nazarzadeh, P.; Mohebali, M.

    2016-06-01

    The half-lives of deformed cluster emission have been calculated by using the proximity potential Christensen-Winther potential form (1976) and compared to the experimental data and theoretical results of the liquid drop model. Also the calculated results have been compared to the results obtained by a universal form of the proximity potential for spherical cluster emission. The results show reasonable agreement with the experimental data.

  5. Temperature dependent effective potential method for accurate free energy calculations of solids

    NASA Astrophysics Data System (ADS)

    Hellman, Olle; Steneteg, Peter; Abrikosov, I. A.; Simak, S. I.

    2013-03-01

    We have developed a thorough and accurate method of determining anharmonic free energies, the temperature dependent effective potential technique (TDEP). It is based on ab initio molecular dynamics followed by a mapping onto a model Hamiltonian that describes the lattice dynamics. The formalism and the numerical aspects of the technique are described in detail. A number of practical examples are given, and results are presented, which confirm the usefulness of TDEP within ab initio and classical molecular dynamics frameworks. In particular, we examine from first principles the behavior of force constants upon the dynamical stabilization of the body centered phase of Zr, and show that they become more localized. We also calculate the phase diagram for 4He modeled with the Aziz potential and obtain results which are in favorable agreement both with respect to experiment and established techniques.

  6. Nucleon-nucleon potential calculated for the cloudy bag and related to effective OBE parameters

    NASA Astrophysics Data System (ADS)

    Kahler, Richard

    1995-05-01

    The nucleon nucleon potential is calculated in the context of the cloudy bag model. The one pion exchange diagram is used to determine the pion quark coupling. Two pion exchange diagrams, box and crossed box, including the delta isobar, produce additional central, spin-spin, and tensor couplings. Using non-relativistic approximations to the OBE model potentials, the two pion exchange contributions are related to the masses and couplings of σ, ω, δ, and ρ mesons. It is found that the two pion exchange generates essentially all of the σ contribution required by the OBE model, that there is a large δ contribution at a lower mass than expected, and that the ω and ρ contributions, while significant, are not able to provide what is required by the OBE model.

  7. Properties of metastable alkaline-earth-metal atoms calculated using an accurate effective core potential

    SciTech Connect

    Santra, Robin; Christ, Kevin V.; Greene, Chris H.

    2004-04-01

    The first three electronically excited states in the alkaline-earth-metal atoms magnesium, calcium, and strontium comprise the (nsnp){sup 3}P{sub J}{sup o}(J=0,1,2) fine-structure manifold. All three states are metastable and are of interest for optical atomic clocks as well as for cold-collision physics. An efficient technique--based on a physically motivated potential that models the presence of the ionic core--is employed to solve the Schroedinger equation for the two-electron valence shell. In this way, radiative lifetimes, laser-induced clock shifts, and long-range interaction parameters are calculated for metastable Mg, Ca, and Sr.

  8. On the use of effective core potentials in the calculation of magnetic properties, such as magnetizabilites and magnetic shieldings.

    PubMed

    van Wüllen, Christoph

    2012-03-21

    State-of-the art effective core potentials (ECPs) that replace electrons of inner atomic cores involve non-local potentials. If such an effective core potential is added to the Hamiltonian of a system in a magnetic field, the resulting Hamiltonian is not gauge invariant. This means, magnetic properties such as magnetisabilities and magnetic shieldings (or magnetic susceptibilities and nuclear magnetic resonance chemical shifts) calculated with different gauge origins are different even for exact solutions of the Schrödinger equation. It is possible to restore gauge invariance of the Hamiltonian by adding magnetic field dependent terms arising from the effective core potential. Numerical calculations on atomic and diatomic model systems (potassium mono-cation and potassium dimer) clearly demonstrate that the standard effective core potential Hamiltonian violates gauge invariance, and this affects the calculation of magnetisabilities more strongly than the calculation of magnetic shieldings. The modified magnetic field dependent effective core potential Hamiltonian is gauge invariant, and therefore it is the correct starting point for distributed gauge origin methods. The formalism for gauge including atomic orbitals (GIAO) and individual gauge for localized orbitals methods is worked out. ECP GIAO results for the potassium dimer are presented. The new method performs much better than a previous ECP GIAO implementation that did not account for the non-locality of the potential. For magnetic shieldings, deviations are clearly seen, but they amount to few ppm only. For magnetisabilities, our new ECP GIAO implementation is a major improvement, as demonstrated by the comparison of all-electron and ECP results. PMID:22443751

  9. Calculate and Plot Complex Potential

    Energy Science and Technology Software Center (ESTSC)

    1998-05-05

    SOLUPLOT is a program designed to calculate and plot complex potential, pH diagrams and log oxygen activity, pH diagrams for aqueous chemical syatems, considering speciation of ligands, from free energy and thermodynamic activity data. These diagrams, commonly referred to as Eh-pH and ao2-pH diagrams, respectively, define areas of predominance in Eh-pH diagrams or ao2-pH space for chemical species of a chemical system at equilibrium. Over an area of predominance, one predominant species is at greatermore » activity than the other species of the system considered. The diagram axes, pH (a measure of hydrogen ion activity) and either Eh or log ao2 (measures of a tendency toward either oxidation or reduction) , are paremeters commonly applied in describing the chemistry of aqueous systems.« less

  10. Molecular wave function and effective adiabatic potentials calculated by extended multi-configuration time-dependent Hartree-Fock method

    NASA Astrophysics Data System (ADS)

    Kato, Tsuyoshi; Ide, Yoshihiro; Yamanouchi, Kaoru

    2015-12-01

    We first calculate the ground-state molecular wave function of 1D model H2 molecule by solving the coupled equations of motion formulated in the extended multi-configuration time-dependent Hartree-Fock (MCTDHF) method by the imaginary time propagation. From the comparisons with the results obtained by the Born-Huang (BH) expansion method as well as with the exact wave function, we observe that the memory size required in the extended MCTDHF method is about two orders of magnitude smaller than in the BH expansion method to achieve the same accuracy for the total energy. Second, in order to provide a theoretical means to understand dynamical behavior of the wave function, we propose to define effective adiabatic potential functions and compare them with the conventional adiabatic electronic potentials, although the notion of the adiabatic potentials is not used in the extended MCTDHF approach. From the comparison, we conclude that by calculating the effective potentials we may be able to predict the energy differences among electronic states even for a time-dependent system, e.g., time-dependent excitation energies, which would be difficult to be estimated within the BH expansion approach.

  11. Molecular wave function and effective adiabatic potentials calculated by extended multi-configuration time-dependent Hartree-Fock method

    SciTech Connect

    Kato, Tsuyoshi; Ide, Yoshihiro; Yamanouchi, Kaoru

    2015-12-31

    We first calculate the ground-state molecular wave function of 1D model H{sub 2} molecule by solving the coupled equations of motion formulated in the extended multi-configuration time-dependent Hartree-Fock (MCTDHF) method by the imaginary time propagation. From the comparisons with the results obtained by the Born-Huang (BH) expansion method as well as with the exact wave function, we observe that the memory size required in the extended MCTDHF method is about two orders of magnitude smaller than in the BH expansion method to achieve the same accuracy for the total energy. Second, in order to provide a theoretical means to understand dynamical behavior of the wave function, we propose to define effective adiabatic potential functions and compare them with the conventional adiabatic electronic potentials, although the notion of the adiabatic potentials is not used in the extended MCTDHF approach. From the comparison, we conclude that by calculating the effective potentials we may be able to predict the energy differences among electronic states even for a time-dependent system, e.g., time-dependent excitation energies, which would be difficult to be estimated within the BH expansion approach.

  12. Efficient methods for including quantum effects in Monte Carlo calculations of large systems: Extension of the displaced points path integral method and other effective potential methods to calculate properties and distributions

    NASA Astrophysics Data System (ADS)

    Mielke, Steven L.; Dinpajooh, Mohammadhasan; Siepmann, J. Ilja; Truhlar, Donald G.

    2013-01-01

    We present a procedure to calculate ensemble averages, thermodynamic derivatives, and coordinate distributions by effective classical potential methods. In particular, we consider the displaced-points path integral (DPPI) method, which yields exact quantal partition functions and ensemble averages for a harmonic potential and approximate quantal ones for general potentials, and we discuss the implementation of the new procedure in two Monte Carlo simulation codes, one that uses uncorrelated samples to calculate absolute free energies, and another that employs Metropolis sampling to calculate relative free energies. The results of the new DPPI method are compared to those from accurate path integral calculations as well as to results of two other effective classical potential schemes for the case of an isolated water molecule. In addition to the partition function, we consider the heat capacity and expectation values of the energy, the potential energy, the bond angle, and the OH distance. We also consider coordinate distributions. The DPPI scheme performs best among the three effective potential schemes considered and achieves very good accuracy for all of the properties considered. A key advantage of the effective potential schemes is that they display much lower statistical sampling variances than those for accurate path integral calculations. The method presented here shows great promise for including quantum effects in calculations on large systems.

  13. A self-consistent, microenvironment modulated screened coulomb potential approximation to calculate pH-dependent electrostatic effects in proteins.

    PubMed

    Mehler, E L; Guarnieri, F

    1999-07-01

    An improved approach is presented for calculating pH-dependent electrostatic effects in proteins using sigmoidally screened Coulomb potentials (SCP). It is hypothesized that a key determinant of seemingly aberrant behavior in pKa shifts is due to the properties of the unique microenvironment around each residue. To help demonstrate this proposal, an approach is developed to characterize the microenvironments using the local hydrophobicity/hydrophilicity around each residue of the protein. The quantitative characterization of the microenvironments shows that the protein is a complex mosaic of differing dielectric regions that provides a physical basis for modifying the dielectric screening functions: in more hydrophobic microenvironments the screening decreases whereas the converse applies to more hydrophilic regions. The approach was applied to seven proteins providing more than 100 measured pKa values and yielded a root mean square deviation of 0.5 between calculated and experimental values. The incorporation of the local hydrophobicity characteristics into the algorithm allowed the resolution of some of the more intractable problems in the calculation of pKa. Thus, the divergent shifts of the pKa of Glu-35 and Asp-66 in hen egg white lysozyme, which are both about 90% buried, was correctly predicted. Mechanistically, the divergence occurs because Glu-35 is in a hydrophobic microenvironment, while Asp-66 is in a hydrophilic microenvironment. Furthermore, because the calculation of the microenvironmental effects takes very little CPU time, the computational speed of the SCP formulation is conserved. Finally, results from different crystal structures of a given protein were compared, and it is shown that the reliability of the calculated pKa values is sufficient to allow identification of conformations that may be more relevant for the solution structure. PMID:10388736

  14. Towards a spectroscopically accurate set of potentials for heavy hydride laser cooling candidates: Effective core potential calculations of BaH

    NASA Astrophysics Data System (ADS)

    Moore, Keith; McLaughlin, Brendan M.; Lane, Ian C.

    2016-04-01

    BaH (and its isotopomers) is an attractive molecular candidate for laser cooling to ultracold temperatures and a potential precursor for the production of ultracold gases of hydrogen and deuterium. The theoretical challenge is to simulate the laser cooling cycle as reliably as possible and this paper addresses the generation of a highly accurate ab initio 2Σ+ potential for such studies. The performance of various basis sets within the multi-reference configuration-interaction (MRCI) approximation with the Davidson correction is tested and taken to the Complete Basis Set (CBS) limit. It is shown that the calculated molecular constants using a 46 electron effective core-potential and even-tempered augmented polarized core-valence basis sets (aug-pCVnZ-PP, n = 4 and 5) but only including three active electrons in the MRCI calculation are in excellent agreement with the available experimental values. The predicted dissociation energy De for the X2Σ+ state (extrapolated to the CBS limit) is 16 895.12 cm-1 (2.094 eV), which agrees within 0.1% of a revised experimental value of <16 910.6 cm-1, while the calculated re is within 0.03 pm of the experimental result.

  15. Model potential calculations of lithium transitions.

    NASA Technical Reports Server (NTRS)

    Caves, T. C.; Dalgarno, A.

    1972-01-01

    Semi-empirical potentials are constructed that have eigenvalues close in magnitude to the binding energies of the valence electron in lithium. The potentials include the long range polarization force between the electron and the core. The corresponding eigenfunctions are used to calculate dynamic polarizabilities, discrete oscillator strengths, photoionization cross sections and radiative recombination coefficients. A consistent application of the theory imposes a modification on the transition operator, but its effects are small for lithium. The method presented can be regarded as a numerical generalization of the widely used Coulomb approximation.

  16. Global Calculations Using Potential Functions and Effective Hamiltonian Models for Vibration-Rotation Spectroscopy and Dynamics of Small Polyatomic Molecules

    NASA Astrophysics Data System (ADS)

    Tuyterev, Vladimir

    2010-06-01

    It has become increasingly common to use accurate potential energy surfaces and dipole moment surfaces for predictions and assignment of high-resolution vibration-rotation molecular spectra. These surfaces are obtained either from high-level ab initio electronic structure calculations or from a direct fit to experimental spectroscopic data. The talk will continue a discussion of some recent advances in the domain of the "potentiology". The role of basis extrapolations, of the Born-Oppenheimer breakdown corrections , in particular for very highly excited vibration states will be considered. As effective polyad Hamiltonians and band transition moment operators are still widely used for data reductions in high-resolutions molecular spectroscopy, experimental spectra analyses invoke a need for accurate methods of building physically meaningful effective models from ab initio surfaces. This involves predictions for various spectroscopic constants, including vibration dependence of rotational and centrifugal distortion and resonance coupling parameters. Topics planned for discussion include: high-order Contact Transformations of rovibrational Hamiltonians and of the dipole moment for small polyatomic molecules; convergence issues; the role of the anharmonicity in a potential energy function and of resonance couplings on the normal mode mixing and on vib-rot assignments with application to high energy vibration levels of SO_2 and to ozone near the dissociation limit; intensity anomalies in H_2S / HDS / D_2S spectra, relation of the shape of ab initio dipole moment surfaces with a "mystery" of nearly vanishing symmetry allowed bands. A full account for symmetry properties requires efficient theoretical tools for transformations of molecular Hamiltonians such as irreducible tensor formalism, applications using phosphine and methane potentials will be discussed. Both potential functions and effective polyad Hamiltonians allow studying changes in quasi-classical vibration

  17. Effect of charge distribution over a chlorophyll dimer on the redox potential of P680 in photosystem II as studied by density functional theory calculations.

    PubMed

    Takahashi, Ryouta; Hasegawa, Koji; Noguchi, Takumi

    2008-06-17

    The effect of charge distribution over a chlorophyll dimer on the redox potential of P680 in photosystem II was studied by density functional theory calculations using the P680 coordinates in the X-ray structure. From the calculated ionization potentials of the dimer and the monomeric constituents, the decrease in the redox potential by charge delocalization over the dimer was estimated to be approximately 140 mV. Such charge delocalization was previously observed in the isolated D1-D2-Cyt b 559 complexes, whereas the charge was primarily localized on P D1 in the core complexes. The calculated potential decrease of approximately 140 mV can explain the inhibition of Y Z oxidation in the former complexes and in turn implies that the charge localization on P D1 upon formation of the core complex increases the P680 potential to the level necessary for water oxidation. PMID:18500822

  18. The effect of the time interval used to calculate mean wind velocity on the calculated drift potential, relative drift potential, and resultant drift direction for sands from three deserts in northern China

    NASA Astrophysics Data System (ADS)

    Zhang, Zhengcai; Dong, Zhibao; Zhao, Aiguo

    2016-01-01

    Wind is the power behind many erosion processes and is responsible for many of the characteristics of arid zone geomorphology. Wind velocity is a key factor in determining the potential sand transport, but the nature of the wind velocity data can strongly affect assessments of the risk of blowing sand. In this study, we obtained real-time wind velocity data in a region of the Tengger Desert with shifting sands, in the Badain Jaran Desert, and in the Madoi desertification land, with the data obtained at 1-min intervals, and used the data to determine the influence of how the wind velocity was calculated (mean versus mid-point values and the averaging time used to calculate these values) on sand drift potential. In the three regions, for both the mean and the mid-point wind velocities, the estimated drift potential decreased with increasing averaging time. The relationships between velocities calculated using the different averaging time intervals and the value calculated using a 1-min interval could be expressed as linear functions. The drift potential calculated using the mid-point wind velocity was larger than that calculated using the mean wind velocity.

  19. Relativistic density-functional theory with the optimized effective potential and self-interaction correction: Application to atomic structure calculations (Z=2-106)

    NASA Astrophysics Data System (ADS)

    Tong, Xiao-Min; Chu, Shih-I.

    1998-02-01

    We present a self-interaction-free relativistic density-functional theory (DFT). The theory is based on the extension of our recent nonrelativistic DFT treatment with optimized effective potential (OEP) and self-interaction correction (SIC) [Phys. Rev. A 55, 3406 (1997)] to the relativistic domain. Such a relativistic OEP-SIC procedure yields an orbital-independent single-particle local potential with proper long-range Coulombic (-1/r) behavior. The method is applied to the ground-state energy calculations for atoms with Z=2-106. A comparison with the corresponding nonrelativistic OEP-SIC calculations and other relativistic calculations is made. It is shown that the ionization potentials (obtained from the highest occupied orbital energies) and individual orbital binding energies determined by the present relativistic OEP-SIC method agree well with the experimental data across the Periodic Table.

  20. Relevance of Hydrodynamic Effects for the Calculation of Outer Surface Potential of Biological Membrane Using Electrophoretic Data.

    PubMed

    Silva, Izan M; Castro, Maria Clícia S; Silva, Dilson; Cortez, Célia M

    2016-06-01

    In this paper, we present the results of a study on the influence of hydrodynamic effects on the surface potentials of the erythrocyte membrane, comparing two different models formulated to simulate the electrophoretic movement of a biological cell: the classical Helmholtz-Smoluchowski model and a model presented by Hsu et al. (1996). This model considers hydrodynamic effects to describe the distribution of the fluid velocity. The electric potential equation was obtained from the non-linear Poisson-Boltzmann equation, considering the spatial distribution of electrical charges fixed in glycocalyx and cytoplasmic proteins, as well as electrolyte charges and ones fixed on the surfaces of lipidic bilayer. Our results show that the Helmholtz-Smoluchowski model is not able to reflect the real forces responsible to the electrophoretic behavior of cell, because it does not take account the hydrodynamic effects of glycocalyx. This charged network that covers cellular surface constitutes a complex physical system whose electromechanical characteristics cannot be neglected. Then, supporting the hypothesis of other authors, we suggest that, in electrophoretic motion analyses of cells, the classical model represents a limiting case of models that take into account hydrodynamic effects to describe the velocity distribution of fluid. PMID:27276378

  1. Recommended improvements to the DS02 dosimetry system's calculation of organ doses and their potential advantages for the Radiation Effects Research Foundation.

    PubMed

    Cullings, Harry M

    2012-03-01

    The Radiation Effects Research Foundation (RERF) uses a dosimetry system to calculate radiation doses received by the Japanese atomic bomb survivors based on their reported location and shielding at the time of exposure. The current system, DS02, completed in 2003, calculates detailed doses to 15 particular organs of the body from neutrons and gamma rays, using new source terms and transport calculations as well as some other improvements in the calculation of terrain and structural shielding, but continues to use methods from an older system, DS86, to account for body self-shielding. Although recent developments in models of the human body from medical imaging, along with contemporary computer speed and software, allow for improvement of the calculated organ doses, before undertaking changes to the organ dose calculations, it is important to evaluate the improvements that can be made and their potential contribution to RERF's research. The analysis provided here suggests that the most important improvements can be made by providing calculations for more organs or tissues and by providing a larger series of age- and sex-specific models of the human body from birth to adulthood, as well as fetal models. PMID:22262817

  2. Rotational excitation of symmetric top molecules by collisions with atoms: Close coupling, coupled states, and effective potential calculations for NH3-He

    NASA Technical Reports Server (NTRS)

    Green, S.

    1976-01-01

    The formalism for describing rotational excitation in collisions between symmetric top rigid rotors and spherical atoms is presented both within the accurate quantum close coupling framework and also the coupled states approximation of McGuire and Kouri and the effective potential approximation of Rabitz. Calculations are reported for thermal energy NH3-He collisions, treating NH3 as a rigid rotor and employing a uniform electron gas (Gordon-Kim) approximation for the intermolecular potential. Coupled states are found to be in nearly quantitative agreement with close coupling results while the effective potential method is found to be at least qualitatively correct. Modifications necessary to treat the inversion motion in NH3 are discussed.

  3. Calculating the potential to emit particulate matter

    SciTech Connect

    Vaart, D.R. van der

    1996-09-01

    As the implementation of the 1990 amendments to the Clean Air Act, and Title V in particular, continues, questions regarding the calculation of a facility`s potential to emit continue to surface. The US Environmental Protection Agency (EPA) has provided limited guidance decisions, although many are still being made during Title V implementation. This paper discusses what is meant by PM-10 and the validity of using sieve analysis in estimating particulate emissions. Title V of the Clean Air Amendments of 1990, and the accompanying regulations in 40 CFR Part 70, defines a major source subject to Title V by calculating its potential emissions of all regulated pollutants, both criteria and hazardous air pollutants. For PM, the threshold emission rate is 100 tons per year (tpy) for applicability to Title V. Much discussion has ensued regarding a definition of PM for the purpose of determining a facility`s potential to emit. Recently, EPA provided guidance which indicated that only PM-10 should be considered for making this determination although many states regulate larger particles through their state implementation plan (SIP) as a surrogate for PM-10.

  4. On calculating the potential vorticity flux

    SciTech Connect

    Hsu, Pei-Chun; Diamond, P. H.

    2015-03-15

    We discuss and compare different approaches to calculating the dynamics of anisotropic flow structure formation in quasi two-dimensional turbulence based on potential vorticity (PV) transport in real space. The general structure of the PV flux in the relaxation processes is deduced non-perturbatively. The transport coefficients of the PV flux are then systematically calculated using perturbation theory. We develop two non-perturbative relaxation models: the first is a mean field theory for the dynamics of minimum enstrophy relaxation based on the requirement that the mean flux of PV dissipates total potential enstrophy but conserves total fluid kinetic energy. The results show that the structure of PV flux has the form of a sum of a positive definite hyper-viscous and a negative or positive viscous flux of PV. Turbulence spreading is shown to be related to PV mixing via the link of turbulence energy flux to PV flux. In the relaxed state, the ratio of the PV gradient to zonal flow velocity is homogenized. This homogenized quantity sets a constraint on the amplitudes of PV and zonal flow in the relaxed state. The second relaxation model is derived from symmetry principles alone. The form of PV flux contains a nonlinear convective term in addition to viscous and hyper-viscous terms. For both cases, the transport coefficients are calculated using perturbation theory. For a broad turbulence spectrum, a modulational calculation of the PV flux gives both a negative viscosity and a positive hyper-viscosity. For a narrow turbulence spectrum, the result of a parametric instability analysis shows that PV transport is also convective. In both relaxation and perturbative analyses, it is shown that turbulent PV transport is sensitive to flow structure, and the transport coefficients are nonlinear functions of flow shear.

  5. Singularity embedding method in potential flow calculations

    NASA Technical Reports Server (NTRS)

    Jou, W. H.; Huynh, H.

    1982-01-01

    The so-called H-type mesh is used in a finite-element (or finite-volume) calculation of the potential flow past an airfoil. Due to coordinate singularity at the leading edge, a special singular trial function is used for the elements neighboring the leading edge. The results using the special singular elements are compared to those using the regular elements. It is found that the unreasonable pressure distribution obtained by the latter is removed by the embedding of the singular element. Suggestions to extend the present method to transonic cases are given.

  6. Effective one-body potential fitted for many-body interactions associated with a Jastrow function: application to the quantum Monte Carlo calculations

    NASA Astrophysics Data System (ADS)

    Umezawa, Naoto; Austin, Brian; Lester, William A., Jr.

    2009-03-01

    An efficient method of optimizing a Slater determinant, D, in the Jastrow-Slater-type wave function, FD, is suggested. Here, the so-called transcorrelated Hamiltonian, 1F H F, which is a similarity transformation of the usual Hamiltonian of an electronic system with respect to a Jastrow function F, is fitted to an effective Hamiltonian, Heff= ∑i^N ( -12 2̂i+ v(ri) ), in which all the electron-electron and electron-neucleus interactions are represented by a one-body potential, v(r). A single-particle Schr"odinger equation is then solved by using v(r) to determine the orbitals, of which the Slater determinant consists. The obtained orbitals improve the atomic total energies in the variational Monte Carlo calculations compared to those given by the density-functional-based orbitals. Advantages of using the optimized orbitals in the diffusion Monte Carlo calculations are also discussed.

  7. Relativistic Density Functional Theory with Optimized Effective Potential and Self-Interaction Correction : Application to Atomic Structure Calculations ( Z = 2 to 106)

    NASA Astrophysics Data System (ADS)

    Tong, X. M.; Chu, S. I.

    1998-05-01

    We introduce a self-interaction-free relativistic density functional theory (DFT) for the treatment of both the static and dynamical properties of many-electron atoms (X.M. Tong and S.I. Chu, Phys. Rev. A57), 855 (1998).. The theory is based on the extension of our recent development of non-relativistic DFT treatment (X.M. Tong and S.I. Chu, Phys. Rev. A55), 3406 (1997). with optimized effective potential (OEP) and self-interaction-corrction (SIC) to the relativistic domain. The relativistic OEP/SIC procedure yields orbital-independent single- particle local potential with proper long-range Coulombic (-1/r) behavior and is capable of providing accurate description of the ground, excited, and autoionizing states. The method is applied to the atomic structure calculations of atoms with Z = 2 to 106. Good agreement with the experimental data for both the ionization potentials (obtained from the highest occupied orbital energies) and individual orbital binding energies is obtained across the periodic table. To our knowledage, this is the first DFT calculation that has achieved such a quantitative accuracy. Detailed results will be presented.

  8. Ground-state properties of alkali dimers and their cations (including the elements Li, Na, and K) from ab initio calculations with effective core polarization potentials

    NASA Astrophysics Data System (ADS)

    Müller, Wolfgang; Meyer, Wilfried

    1984-04-01

    Extensive all-electron SCF and valence CI calculations are presented for alkali dimer systems with consideration of intershell correlation effects by use of an effective core polarization potential (CPP), which contains only a single adjustable atomic parameter. High accuracy is obtained for the ground-state spectroscopic constants of the studied molecules. The maximum deviations from accurate experimental data are as follows: 1% or 0.03 Å for Re, 2% or 100 cm-1 for De, 0.5% or 1 cm-1 for ωe, and 0.2% or 100 cm-1 for ionization energies. For experimentally uncertain or unknown values reliable predictions can thus be made. The calculated dipole moments for LiK and NaK agree with experiment to within 0.1%, but for LiNa we obtain a deviation of 8% or 0.036 D. An analysis of molecular core polarization contributions reveals the reasons for some systematic defects in previous pseudopotential calculations.

  9. Relativistically corrected nuclear magnetic resonance chemical shifts calculated with the normalized elimination of the small component using an effective potential-NMR chemical shifts of molybdenum and tungsten

    NASA Astrophysics Data System (ADS)

    Filatov, Michael; Cremer, Dieter

    2003-07-01

    A new method for relativistically corrected nuclear magnetic resonance (NMR) chemical shifts is developed by combining the individual gauge for the localized orbital approach for density functional theory with the normalized elimination of a small component using an effective potential. The new method is used for the calculation of the NMR chemical shifts of 95Mo and 183W in various molybdenum and tungsten compounds. It is shown that quasirelativistic corrections lead to an average improvement of calculated NMR chemical shift values by 300 and 120 ppm in the case of 95Mo and 183W, respectively, which is mainly due to improvements in the paramagnetic contributions. The relationship between electronic structure of a molecule and the relativistic paramagnetic corrections is discussed. Relativistic effects for the diamagnetic part of the magnetic shielding caused by a relativistic contraction of the s,p orbitals in the core region concern only the shielding values, however, have little consequence for the shift values because of the large independence from electronic structure and a cancellation of these effects in the shift values. It is shown that the relativistic corrections can be improved by level shift operators and a B3LYP hybrid functional, for which Hartree-Fock exchange is reduced to 15%.

  10. Two-component relativistic density-functional calculations of the dimers of the halogens from bromine through element 117 using effective core potential and all-electron methods.

    PubMed

    Mitin, Alexander V; van Wüllen, Christoph

    2006-02-14

    A two-component quasirelativistic Hamiltonian based on spin-dependent effective core potentials is used to calculate ionization energies and electron affinities of the heavy halogen atom bromine through the superheavy element 117 (eka-astatine) as well as spectroscopic constants of the homonuclear dimers of these atoms. We describe a two-component Hartree-Fock and density-functional program that treats spin-orbit coupling self-consistently within the orbital optimization procedure. A comparison with results from high-order Douglas-Kroll calculations--for the superheavy systems also with zeroth-order regular approximation and four-component Dirac results--demonstrates the validity of the pseudopotential approximation. The density-functional (but not the Hartree-Fock) results show very satisfactory agreement with theoretical coupled cluster as well as experimental data where available, such that the theoretical results can serve as an estimate for the hitherto unknown properties of astatine, element 117, and their dimers. PMID:16483205

  11. Extreme tunnelling in methylamine dehydrogenase revealed by hybrid QM/MM calculations: potential energy surface profile for methylamine and ethanolamine substrates and kinetic isotope effect values

    NASA Astrophysics Data System (ADS)

    Tresadern, Gary; Wang, Hong; Faulder, Paul F.; Burton, Neil A.; Hillier, Ian H.

    The rate-determining proton transfer step in the amine reduction reaction catalysed by the enzyme methylamine dehydrogenase has been studied using a hybrid quantum mechanical/molecular mechanical (QM/MM) model. Variational transition state theory, combined with multidimensional tunnelling corrections, has been employed to calculate reaction rate constants, and hence deuterium kinetic isotope effects (KIE). To render these calculations computationally feasible, the electronic structure was described using a PM3 method with specific reaction parameters obtained by a fit to energetics obtained at a high level for a small model system. Compared to the use of standard parameters, these revised parameters result in a considerable improvement in the predicted KIE values and activation energy. For both methylamine and ethanolamine substrates, through-barrier, rather than over-barrier, motion is found to dominate with KIE values that are large and close to the experimental values. A major difference between the two substrates is that, for ethanolamine, different hydrogen bonding structures involving the substrate hydroxyl are possible, leading to very different potential energy surfaces with KIE values covering a considerable range. We speculate that this is the origin of the differing temperature behaviour observed for the KIEs of the two substrates.

  12. Calculations of 21 Λ-S and 42 Ω states of BC molecule: Potential energy curves, spectroscopic parameters and spin-orbit coupling effect.

    PubMed

    Xing, Wei; Shi, Deheng; Sun, Jinfeng; Zhu, Zunlue

    2016-01-15

    The potential energy curves (PECs) were calculated for the 42 Ω states generated from the 21 Λ-S states (X(4)Σ(-), A(4)Π, B(4)Σ(-), a(2)Π, b(2)Σ(-), c(2)Δ, d(2)Σ(+), e(2)Π, 3(2)Π, 4(2)Π, 5(2)Π, 2(2)Σ(-), 3(2)Σ(-), 2(2)Σ(+), 3(2)Σ(+), 2(2)Δ, 3(2)Δ, 1(4)Σ(+), 2(4)Π, 1(4)Δ and 1(2)Φ), which originated from the lowest two dissociation channels, B((2)Pu)+C((3)Pg) and B((2)Pu)+C((1)Dg), of the BC molecule. The PECs were calculated for internuclear separations from 0.08 to 1.10 nm using the CASSCF method, which was followed by the icMRCI approach with the aug-cc-pV6Z basis set. Of these 21 Λ-S states, the e(2)Π, 2(2)Δ, 2(2)Σ(-), 4(2)Π, 1(2)Φ and 3(2)Δ possess the double wells. The A(4)Π, a(2)Π, c(2)Δ, 2(4)Π, 4(2)Π, 5(2)Π, 1(4)Δ and 1(2)Φ states are inverted with the spin-orbit coupling (SOC) effect taken into account. The first well of e(2)Π state and the second well of 4(2)Π and 2(2)Δ states do not have any vibrational states whether with or without the SOC effect included. All the Λ-S and Ω states involved in this paper are bound states. Scalar relativistic correction was included by the third-order Douglas-Kroll Hamiltonian approximation at the level of an aug-cc-pV5Z basis set. Core-valence correlation correction was included at the level of an aug-cc-pCV5Z basis set. The SOC effect was accounted for by the state interaction method with the Breit-Pauli Hamiltonian and the all-electron cc-pCV5Z basis set. The PECs of all the states were extrapolated to the complete basis set limit. The spectroscopic parameters were obtained. The vibrational properties of several Λ-S and Ω states with the relatively shallow wells were evaluated. The SOC effect on the spectroscopic parameters is not obvious for almost all the states. The spectroscopic properties reported in this paper can be expected to be reliably predicted ones. PMID:26476070

  13. An assessment of a modified potential flow code for calculating the effect of small geometric change on the pressure and forces of supercritical airfoils

    NASA Technical Reports Server (NTRS)

    Hicks, R. M.

    1982-01-01

    Wind-tunnel test data for two closely related supercritical airfoils were compared with calculations obtained from a nonconservative, potential flow code over a Mach number range from 0.20 to 0.80. The potential flow code includes an iterated, integral boundary-layer correction. The theoretical pressure distributions correlated more closely with the experimental pressure distributions when the flow was entirely subsonic or subsonic with a small supersonic zone than when the flow contained a large supersonic zone. The predicted drag level was below the experimental values at nearly all test conditions and the difference in drag level for the two airfoils was not accurately predicted.

  14. Design and development of a new micro-beam treatment planning system: effectiveness of algorithms of optimization and dose calculations and potential of micro-beam treatment.

    PubMed

    Tachibana, Hidenobu; Kojima, Hiroyuki; Yusa, Noritaka; Miyajima, Satoshi; Tsuda, Akihisa; Yamashita, Takashi

    2012-07-01

    A new treatment planning system (TPS) was designed and developed for a new treatment system, which consisted of a micro-beam-enabled linac with robotics and a real-time tracking system. We also evaluated the effectiveness of the implemented algorithms of optimization and dose calculations in the TPS for the new treatment system. In the TPS, the optimization procedure consisted of the pseudo Beam's-Eye-View method for finding the optimized beam directions and the steepest-descent method for determination of beam intensities. We used the superposition-/convolution-based (SC-based) algorithm and Monte Carlo-based (MC-based) algorithm to calculate dose distributions using CT image data sets. In the SC-based algorithm, dose density scaling was applied for the calculation of inhomogeneous corrections. The MC-based algorithm was implemented with Geant4 toolkit and a phase-based approach using a network-parallel computing. From the evaluation of the TPS, the system can optimize the direction and intensity of individual beams. The accuracy of the dose calculated by the SC-based algorithm was less than 1% on average with the calculation time of 15 s for one beam. However, the MC-based algorithm needed 72 min for one beam using the phase-based approach, even though the MC-based algorithm with the parallel computing could decrease multiple beam calculations and had 18.4 times faster calculation speed using the parallel computing. The SC-based algorithm could be practically acceptable for the dose calculation in terms of the accuracy and computation time. Additionally, we have found a dosimetric advantage of proton Bragg peak-like dose distribution in micro-beam treatment. PMID:22544809

  15. Sequence and conformation effects on ionization potential and charge distribution of homo-nucleobase stacks using M06-2X hybrid density functional theory calculations

    PubMed Central

    Rooman, Marianne; Wintjens, René

    2013-01-01

    DNA is subject to oxidative damage due to radiation or by-products of cellular metabolism, thereby creating electron holes that migrate along the DNA stacks. A systematic computational analysis of the dependence of the electronic properties of nucleobase stacks on sequence and conformation was performed here, on the basis of single- and double-stranded homo-nucleobase stacks of 1–10 bases or 1–8 base pairs in standard A-, B-, and Z-conformation. First, several levels of theory were tested for calculating the vertical ionization potentials of individual nucleobases; the M06-2X/6-31G* hybrid density functional theory method was selected by comparison with experimental data. Next, the vertical ionization potential, and the Mulliken charge and spin density distributions were calculated and considered on all nucleobase stacks. We found that (1) the ionization potential decreases with the number of bases, the lowest being reached by Gua≡Cyt tracts; (2) the association of two single strands into a double-stranded tract lowers the ionization potential significantly (3) differences in ionization potential due to sequence variation are roughly three times larger than those due to conformational modifications. The charge and spin density distributions were found (1) to be located toward the 5′-end for single-stranded Gua-stacks and toward the 3′-end for Cyt-stacks and basically delocalized over all bases for Ade- and Thy-stacks; (2) the association into double-stranded tracts empties the Cyt- and Thy-strands of most of the charge and all the spin density and concentrates them on the Gua- and Ade-strands. The possible biological implications of these results for transcription are discussed. PMID:23582046

  16. Sequence and conformation effects on ionization potential and charge distribution of homo-nucleobase stacks using M06-2X hybrid density functional theory calculations.

    PubMed

    Rooman, Marianne; Wintjens, René

    2014-04-01

    DNA is subject to oxidative damage due to radiation or by-products of cellular metabolism, thereby creating electron holes that migrate along the DNA stacks. A systematic computational analysis of the dependence of the electronic properties of nucleobase stacks on sequence and conformation was performed here, on the basis of single- and double-stranded homo-nucleobase stacks of 1-10 bases or 1-8 base pairs in standard A-, B-, and Z-conformation. First, several levels of theory were tested for calculating the vertical ionization potentials of individual nucleobases; the M06-2X/6-31G* hybrid density functional theory method was selected by comparison with experimental data. Next, the vertical ionization potential, and the Mulliken charge and spin density distributions were calculated and considered on all nucleobase stacks. We found that (1) the ionization potential decreases with the number of bases, the lowest being reached by Gua≡Cyt tracts; (2) the association of two single strands into a double-stranded tract lowers the ionization potential significantly (3) differences in ionization potential due to sequence variation are roughly three times larger than those due to conformational modifications. The charge and spin density distributions were found (1) to be located toward the 5'-end for single-stranded Gua-stacks and toward the 3'-end for Cyt-stacks and basically delocalized over all bases for Ade- and Thy-stacks; (2) the association into double-stranded tracts empties the Cyt- and Thy-strands of most of the charge and all the spin density and concentrates them on the Gua- and Ade-strands. The possible biological implications of these results for transcription are discussed. PMID:23582046

  17. Surface potentials and the calculated selectivity of ion channels.

    PubMed Central

    Miedema, Henk

    2002-01-01

    Ion channels catalyze the transport of ions across biological membranes. A proper understanding of ion-channel functioning is essential to our knowledge of cell physiology, and, in this context, ion-channel selectivity is a key concept. The extent to which a channel permeates two ion species, a and b, is expressed by the permeability ratio, P(a)/P(b). This paper addresses a complication in the calculation of P(a)/P(b) that is related to the existence of surface potentials (psi) and that so far has not been fully appreciated. This paper shows the rather surprising effect of psi on the calculated P(a)/P(b) of a channel that is permeable to two ion species of different valence. If we ignore psi, we conclude, for instance, P(a) > P(b). If we implement psi in the calculation of P(a)/P(b), we may, however, conclude exactly the reverse, i.e., P(a) < P(b). Because electrostatic potentials arise at the surface of essentially all biological membranes, this paper argues for a more critical evaluation of ion channel selectivity measurements. PMID:11751304

  18. The Calculation of the Electrostatic Potential of Infinite Charge Distributions

    ERIC Educational Resources Information Center

    Redzic, Dragan V.

    2012-01-01

    We discuss some interesting aspects in the calculation of the electrostatic potential of charge distributions extending to infinity. The presentation is suitable for the advanced undergraduate level. (Contains 3 footnotes.)

  19. Potential Dependence of Electrochemical Barriers from ab Initio Calculations.

    PubMed

    Chan, Karen; Nørskov, Jens K

    2016-05-01

    We present a simple and computationally efficient method to determine the potential dependence of the activation energies for proton-electron transfer from a single ab initio barrier calculation. We show that the potential dependence of the activation energy is given by the partial charge transferred at the transition state. The method is evaluated against the potential dependence determined explicitly through multiple calculations at varying potential. We show that the transfer coefficient is given by the charge transferred from the initial to transition state, which has significant implications for electrochemical kinetics. PMID:27088442

  20. Calculation of molecular free energies in classical potentials

    NASA Astrophysics Data System (ADS)

    Farhi, Asaf; Singh, Bipin

    2016-02-01

    Free energies of molecules can be calculated by quantum chemistry computations or by normal mode classical calculations. However, the first can be computationally impractical for large molecules and the second is based on the assumption of harmonic dynamics. We present a novel, accurate and complete calculation of molecular free energies in standard classical potentials. In this method we transform the molecule by relaxing potential terms which depend on the coordinates of a group of atoms in that molecule and calculate the free energy difference associated with the transformation. Then, since the transformed molecule can be treated as non-interacting systems, the free energy associated with these atoms is analytically or numerically calculated. This two-step calculation can be applied to calculate free energies of molecules or free energy difference between (possibly large) molecules in a general environment. We demonstrate the method in free energy calculations for methanethiol and butane molecules in vacuum and solvent. We suggest the potential application of free energy calculation of chemical reactions in classical molecular simulations.

  1. Comparison of the Effective Fragment Potential Method with Symmetry-Adapted Perturbation Theory in the Calculation of Intermolecular Energies for Ionic Liquids.

    PubMed

    Tan, Samuel Y S; Izgorodina, Ekaterina I

    2016-06-14

    The effective fragment potential (EFP) method that decomposes the interaction energy as a sum of the five fundamental forces-electrostatic, exchange-repulsion, polarization, dispersion, and charge transfer-was applied to a large test set of ionic liquid ion pairs and compared against the state-of-the-art method, Symmetry-Adapted Perturbation Theory (SAPT). The ion pairs include imidazolium and pyrrolidinium cations combined with anions that are routinely used in the field of ionic liquids. The aug-cc-pVDZ, aug-cc-pVTZ, and 6-311++G(d,p) basis sets were used for EFP, while SAPT2+3/aug-cc-pVDZ provided the benchmark energies. Differences between the two methods were found to be large, and strongly dependent on the anion type. For the aug-cc-pVTZ basis set, which produced the least errors, average relative errors were between 2.3% and 18.4% for pyrrolidinium ion pairs and between 2.1% and 27.7% for imidazolium ion pairs for each individual energetic component (excluding charge transfer), as well as the total interaction energy. Charge transfer gave the largest relative errors: 56% and 63% on average for pyrrolidinium- and imidazolium-based ion pairs, respectively. Scaling of the EFP components against SAPT2+3 showed improvement for polarization (induction) and dispersion terms, thus indicating potential for the development of cost-effective alternatives for intermolecular induction and dispersion potentials for ionic liquids. PMID:27116302

  2. Accurate ionization potential of semiconductors from efficient density functional calculations

    NASA Astrophysics Data System (ADS)

    Ye, Lin-Hui

    2016-07-01

    Despite its huge successes in total-energy-related applications, the Kohn-Sham scheme of density functional theory cannot get reliable single-particle excitation energies for solids. In particular, it has not been able to calculate the ionization potential (IP), one of the most important material parameters, for semiconductors. We illustrate that an approximate exact-exchange optimized effective potential (EXX-OEP), the Becke-Johnson exchange, can be used to largely solve this long-standing problem. For a group of 17 semiconductors, we have obtained the IPs to an accuracy similar to that of the much more sophisticated G W approximation (GWA), with the computational cost of only local-density approximation/generalized gradient approximation. The EXX-OEP, therefore, is likely as useful for solids as for finite systems. For solid surfaces, the asymptotic behavior of the vx c has effects similar to those of finite systems which, when neglected, typically cause the semiconductor IPs to be underestimated. This may partially explain why standard GWA systematically underestimates the IPs and why using the same GWA procedures has not been able to get an accurate IP and band gap at the same time.

  3. Generalized transition state theory calculations for the reactions D+H2 and H+D2 using an accurate potential energy surface: Explanation of the kinetic isotope effect

    NASA Astrophysics Data System (ADS)

    Garrett, Bruce C.; Truhlar, Donald G.

    1980-03-01

    Rate constants are calculated for the reactions D+H2→DH+H and H+D2→HD+D and compared to measured values. An accurate potential energy surface, based on the ab initio calculations of Liu and Siegbahn, was used. Rates were calculated using both conventional transition state theory and canonical variational theory. In the former, the generalized transition state dividing surface is located at the saddle point; in the latter it is located to maximize the generalized free energy of activation. We show that, in the absence of tunneling corrections, locating the generalized-transition-state dividing surface variationally has an important quantitative effect on the predicted rate constants for these systems and that, when tunneling is included, most of the effect of using a better dividing surface can be included in conventional transition state theory for these systems by using a consistent transmission coefficient for quantal scattering by the vibrationally adiabatic potential energy curve. Tunneling effects are important for these reactions even for temperatures larger than 400 K. We show how to separate classical recrossing effects from quantal corrections on reaction-coordinate motion in both the transmission coefficients and the kinetic isotope effects. Our most complete calculations are in excellent agreement with most of the measured rate constants and kinetic isotope effects.

  4. Calculation of the potentials and 3D electric fields in a proton decay detector

    SciTech Connect

    Lari, R.J.; Dawson, J.W.; Turner, L.R.

    1987-01-01

    An electrostatic detector for measuring the lifetime of the proton has been modeled in three dimensions. Linear hexahedral finite elements were used and the potential obtained at all nodes. The three components of the electric field were calculated and used to determine field lines, calculate drift fields and drift times. Effective aperture calculations agreed with the measurements.

  5. Calculating Standard Reduction Potentials of [4Fe–4S] Proteins

    SciTech Connect

    Perrin, Bradley S.; Niu, Shuqiang; Ichiye, Toshiko

    2013-03-15

    The oxidation–reduction potentials of electron transfer proteins determine the driving forces for their electron transfer reactions. Although the type of redox site determines the intrinsic energy required to add or remove an electron, the electrostatic interaction energy between the redox site and its surrounding environment can greatly shift the redox potentials. Here, a method for calculating the reduction potential versus the standard hydrogen electrode, E°, of a metalloprotein using a combinatio of density functional theory and continuum electrostatics is presented. This work focuses on the methodology for the continuum electrostatics calculations, including various factors that may affect the accuracy. The calculations are demonstrated using crystal structures of six homologous HiPIPs, which give E° that are in excellent agreement with experimental results.

  6. New approach to calculating the potential energy of colliding nuclei

    SciTech Connect

    Kurmanov, R. S.; Kosenko, G. I.

    2014-12-15

    The differential method proposed by the present authors earlier for the reduction of volume integrals in calculating the potential energy of a compound nucleus is generalized to the case of two interacting nuclei. The Coulomb interaction energy is obtained for the cases of a sharp and a diffuse boundary of nuclei, while the nuclear interaction energy is found only for nuclei with a sharp boundary, the finiteness of the nuclear-force range being taken into account. The present method of calculations permits reducing the time it takes to compute the potential energy at least by two orders of magnitude.

  7. Computer programs for calculating potential flow in propulsion system inlets

    NASA Technical Reports Server (NTRS)

    Stockman, N. O.; Button, S. L.

    1973-01-01

    In the course of designing inlets, particularly for VTOL and STOL propulsion systems, a calculational procedure utilizing three computer programs evolved. The chief program is the Douglas axisymmetric potential flow program called EOD which calculates the incompressible potential flow about arbitrary axisymmetric bodies. The other two programs, original with Lewis, are called SCIRCL AND COMBYN. Program SCIRCL generates input for EOD from various specified analytic shapes for the inlet components. Program COMBYN takes basic solutions output by EOD and combines them into solutions of interest, and applies a compressibility correction.

  8. Assessment of the accuracy of shape-consistent relativistic effective core potentials using multireference spin-orbit configuration interaction singles and doubles calculations of the ground and low-lying excited states of U(4+) and U(5+).

    PubMed

    Beck, Eric V; Brozell, Scott R; Blaudeau, Jean-Philippe; Burggraf, Larry W; Pitzer, Russell M

    2009-11-12

    Multireference spin-orbit configuration interaction calculations were used to determine the accuracy of 60-, 68-, and 78-electron shape-consistent relativistic effective core potentials (RECPs) for uranium V and VI ground and low-lying excited states. Both 5f(n) and (5f6d)(n), (n = 1, 2) reference spaces were investigated using correlation-consistent double-zeta quality basis sets. Accuracy was assessed against gas-phase experimental spectra. The 68-electron RECP calculations yielded low relative and rms errors and predicted the empirical ordering of states most consistently. PMID:19888778

  9. Fully relativistic multiple scattering calculations for general potentials

    NASA Astrophysics Data System (ADS)

    Ebert, H.; Braun, J.; Ködderitzsch, D.; Mankovsky, S.

    2016-02-01

    The formal basis for fully relativistic Korringa-Kohn-Rostoker (KKR) or multiple scattering calculations for the electronic Green function in case of a general potential is discussed. Simple criteria are given to identify situations that require to distinguish between right- and left-hand-side solutions to the Dirac equation when setting up the electronic Green function. In addition, various technical aspects of an implementation of the relativistic KKR for general local and nonlocal potentials will be discussed.

  10. Effect of soil type patterns on the variability of bare soil evaporation within a field: comparison of eddy covariance measurements with potential and actual evaporation calculations

    NASA Astrophysics Data System (ADS)

    Vanderborght, J.; Graf, A.; Steenpass, C.; Scharnagl, B.; Prolingheuer, N.; Herbst, M.; Vereecken, H.

    2009-12-01

    Bare soil evaporation was measured with the eddy-covariance method at the Selhausen field site. The site has a distinct gradient in soil texture with a considerably higher stone content at the upper part of the field. Because of this gradient, a spatial variation in evaporation fluxes in the field is expected. Because of the higher stone content at the upper part of the field, it is expected that the water that is stored in the soil surface layer and can be evaporated at a maximal evaporation rate, which is determined by the energy that is available for evaporation, is considerable smaller in the upper than in the lower part of the field. We investigated whether this hypothesis is supported by eddy covariance (EC) measurements of the evaporation fluxes at the field site. The EC measurements were combined with a footprint model that predicts the location of the soil surface that contributes to the measured evaporation flux. In this way, evaporation measurements of the two parts of the field site could be distinguished. However, since only one EC station was available, simultaneous evaporation measurements for the two field parts were not available. As a consequence, the datasets of measurements had to be interpreted and put into context of the meteorological and soil hydrological conditions. The potential evapotranspiration was calculated using the FAO method (Allen et al., 1998) to represent the meteorological conditions whereas a simple soil evaporation model (Boesten and Stroosnijder, 1986) was used to represent the influence of the precipitation and soil hydrological conditions on the actual evaporation rate. Since different soil parameters were required to describe the evaporation measurements for the upper and lower part of the plot, our starting hypothesis that more water is evaporated in the lower part of the field could be confirmed. Allen, R. G., L. S. Pereira, D. Raes, and M. Smith (1998), Crop evapotranspiration: Guidelines for computing crop water

  11. Calculation of interaction-induced spectra using complex absorbing potentials

    SciTech Connect

    Gustafsson, Magnus; Antipov, Sergey V.

    2010-10-29

    A complex absorbing potential method is implemented for calculation of collision-induced spectra. The scheme provides a way to avoid the integration of the Schroedinger equation to very large separations of the collisional pair. The method is tested by reproducing a previously computed absorption spectrum for H-He at two different temperatures.

  12. Numerical Green's functions in optical potential calculations for positron scattering from argon and neon

    NASA Technical Reports Server (NTRS)

    Bartschat, K.; Mceachran, R. P.; Stauffer, A. D.

    1990-01-01

    An optical potential method was applied to the calculation of positron scattering from the noble gases in order to determine the effect of open excitation channels on the shape of differential scattering cross sections.

  13. Off disk-center potential field calculations using vector magnetograms

    NASA Technical Reports Server (NTRS)

    Venkatakrishnan, P.; Gary, G. Allen

    1989-01-01

    A potential field calculation for off disk-center vector magnetograms that uses all the three components of the measured field is investigated. There is neither any need for interpolation of grid points between the image plane and the heliographic plane nor for an extension or a truncation to a heliographic rectangle. Hence, the method provides the maximum information content from the photospheric field as well as the most consistent potential field independent of the viewing angle. The introduction of polarimetric noise produces a less tolerant extrapolation procedure than using the line-of-sight extrapolation, but the resultant standard deviation is still small enough for the practical utility of this method.

  14. Misfit effects in adhesion calculations

    NASA Astrophysics Data System (ADS)

    Schnitker, Jurgen; Srolovitz, David J.

    1998-03-01

    The work of adhesion of bimaterial interfaces is commonly computed using quantum mechanical methods in which the two materials are strained into coherency. There is no relaxation of the coherency by the formation of an array of interfacial misfit dislocations, contrary to what is commonly observed for essentially all systems other than very thin films. In this paper, we investigate the errors introduced into the work of adhesion associated with the assumption of coherency. Series of atomistic simulations in two and three dimensions are performed using a simple Lennard-Jones-type model potential. We demonstrate that the assumption of coherency introduces errors that increase rapidly with misfit (for small misfit) and can easily be of the order of several tens of percent. We trace the source of these errors to the neglect of the elastic fields of misfit dislocations and to the variation in the number of bonds per unit interfacial area with misfit when coherency is assumed. Suggestions are made to minimize and/or correct for this error.

  15. Accurate potential energy curve of the LiH+ molecule calculated with explicitly correlated Gaussian functions.

    PubMed

    Tung, Wei-Cheng; Adamowicz, Ludwik

    2014-03-28

    Very accurate calculations of the ground-state potential energy curve (PEC) of the LiH(+) ion performed with all-electron explicitly correlated Gaussian functions with shifted centers are presented. The variational method is employed. The calculations involve optimization of nonlinear exponential parameters of the Gaussians performed with the aid of the analytical first derivatives of the energy determined with respect to the parameters. The diagonal adiabatic correction is also calculated for each PEC point. The PEC is then used to calculate the vibrational energies of the system. In that calculation, the non-adiabatic effects are accounted for by using an effective vibrational mass obtained by the minimization of the difference between the vibrational energies obtained from the calculations where the Born-Oppenheimer approximation was not assumed and the results of the present calculations. PMID:24697449

  16. Accurate potential energy curve of the LiH+ molecule calculated with explicitly correlated Gaussian functions

    NASA Astrophysics Data System (ADS)

    Tung, Wei-Cheng; Adamowicz, Ludwik

    2014-03-01

    Very accurate calculations of the ground-state potential energy curve (PEC) of the LiH+ ion performed with all-electron explicitly correlated Gaussian functions with shifted centers are presented. The variational method is employed. The calculations involve optimization of nonlinear exponential parameters of the Gaussians performed with the aid of the analytical first derivatives of the energy determined with respect to the parameters. The diagonal adiabatic correction is also calculated for each PEC point. The PEC is then used to calculate the vibrational energies of the system. In that calculation, the non-adiabatic effects are accounted for by using an effective vibrational mass obtained by the minimization of the difference between the vibrational energies obtained from the calculations where the Born-Oppenheimer approximation was not assumed and the results of the present calculations.

  17. Proliferation Potential of Accelerator-Drive Systems: Feasibility Calculations

    SciTech Connect

    Riendeau, C.D.; Moses, D.L.; Olson, A.P.

    1998-11-01

    Accelerator-driven systems for fissile materials production have been proposed and studied since the early 1950s. Recent advances in beam power levels for small accelerators have raised the possibility that such use could be feasible for a potential proliferator. The objective of this study is to review the state of technology development for accelerator-driven spallation neutron sources and subcritical reactors. Energy and power requirements were calculated for a proton accelerator-driven neutron spallation source and subcritical reactors to produce a significant amount of fissile material--plutonium.

  18. Modified method of aerodynamic resistance calculation and its application to potential evapotranspiration estimation

    NASA Astrophysics Data System (ADS)

    Rodný, Marek; Nolz, Reinhard; Novák, Viliam; Hlaváčiková, Hana; Loiskandl, Willibald; Himmelbauer, Margarita

    2016-04-01

    The aim of this study was to present and validate an alternative evapotranspiration calculation procedure that includes specific expression for the aerodynamic resistance. Calculated daily potential evapotranspiration totals were compared to the results of FAO56 procedure application and to the results of measurements taken with a precision weighing lysimeter permanently grown with irrigated, short grass. For the examination period from March 17 through October 31, 2011, it was found that daily potential evapotranspiration estimates obtained by both calculation procedures fitted well to the lysimeter measurements. Potential evapotranspiration daily totals calculated with the use of the proposed aerodynamic resistance calculation procedure gave better results for days with higher evapotranspiration, compared to the FAO56 method. The most important is that the approach based on the proposed alternative aerodynamic resistance could be effectively used even for a wide variety of crops, because it is not limited to any particular crop.

  19. Surface calculations with asymptotically long-ranged potentials in the full-potential linearized augmented plane-wave method

    NASA Astrophysics Data System (ADS)

    Ye, Lin-Hui

    2015-09-01

    Although the supercell method has been widely used for surface calculations, it only works well with short-ranged potentials, but meets difficulty when the potential decays very slowly into the vacuum. Unfortunately, the exact exchange-correlation potential of the density functional theory is asymptotically long ranged, and therefore is not easily handled by use of supercells. This paper illustrates that the authentic slab geometry, another technique for surface calculations, is not affected by this issue: It works equally well with both short- and long-ranged potentials, with the computational cost and the convergence speed being essentially the same. Using the asymptotically long-ranged Becke-Roussel'89 exchange potential as an example, we have calculated six surfaces of various types. We found that accurate potential values can be obtained even in extremely low density regions of more than 100 Å away from the surface. This high performance allows us to explore the asymptotic region, and prove with clean numerical evidence that the Becke-Roussel'89 potential satisfies the correct asymptotic behavior for slab surfaces, as it does for finite systems. Our finding further implies that the Slater component of the exact exchange optimized effective potential is responsible for the asymptotic behavior, not only for jellium slabs, but for slabs of any type. The Becke-Roussel'89 potential may therefore be used to build asymptotically correct model exchange potentials applicable to both finite systems and slab surfaces.

  20. Efficient calculation of potential energy surfaces for the generation of vibrational wave functions.

    PubMed

    Rauhut, Guntram

    2004-11-15

    An automatic procedure for the generation of potential energy surfaces based on high level ab initio calculations is described. It allows us to determine the vibrational wave functions for molecules of up to ten atoms. Speedups in computer time of about four orders of magnitude in comparison to standard implementations were achieved. Effects due to introduced approximations--within the computation of the potential--on fundamental modes obtained from vibrational self-consistent field and vibrational configuration interaction calculations are discussed. Benchmark calculations are provided for formaldehyde and 1,2,5-oxadiazole (furazan). PMID:15538851

  1. A Molecular Full-Potential LMTO Calculation for Copper Clusters

    NASA Astrophysics Data System (ADS)

    Datta, Radhika Prosad; Banerjea, Amitava; Mookerjee, Abhijit; Bhattacharyya, A. K.

    We study the electronic properties of small (10-20 atoms) copper clusters using the newly-developed molecular full-potential linearized muffin-tin orbital two-centre-fit (TCF) method of Methfessel and van Schilfgaarde. The geometric structures of the clusters had earlier been determined by us through simulated annealing using the Equivalent Crystal Theory to compute total energies. We report the variation of the binding energy, as obtained from the TCF calculations, with cluster size and compare these to the binding energies determined, for the same structures, from the ECT. We also show the variation of the HOMO-LUMO gap with cluster size, and the pseudo-density of states for select cluster sizes.

  2. Electronic structure calculations toward new potentially AChE inhibitors

    NASA Astrophysics Data System (ADS)

    de Paula, A. A. N.; Martins, J. B. L.; Gargano, R.; dos Santos, M. L.; Romeiro, L. A. S.

    2007-10-01

    The main purpose of this study was the use of natural non-isoprenoid phenolic lipid of cashew nut shell liquid from Anacardium occidentale as lead material for generating new potentially candidates of acetylcholinesterase inhibitors. Therefore, we studied the electronic structure of 15 molecules derivatives from the cardanol using the following groups: methyl, acetyl, N, N-dimethylcarbamoyl, N, N-dimethylamine, N, N-diethylamine, piperidine, pyrrolidine, and N-benzylamine. The calculations were performed at RHF level using 6-31G, 6-31G(d), 6-31+G(d) and 6-311G(d,p) basis functions. Among the proposed compounds we found that the structures with substitution by acetyl, N, N-dimethylcarbamoyl, N, N-dimethylamine, and pyrrolidine groups were better correlated to rivastigmine indicating possible activity.

  3. Multi-grid calculation of transonic potential flows

    NASA Technical Reports Server (NTRS)

    Caughey, D. A.; Shmilovich, A.

    1985-01-01

    The finite-volume method discussed by Jameson and Caughey (1977), and Caughey and Jameson (1979, 1980) has made it possible to calculate the transonic potential flow past any configuration for which a suitable boundary-conforming coordinate grid can be constructed. However, computations for practical three-dimensional problems have remained quite expensive in terms of the required computer time. The reason for this is primarily related to the large number of grid cells necessary for adequate resolution in these complex three-dimensional problems, taking into account the large number of iterations required to achieve even modest convergence on these fine grids. The present chapter provides a description of work directed at removing this latter difficulty by making use of the multigrid method. Attention is given to finite-volume formulation, multigrid iteration, geometrical aspects, and computed results.

  4. A submerged singularity method for calculating potential flow velocities at arbitrary near-field points

    NASA Technical Reports Server (NTRS)

    Maskew, B.

    1976-01-01

    A discrete singularity method has been developed for calculating the potential flow around two-dimensional airfoils. The objective was to calculate velocities at any arbitrary point in the flow field, including points that approach the airfoil surface. That objective was achieved and is demonstrated here on a Joukowski airfoil. The method used combined vortices and sources ''submerged'' a small distance below the airfoil surface and incorporated a near-field subvortex technique developed earlier. When a velocity calculation point approached the airfoil surface, the number of discrete singularities effectively increased (but only locally) to keep the point just outside the error region of the submerged singularity discretization. The method could be extended to three dimensions, and should improve nonlinear methods, which calculate interference effects between multiple wings, and which include the effects of force-free trailing vortex sheets. The capability demonstrated here would extend the scope of such calculations to allow the close approach of wings and vortex sheets (or vortices).

  5. Efficient calculation of potential energy surfaces for the generation of vibrational wave functions

    NASA Astrophysics Data System (ADS)

    Rauhut, Guntram

    2004-11-01

    An automatic procedure for the generation of potential energy surfaces based on high level ab initio calculations is described. It allows us to determine the vibrational wave functions for molecules of up to ten atoms. Speedups in computer time of about four orders of magnitude in comparison to standard implementations were achieved. Effects due to introduced approximations—within the computation of the potential—on fundamental modes obtained from vibrational self-consistent field and vibrational configuration interaction calculations are discussed. Benchmark calculations are provided for formaldehyde and 1,2,5-oxadiazole (furazan).

  6. Helium-surface interaction potential of Sb(111) from scattering experiments and close-coupling calculations

    NASA Astrophysics Data System (ADS)

    Mayrhofer-Reinhartshuber, M.; Kraus, P.; Tamtögl, A.; Miret-Artés, S.; Ernst, W. E.

    2013-11-01

    Helium atom scattering (HAS) was used to study the antimony Sb(111) surface beyond the hard-wall model. HAS angular distributions and drift spectra show a number of selective adsorption resonance features, which correspond to five bound-state energies for He atoms trapped in the surface-averaged He-Sb(111) potential. As their best representation, a 9-3 potential with a depth of 4.4±0.1 meV was determined. Furthermore, the charge density corrugation of the surface was analyzed using close-coupling calculations. By using a hybrid potential, consisting of a corrugated Morse potential (short range) and a 9-3 potential (long range), a peak-to-peak corrugation of 17% was obtained. The kinematic focusing effects that occurred were in good agreement with surface phonon dispersion curves from already published density functional perturbation theory calculations.

  7. Electron-N/sub 2/ scattering calculations with a parameter-free model polarization potential

    SciTech Connect

    Morrison, M.A.; Saha, B.C.; Gibson, T.L.

    1987-10-15

    We have extended our variationally determined nonadiabatic polarization potential (Gibson and Morrison, Phys. Rev. A 29, 2497 (1984)) to the e-N/sub 2/ system and calculated elastic, total momentum transfer, and rotational excitation cross sections. This model potential, which requires no scaling and contains no adjustable parameters, is presented in tabular and analytic (fitted) form for possible use in future studies. We evaluated the static potential at the near-Hartree-Fock level of accuracy and included exchange effects exactly via the linear algebraic method of Collins and Schneider (Phys. Rev. A 24, 2387 (1981)). Diverse cross sections based on this model are in excellent agreement with existing experiment. We also compare various scattering quantities calculated with our model to prior theoretical results and to newly determined numbers using two other model potentials: a cutoff phenomenological form and the correlation-polarization potential of O'Connell and Lane (Phys. Rev. A 27, 1893 (1983)).

  8. Bonn potential and shell-model calculations for N=126 isotones

    SciTech Connect

    Coraggio, L.; Covello, A.; Gargano, A.; Itaco, N.; Kuo, T. T. S.

    1999-12-01

    We have performed shell-model calculations for the N=126 isotones {sup 210}Po, {sup 211}At, and {sup 212}Rn using a realistic effective interaction derived from the Bonn-A nucleon-nucleon potential by means of a G-matrix folded-diagram method. The calculated binding energies, energy spectra, and electromagnetic properties show remarkably good agreement with the experimental data. The results of this paper complement those of our previous study on neutron hole Pb isotopes, confirming that realistic effective interactions are now able to reproduce with quantitative accuracy the spectroscopic properties of complex nuclei. (c) 1999 The American Physical Society.

  9. NATIONAL EMISSION STANDARDS FOR HAZARDOUS AIR POLLUTANTS (NESHAP) SUBPART H RADIONUCLIDES POTENTIAL TO EMIT CALCULATIONS

    SciTech Connect

    EARLEY JN

    2008-07-23

    This document provides an update of the status of stacks on the Hanford Site and the potential radionuclide emissions, i.e., emissions that could occur with no control devices in place. This review shows the calculations that determined whether the total effective dose equivalent (TEDE) received by the maximum public receptor as a result of potential emissions from any one of these stacks would exceed 0.1 millirem/year. Such stacks require continuous monitoring of the effluent, or other monitoring, to meet the requirements of Washington Administrative code (WAC) 246-247-035(1)(a)(ii) and WAC 246-247-075(1), -(2), and -(6). This revised update reviews the potential-to-emit (PTE) calculations of 31 stacks for Fluor Hanford, Inc. Of those 31 stacks, 11 have the potential to cause a TEDE greater than 0.1 mrem/year.

  10. Effective Potential in Noncommutative BTZ Black Hole

    NASA Astrophysics Data System (ADS)

    Sadeghi, Jafar; Shajiee, Vahid Reza

    2016-02-01

    In this paper, we investigated the noncommutative rotating BTZ black hole and showed that such a space-time is not maximally symmetric. We calculated effective potential for the massive and the massless test particle by geodesic equations, also we showed effect of non-commutativity on the minimum mass of BTZ black hole.

  11. Analytical second derivatives for effective core potentials

    NASA Astrophysics Data System (ADS)

    Breidung, Jürgen; Thiel, Walter; Komornicki, Andrew

    1988-12-01

    Analytical first and second derivatives for effective core potentials are reported. The computational implementation of the derivative formulas makes use of new integral routines which take advantage of the shell concept. Test calculations for H 3SnBr and F 3AsS demonstrate the efficiency of the analytical determination of harmonic force fields using effective core potentials. The spectroscopic constants of the unknown molecule F 3AsS are predicted.

  12. Absolute Binding Free Energy Calculations Using Molecular Dynamics Simulations with Restraining Potentials

    PubMed Central

    Wang, Jiyao; Deng, Yuqing; Roux, Benoît

    2006-01-01

    The absolute (standard) binding free energy of eight FK506-related ligands to FKBP12 is calculated using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. A number of features are implemented to improve the accuracy and enhance the convergence of the calculations. First, the absolute binding free energy is decomposed into sequential steps during which the ligand-surrounding interactions as well as various biasing potentials restraining the translation, orientation, and conformation of the ligand are turned “on” and “off.” Second, sampling of the ligand conformation is enforced by a restraining potential based on the root mean-square deviation relative to the bound state conformation. The effect of all the restraining potentials is rigorously unbiased, and it is shown explicitly that the final results are independent of all artificial restraints. Third, the repulsive and dispersive free energy contribution arising from the Lennard-Jones interactions of the ligand with its surrounding (protein and solvent) is calculated using the Weeks-Chandler-Andersen separation. This separation also improves convergence of the FEP/MD calculations. Fourth, to decrease the computational cost, only a small number of atoms in the vicinity of the binding site are simulated explicitly, while all the influence of the remaining atoms is incorporated implicitly using the generalized solvent boundary potential (GSBP) method. With GSBP, the size of the simulated FKBP12/ligand systems is significantly reduced, from ∼25,000 to 2500. The computations are very efficient and the statistical error is small (∼1 kcal/mol). The calculated binding free energies are generally in good agreement with available experimental data and previous calculations (within ∼2 kcal/mol). The present results indicate that a strategy based on FEP/MD simulations of a reduced GSBP atomic model sampled with conformational, translational, and orientational restraining

  13. Perturbation method to calculate the interaction potentials and electronic excitation spectra of atoms in He nanodroplets.

    PubMed

    Callegari, Carlo; Ancilotto, Francesco

    2011-06-30

    A method is proposed for the calculation of potential energy curves and related electronic excitation spectra of dopant atoms captured in/on He nanodroplets and is applied to alkali metal atoms. The method requires knowledge of the droplet density distribution at equilibrium (here calculated within a bosonic-He density functional approach) and of a set of valence electron orbitals of the bare dopant atom (here calculated by numeric solution of the Schrödinger equation in a suitably parametrized model potential). The electron-helium interaction is added as a perturbation, and potential energy curves are obtained by numeric diagonalization of the resulting Hamiltonian as a function of an effective coordinate z(A) (here the distance between the dopant atom and center of mass of the droplet, resulting in a pseudodiatomic potential). Excitation spectra are calculated for Na in the companion paper as the Franck-Condon factors between the v = 0 vibrational state in the ground electronic state and excited states of the pseudodiatomic molecule. They agree well with available experimental data, even for highly excited states where a more traditional approach fails. PMID:21434657

  14. Model potential calculation of the thermal donor energy spectrum in silicon

    NASA Astrophysics Data System (ADS)

    Chen, C. S.; Schroder, D. K.

    1988-06-01

    The two-parameter model potential originally proposed by Ning and Sah [Phys. Rev. B 4, 3468 (1971)] for calculating the ground-state energies of group V and group VI impurities in silicon is extended to the variational calculation of the thermal donor ionization energies. In the multivalley effective mass approximation, the theoretical results are in excellent agreement with the reported experimental data. This provides additional evidence for the assumption that thermal donors consist of five to thirteen oxygen atoms, as first proposed by Ourmazd, Schröter, and Bourret [J. Appl. Phys. 56, 1670 (1984)].

  15. Model potential calculation of the thermal donor energy spectrum in silicon

    SciTech Connect

    Chen, C.S.; Schroder, D.K.

    1988-06-15

    The two-parameter model potential originally proposed by Ning and Sah (Phys. Rev. B 4, 3468 (1971)) for calculating the ground-state energies of group V and group VI impurities in silicon is extended to the variational calculation of the thermal donor ionization energies. In the multivalley effective mass approximation, the theoretical results are in excellent agreement with the reported experimental data. This provides additional evidence for the assumption that thermal donors consist of five to thirteen oxygen atoms, as first proposed by Ourmazd, Schroeter, and Bourret (J. Appl. Phys. 56, 1670 (1984)).

  16. The use of model potentials in molecular calculations. II

    NASA Astrophysics Data System (ADS)

    Sakai, Y.; Huzinaga, S.

    1982-03-01

    The model potential method is applied to CO, HCl, P2, Cl2, SH2, Cu2, Br2, Ni(CO)4, and Pd(CO)4. The results are generally very satisfactory. Reduction of computing cost is substantial for molecules containing heavy atoms.

  17. Improved finite difference schemes for transonic potential calculations

    NASA Technical Reports Server (NTRS)

    Hafez, M.; Osher, S.; Whitlow, W., Jr.

    1984-01-01

    Engquist and Osher (1980) have introduced a finite difference scheme for solving the transonic small disturbance equation, taking into account cases in which only compression shocks are admitted. Osher et al. (1983) studied a class of schemes for the full potential equation. It is proved that these schemes satisfy a new discrete 'entropy inequality' which rules out expansion shocks. However, the conducted analysis is restricted to steady two-dimensional flows. The present investigation is concerned with the adoption of a heuristic approach. The full potential equation in conservation form is solved with the aid of a modified artificial density method, based on flux biasing. It is shown that, with the current scheme, expansion shocks are not possible.

  18. Comparison of Electron Elastic-Scattering Cross Sections Calculated from Two Commonly Used Atomic Potentials

    NASA Astrophysics Data System (ADS)

    Jablonski, A.; Salvat, F.; Powell, C. J.

    2004-06-01

    We have analyzed differential cross sections (DCSs) for the elastic scattering of electrons by neutral atoms that have been derived from two commonly used atomic potentials: the Thomas-Fermi-Dirac (TFD) potential and the Dirac-Hartree-Fock (DHF) potential. DCSs from the latter potential are believed to be more accurate. We compared DCSs for six atoms (H, Al, Ni, Ag, Au, and Cm) at four energies (100, 500, 1000, and 10 000 eV) from two databases issued by the National Institute of Standards and Technology in which DCSs had been obtained from the TFD and DHF potentials. While the DCSs from the two potentials had similar shapes and magnitudes, there can be pronounced deviations (up to 70%) for small scattering angles for Al, Ag, Au, and Cm. In addition, there were differences of up to 400% at scattering angles for which there were deep minima in the DCSs; at other angles, the differences were typically less than 20%. The DCS differences decreased with increasing electron energy. DCSs calculated from the two potentials were compared with measured DCSs for six atoms (He, Ne, Ar, Kr, Xe, and Hg) at energies between 50 eV and 3 keV. For Ar, the atom for which experimental data are available over the largest energy range there is good agreement between the measured DCSs and those calculated from the TFD and DHF potentials at 2 and 3 keV, but the experimental DCSs agree better with the DCSs from the DHF potential at lower energies. A similar trend is found for the other atoms. At energies less than about 1 keV, there are increasing differences between the measured DCSs and the DCSs calculated from the DHF potential. These differences were attributed to the neglect of absorption and polarizability effects in the calculations. We compare transport cross sections for H, Al, Ni, Ag, Au, and Cm obtained from the DCSs for each potential. For energies between 200 eV and 1 keV, the largest differences are about 20% (for H, Au, and Cm); at higher energies, the differences are

  19. Potential theoretic methods for far field sound radiation calculations

    NASA Technical Reports Server (NTRS)

    Hariharan, S. I.; Stenger, Edward J.; Scott, J. R.

    1995-01-01

    In the area of computational acoustics, procedures which accurately predict the far-field sound radiation are much sought after. A systematic development of such procedures are found in a sequence of papers by Atassi. The method presented here is an alternate approach to predicting far field sound based on simple layer potential theoretic methods. The main advantages of this method are: it requires only a simple free space Green's function, it can accommodate arbitrary shapes of Kirchoff surfaces, and is readily extendable to three-dimensional problems. Moreover, the procedure presented here, though tested for unsteady lifting airfoil problems, can easily be adapted to other areas of interest, such as jet noise radiation problems. Results are presented for lifting airfoil problems and comparisons are made with the results reported by Atassi. Direct comparisons are also made for the flat plate case.

  20. On the calculation of potential flow about a body in an unbounded fluid

    NASA Astrophysics Data System (ADS)

    Noblesse, F.; Triantafyllou, G.

    1980-09-01

    This study is concerned with the problem of calculating potential flow about a nonlifting body in an unbounded fluid. Several simple explicit approximations for the velocity potential are obtained and investigated numerically. Results of calculations are presented for the simple cases of potential flows due to translations of ellipsoids and ogives.

  1. Piezoelectric constants for ZnO calculated using classical polarizable core-shell potentials

    NASA Astrophysics Data System (ADS)

    Dai, Shuangxing; Dunn, Martin L.; Park, Harold S.

    2010-11-01

    We demonstrate the feasibility of using classical atomistic simulations, i.e. molecular dynamics and molecular statics, to study the piezoelectric properties of ZnO using core-shell interatomic potentials. We accomplish this by reporting the piezoelectric constants for ZnO as calculated using two different classical interatomic core-shell potentials: that originally proposed by Binks and Grimes (1994 Solid State Commun. 89 921-4), and that proposed by Nyberg et al (1996 J. Phys. Chem. 100 9054-63). We demonstrate that the classical core-shell potentials are able to qualitatively reproduce the piezoelectric constants as compared to benchmark ab initio calculations. We further demonstrate that while the presence of the shell is required to capture the electron polarization effects that control the clamped ion part of the piezoelectric constant, the major shortcoming of the classical potentials is a significant underprediction of the clamped ion term as compared to previous ab initio results. However, the present results suggest that overall, these classical core-shell potentials are sufficiently accurate to be utilized for large scale atomistic simulations of the piezoelectric response of ZnO nanostructures.

  2. Calculations supporting evaluation of potential environmental standards for Yucca Mountain

    SciTech Connect

    Duguid, J.O.; Andrews, R.W.; Brandstetter, E.; Dale, T.F.; Reeves, M.

    1994-04-01

    The Energy Policy Act of 1992, Section 801 (US Congress, 1992) provides for the US Environmental Protection Agency (EPA) to contract the National Academy of Sciences (NAS) to conduct a study and provide findings and recommendations on reasonable standards for the disposal of high-level wastes at the Yucca Mountain site. The NAS study is to provide findings and recommendations which include, among other things, whether a health-based standard based on dose to individual members of the public from releases to the accessible environment will provide a reasonable standard for the protection of the health and safety of the public. The EPA, based upon and consistent with the findings and recommendations of the NAS, is required to promulgate standards for protection of the public from releases from radioactive materials stored or disposed of in a repository at the Yucca Mountain site. This document presents a number of different ``simple`` analyses of undisturbed repository performance that are intended to provide input to those responsible for setting appropriate environmental standards for a potential repository at the Yucca Mountain site in Nevada. Each of the processes included in the analyses has been simplified to capture the primary significance of that process in containing or isolating the waste from the biosphere. In these simplified analyses, the complex waste package interactions were approximated by a simple waste package ``failure`` distribution which is defined by the initiation and rate of waste package ``failures``. Similarly, releases from the waste package and the engineered barrier system are controlled by the very near field environment and the presence and rate of advective and diffusive release processes. Release was approximated by either a simple alteration-controlled release for the high solubility radionuclides and either a diffusive or advective-controlled release for the solubility-limited radionuclides.

  3. Accurate potential energy curve of the LiH{sup +} molecule calculated with explicitly correlated Gaussian functions

    SciTech Connect

    Tung, Wei-Cheng; Adamowicz, Ludwik

    2014-03-28

    Very accurate calculations of the ground-state potential energy curve (PEC) of the LiH{sup +} ion performed with all-electron explicitly correlated Gaussian functions with shifted centers are presented. The variational method is employed. The calculations involve optimization of nonlinear exponential parameters of the Gaussians performed with the aid of the analytical first derivatives of the energy determined with respect to the parameters. The diagonal adiabatic correction is also calculated for each PEC point. The PEC is then used to calculate the vibrational energies of the system. In that calculation, the non-adiabatic effects are accounted for by using an effective vibrational mass obtained by the minimization of the difference between the vibrational energies obtained from the calculations where the Born-Oppenheimer approximation was not assumed and the results of the present calculations.

  4. On the consistent definition of spin-orbit effects calculated by relativistic effective core potentials with one-electron spin-orbit operators: Comparison of spin-orbit effects for Tl, TlH, TlH3, PbH2, and PbH4

    NASA Astrophysics Data System (ADS)

    Han, Young-Kyu; Bae, Cheolbeom; Lee, Yoon Sup

    1999-05-01

    The spin-orbit effects for Tl, TlH, TlH3, PbH2, and PbH4 are evaluated by two-component calculations using several relativistic effective core potentials (RECP) with one-electron spin-orbit operators. The used RECPs are shape-consistent RECPs derived by Wildman et al. [J. Chem. Phys. 107, 9975 (1997)] and three sets of energy-consistent (or adjusted) RECPs published by Schwerdtfeger et al. [Phys. Scr. 36, 453 (1987); J. Chem. Phys. 90, 762 (1989)], Küchle et al. [Mol. Phys. 74, 1245 (1991)], and Leininger et al. [Chem. Phys. 217, 19 (1997)]. The shape-consistent RECP results are in very good agreement with the Küchle et al. energy-consistent RECP results for all the molecules studied here and all-electron results for TlH. The RECPs of Schwerdtfeger et al. and Leininger et al. seem to provide qualitatively different spin-orbit effects. If one defines spin-free RECP as the potential average of the corresponding two-component RECP, all RECPs give very similar spin-orbit effects for all the cases. Most of the discrepancies of molecular spin-orbit effects among various RECPs reported in the literature may originate from different definitions of RECPs with or without a spin-orbit term and not from the inherent difference in spin-orbit operators.

  5. Effect of orientation anisotropy on calculating effective electrical conductivities

    NASA Astrophysics Data System (ADS)

    Myles, Timothy D.; Peracchio, Aldo A.; Chiu, Wilson K. S.

    2014-05-01

    This paper develops an analytical effective medium theory (EMT) equation for calculating the effective conductivity of a mixture based on Maxwell's and Maxwell-Garnett's theories, extended to higher volume fractions using Bruggeman's unsymmetrical treatment (BUT), with a long term goal of extending the treatment to mixtures more representative of real materials in order to calculate their effective electrical conductivity. The development accounts for spheroid shaped inclusions of varying degrees of anisotropic orientation. The orientation is described by the introduction of a distribution function. Two methodologies valid for the inclusion dilute limit were used to evaluate the effective conductivity: one based on Maxwell's far field approach, and the other based on the Maxwell-Garnett in the matrix approach. It was found that while the dilute limit equations for the effective conductivity were different, the final EMT equations derived by applying BUT collapsed to the same formula which was generalized for anisotropic orientation based on the distribution function presented.

  6. Ab initio calculation of ionization potential and electron affinity in solid-state organic semiconductors

    NASA Astrophysics Data System (ADS)

    Kang, Youngho; Jeon, Sang Ho; Cho, Youngmi; Han, Seungwu

    2016-01-01

    We investigate the vertical ionization potential (IP) and electron affinity (EA) of organic semiconductors in the solid state that govern the optoelectrical property of organic devices using a fully ab initio way. The present method combines the density functional theory and many-body perturbation theory based on G W approximations. To demonstrate the accuracy of this approach, we carry out calculations on several prototypical organic molecules. Since IP and EA depend on the molecular orientation at the surface, the molecular geometry of the surface is explicitly considered through the slab model. The computed IP and EA are in reasonable and consistent agreements with spectroscopic data on organic surfaces with various molecular arrangements. However, the transport gaps are slightly underestimated in calculations, which can be explained by different screening effects between surface and bulk regions.

  7. Potential effects of gallium on cladding materials

    SciTech Connect

    Wilson, D.F.; Beahm, E.C.; Besmann, T.M.; DeVan, J.H.; DiStefano, J.R.; Gat, U.; Greene, S.R.; Rittenhouse, P.L.; Worley, B.A.

    1997-10-01

    This paper identifies and examines issues concerning the incorporation of gallium in weapons derived plutonium in light water reactor (LWR) MOX fuels. Particular attention is given to the more likely effects of the gallium on the behavior of the cladding material. The chemistry of weapons grade (WG) MOX, including possible consequences of gallium within plutonium agglomerates, was assessed. Based on the calculated oxidation potentials of MOX fuel, the effect that gallium may have on reactions involving fission products and possible impact on cladding performance were postulated. Gallium transport mechanisms are discussed. With an understanding of oxidation potentials and assumptions of mechanisms for gallium transport, possible effects of gallium on corrosion of cladding were evaluated. Potential and unresolved issues and suggested research and development (R and D) required to provide missing information are presented.

  8. Calculation of Water Drop Trajectories to and About Arbitrary Three-Dimensional Bodies in Potential Airflow

    NASA Technical Reports Server (NTRS)

    Norment, H. G.

    1980-01-01

    Calculations can be performed for any atmospheric conditions and for all water drop sizes, from the smallest cloud droplet to large raindrops. Any subsonic, external, non-lifting flow can be accommodated; flow into, but not through, inlets also can be simulated. Experimental water drop drag relations are used in the water drop equations of motion and effects of gravity settling are included. Seven codes are described: (1) a code used to debug and plot body surface description data; (2) a code that processes the body surface data to yield the potential flow field; (3) a code that computes flow velocities at arrays of points in space; (4) a code that computes water drop trajectories from an array of points in space; (5) a code that computes water drop trajectories and fluxes to arbitrary target points; (6) a code that computes water drop trajectories tangent to the body; and (7) a code that produces stereo pair plots which include both the body and trajectories. Code descriptions include operating instructions, card inputs and printouts for example problems, and listing of the FORTRAN codes. Accuracy of the calculations is discussed, and trajectory calculation results are compared with prior calculations and with experimental data.

  9. Interface effects on calculated defect levels for oxide defects

    NASA Astrophysics Data System (ADS)

    Edwards, Arthur; Barnaby, Hugh; Schultz, Peter; Pineda, Andrew

    2014-03-01

    Density functional theory (DFT) has had impressive recent success predicting defect levels in insulators and semiconductors [Schultz and von Lillienfeld, 2009]. Such success requires care in accounting for long-range electrostatic effects. Recently, Komsa and Pasquarello have started to address this problem in systems with interfaces. We report a multiscale technique for calculating electrostatic energies for charged defects in oxide of the metal-oxide-silicon (MOS) system, but where account is taken of substrate doping density, oxide thickness, and gate bias. We use device modeling to calculate electric fields for a point charge a fixed distance from the interface, and used the field to numerically calculate the long-range electrostatic interactions. We find, for example, that defect levels in the oxide do depend on both the magnitude and the polarity the substrate doping density. Furthermore, below 20 Å, oxide thickness also has significant effects. So, transferring results directly from bulk calculations leads to inaccuracies up to 0.5 eV- half of the silicon band gap. We will present trends in defect levels as a function of device parameters. We show that these results explain previous experimental results, and we comment on their potential impact on models for NBTI. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the United States Department of Energy's National Nuclear Security Administration under co.

  10. The ozone depletion potentials on halocarbons: Their dependence of calculation assumptions

    NASA Technical Reports Server (NTRS)

    Karol, Igor L.; Kiselev, Andrey A.

    1994-01-01

    The concept of Ozone Depletion Potential (ODP) is widely used in the evaluation of numerous halocarbons and of their replacement effects on ozone, but the methods, assumptions and conditions used in ODP calculations have not been analyzed adequately. In this paper a model study of effects on ozone of the instantaneous releases of various amounts of CH3CCl3 and of CHF2Cl (HCFC-22) for several compositions of the background atmosphere are presented, aimed at understanding connections of ODP values with the assumptions used in their calculations. To facilitate the ODP computation in numerous versions for the long time periods after their releases, the above rather short-lived gases and the one-dimensional radiative photochemical model of the global annually averaged atmospheric layer up to 50 km height are used. The variation of released gas global mass from 1 Mt to 1 Gt leads to ODP value increase with its stabilization close to the upper bound of this range in the contemporary atmosphere. The same variations are analyzed for conditions of the CFC-free atmosphere of 1960's and for the anthropogenically loaded atmosphere in the 21st century according to the known IPCC 'business as usual' scenario. Recommendations for proper ways of ODP calculations are proposed for practically important cases.

  11. Calculators in the Mathematics Curriculum: Effects and Changes.

    ERIC Educational Resources Information Center

    Rabe, Rebecca Moore

    The purpose of this paper was to determine the effects of calculators in mathematics classes and to assess proposed curriculum revisions related to calculators. Twenty-six calculator studies and other selected sources were reviewed and annotated. Major conclusions of the study include: (1) calculator use has produced significant gains in…

  12. Interference effects in potential wells

    NASA Astrophysics Data System (ADS)

    Mullin, W. J.; Laloë, F.

    2015-05-01

    We propose using an array of potential wells as an interferometer in which the beam splitters are provided by tunneling during an appropriate time through the barrier between wells. This arrangement allows demonstration of generalized Hong-Ou-Mandel effects with multiple particles traversing one or several beam splitters. Other interferometer effects can occur, including a violation of the Bell-Clauser-Horne-Shimony-Holt form of the Bell inequality. With interactions, one sees various effects, including so-called fermionization, collective tunneling, and self-trapping.

  13. Quantum-Mechanical Calculation of Ionization-Potential Lowering in Dense Plasmas

    NASA Astrophysics Data System (ADS)

    Son, Sang-Kil; Thiele, Robert; Jurek, Zoltan; Ziaja, Beata; Santra, Robin

    2014-07-01

    The charged environment within a dense plasma leads to the phenomenon of ionization-potential depression (IPD) for ions embedded in the plasma. Accurate predictions of the IPD effect are of crucial importance for modeling atomic processes occurring within dense plasmas. Several theoretical models have been developed to describe the IPD effect, with frequently discrepant predictions. Only recently, first experiments on IPD in Al plasma have been performed with an x-ray free-electron laser, where their results were found to be in disagreement with the widely used IPD model by Stewart and Pyatt. Another experiment on Al, at the Orion laser, showed disagreement with the model by Ecker and Kröll. This controversy shows a strong need for a rigorous and consistent theoretical approach to calculate the IPD effect. Here, we propose such an approach: a two-step Hartree-Fock-Slater model. With this parameter-free model, we can accurately and efficiently describe the experimental Al data and validate the accuracy of standard IPD models. Our model can be a useful tool for calculating atomic properties within dense plasmas with wide-ranging applications to studies on warm dense matter, shock experiments, planetary science, inertial confinement fusion, and nonequilibrium plasmas created with x-ray free-electron lasers.

  14. Interpolation effects in tabulated interatomic potentials

    NASA Astrophysics Data System (ADS)

    Wen, M.; Whalen, S. M.; Elliott, R. S.; Tadmor, E. B.

    2015-10-01

    Empirical interatomic potentials are widely used in atomistic simulations due to their ability to compute the total energy and interatomic forces quickly relative to more accurate quantum calculations. The functional forms in these potentials are sometimes stored in a tabulated format, as a collection of data points (argument-value pairs), and a suitable interpolation (often spline-based) is used to obtain the function value at an arbitrary point. We explore the effect of these interpolations on the potential predictions by calculating the quasi-harmonic thermal expansion and finite-temperature elastic constant of a one-dimensional chain compared with molecular dynamics simulations. Our results show that some predictions are affected by the choice of interpolation regardless of the number of tabulated data points. Our results clearly indicate that the interpolation must be considered part of the potential definition, especially for lattice dynamics properties that depend on higher-order derivatives of the potential. This is facilitated by the Knowledgebase of Interatomic Models (KIM) project, in which both the tabulated data (‘parameterized model’) and the code that interpolates them to compute energy and forces (‘model driver’) are stored and given unique citeable identifiers. We have developed cubic and quintic spline model drivers for pair functional type models (EAM, FS, EMT) and uploaded them to the OpenKIM repository (https://openkim.org).

  15. Finite volume calculation of three-dimensional potential flow around a propeller

    NASA Technical Reports Server (NTRS)

    Jou, W.-H.

    1982-01-01

    The finite volume scheme of Jameson (1977) is used to calculate potential flow around a propeller rotating at high speed. An H-type mesh is generated and used successfully in the calculations. A test calculation with a thick blade cross section shows that the present code is capable of computing the propeller flow at the advance Mach number 0.8. The possible physical mechanisms which may play an important role in the propeller aerodynamics are discussed.

  16. Phase-integral calculation of phase shifts for a heavy-ion optical potential

    SciTech Connect

    Linnaeus, S.

    1986-10-01

    Phase shifts for an optical potential representing the nuclear scattering of /sup 18/O by /sup 58/Ni at 60 MeV laboratory energy are calculated by means of an arbitrary-order phase-integral formula taking account of two turning points. The results are found to be in excellent agreement with previously published numerical calculations.

  17. 40 CFR Appendix D to Part 72 - Calculation of Potential Electric Output Capacity

    Code of Federal Regulations, 2011 CFR

    2011-07-01

    ... 40 Protection of Environment 16 2011-07-01 2011-07-01 false Calculation of Potential Electric Output Capacity D Appendix D to Part 72 Protection of Environment ENVIRONMENTAL PROTECTION AGENCY (CONTINUED) AIR PROGRAMS (CONTINUED) PERMITS REGULATION Pt. 72, App. D Appendix D to Part 72—Calculation...

  18. 40 CFR Appendix D to Part 72 - Calculation of Potential Electric Output Capacity

    Code of Federal Regulations, 2010 CFR

    2010-07-01

    ... 40 Protection of Environment 16 2010-07-01 2010-07-01 false Calculation of Potential Electric Output Capacity D Appendix D to Part 72 Protection of Environment ENVIRONMENTAL PROTECTION AGENCY (CONTINUED) AIR PROGRAMS (CONTINUED) PERMITS REGULATION Pt. 72, App. D Appendix D to Part 72—Calculation...

  19. 40 CFR Appendix D to Part 72 - Calculation of Potential Electric Output Capacity

    Code of Federal Regulations, 2013 CFR

    2013-07-01

    ... 40 Protection of Environment 17 2013-07-01 2013-07-01 false Calculation of Potential Electric Output Capacity D Appendix D to Part 72 Protection of Environment ENVIRONMENTAL PROTECTION AGENCY (CONTINUED) AIR PROGRAMS (CONTINUED) PERMITS REGULATION Pt. 72, App. D Appendix D to Part 72—Calculation...

  20. Crystal-field calculations for transition-metal ions by application of an opposing potential

    DOE PAGESBeta

    Zhou, Fei; Aberg, Daniel

    2016-02-16

    We propose a fully ab initio method, the opposing crystal potential (OCP), to calculate the crystal-field parameters of transition-metal impurities in insulator hosts. Through constrained density functional calculations, OCP obtains the constraining Lagrange multipliers, which act as a cancellation potential against the crystal field and lead to spherical d-electron distribution. Furthermore, the method is applied to several insulators doped with Mn4+ and Mn2+ ions and shown to be in good agreement with experiment.

  1. Calculating Second-Order Effects in MOSFET's

    NASA Technical Reports Server (NTRS)

    Benumof, Reuben; Zoutendyk, John A.; Coss, James R.

    1990-01-01

    Collection of mathematical models includes second-order effects in n-channel, enhancement-mode, metal-oxide-semiconductor field-effect transistors (MOSFET's). When dimensions of circuit elements relatively large, effects neglected safely. However, as very-large-scale integration of microelectronic circuits leads to MOSFET's shorter or narrower than 2 micrometer, effects become significant in design and operation. Such computer programs as widely-used "Simulation Program With Integrated Circuit Emphasis, Version 2" (SPICE 2) include many of these effects. In second-order models of n-channel, enhancement-mode MOSFET, first-order gate-depletion region diminished by triangular-cross-section deletions on end and augmented by circular-wedge-cross-section bulges on sides.

  2. Improved calculation of Si sputter yield via first principles derived interatomic potential

    NASA Astrophysics Data System (ADS)

    Hossain, M. Z.; Freund, J. B.; Johnson, H. T.

    2009-04-01

    Silicon sputter yield under medium energy Ar+ ion bombardment is calculated via molecular dynamics, using a highly accurate interatomic potential for Ar-Si interactions derived from first-principles calculations. Unlike the widely used universal repulsive potentials such as the Moliere or ZBL parameterizations, this new potential, referred to as DFT-ArSi, is developed via localized basis density functional theory. Sputter yields for Si obtained with the DFT-ArSi potential at 500 eV and 1 keV incident energies are found to be within 6% and 2% of experimental results, respectively, while errors using existing potentials are typically on the order of 11%. The DFT-ArSi potential differs from existing empirical potentials in the ˜1 Å interatomic separation range which is shown to be the most important range for modeling low-to-medium energy ion bombardment.

  3. Effect of EMP fields on cell membrane potentials

    SciTech Connect

    Gailey, P.C.; Easterly, C.E.

    1993-06-01

    A simple model is presented for cell membrane potentials induced during exposure to electromagnetic pulse (EMP). Using calculated values of internal electric field strength induced during EMP exposure, the model predicts that cell membrane potentials of about 100 mV may be induced for time frames on the order of 10 ns. Possible biological effects of these potentials including electroporation area discussed.

  4. Large magnetocrystalline anisotropy in bilayer transition metal phases from first-principles full-potential calculations

    NASA Astrophysics Data System (ADS)

    Ravindran, P.; Kjekshus, A.; Fjellvåg, H.; James, P.; Nordström, L.; Johansson, B.; Eriksson, O.

    2001-04-01

    The computational framework of this study is based on the local-spin-density approximation with first-principles full-potential linear muffin-tin orbital calculations including orbital polarization (OP) correction. We have studied the magnetic anisotropy for a series of bilayer CuAu(I)-type materials such as FeX, MnX (X=Ni,Pd,Pt), CoPt, NiPt, MnHg, and MnRh in a ferromagnetic state using experimental structural parameters to understand the microscopic origin of magnetic-anisotropy energy (MAE) in magnetic multilayers. Except for MnRh and MnHg, all these phases show perpendicular magnetization. We have analyzed our results in terms of angular momentum-, spin- and site-projected density of states, magnetic-angular-momentum-projected density of states, orbital-moment density of states, and total density of states. The orbital-moment number of states and the orbital-moment anisotropy for FeX (X=Ni,Pd,Pt) are calculated as a function of band filling to study its effect on MAE. The total and site-projected spin and orbital moments for all these systems are calculated with and without OP when the magnetization is along or perpendicular to the plane. The results are compared with available experimental as well as theoretical results. Our calculations show that OP always enhances the orbital moment in these phases and brings them closer to experimental values. The changes in MAE are analyzed in terms of exchange splitting, spin-orbit splitting, and tetragonal distortion/crystal-field splitting. The calculated MAE is found to be in good agreement with experimental values when the OP correction is included. Some of the materials considered here show large magnetic anisotropy of the order of meV. In particular we found that MnPt will have a very large MAE if it could be stabilized in a ferromagnetic configuration. Our analysis indicates that apart from large spin-orbit interaction and exchange interaction from at least one of the constituents, a large crystal-field splitting

  5. Transonic flow analysis for rotors. Part 2: Three-dimensional, unsteady, full-potential calculation

    NASA Technical Reports Server (NTRS)

    Chang, I. C.

    1985-01-01

    A numerical method is presented for calculating the three-dimensional unsteady, transonic flow past a helicopter rotor blade of arbitrary geometry. The method solves the full-potential equations in a blade-fixed frame of reference by a time-marching implicit scheme. At the far-field, a set of first-order radiation conditions is imposed, thus minimizing the reflection of outgoing wavelets from computational boundaries. Computed results are presented to highlight radial flow effects in three dimensions, to compare surface pressure distributions to quasi-steady predictions, and to predict the flow field on a swept-tip blade. The results agree well with experimental data for both straight- and swept-tip blade geometries.

  6. FINITE EXPANSION METHOD FOR THE CALCULATION AND INTERPRETATION OF MOLECULAR ELECTROSTATIC POTENTIALS

    EPA Science Inventory

    Because it is useful to have the molecular electrostatic potential as an element in a complex scheme to assess the toxicity of large molecules, efficient and reliable methods are needed for the calculation and characterization of these potentials. A multicenter multipole expansio...

  7. Comparison of spin-orbit configuration interaction methods employing relativistic effective core potentials for the calculation of zero-field splittings of heavy atoms with a 2Po ground state

    NASA Astrophysics Data System (ADS)

    Buenker, Robert J.; Alekseyev, Aleksey B.; Liebermann, Heinz-Peter; Lingott, Rainer; Hirsch, Gerhard

    1998-03-01

    Computational strategies for the treatment of relativistic effects including spin-orbit coupling at a highly correlated level are compared for a number of heavy atoms: indium, iodine, thallium, and astatine. Initial tests with perturbation theory emphasize the importance of high-energy singly excited configurations which possess large spin-orbit matrix elements with the ground state. A contracted basis consisting of L-S CI eigenfunctions (LSC-SO-CI) is found to give an accurate representation of both spin-perturbed 2Po components as long as key np→pi* singly excited configurations are included. Comparison is made with a more extensive treatment in which all selected configurations of various L-S symmetries form the basis for the multireference-spin-orbit-configuration interaction (MR-SO-CI). Good agreement is obtained with experimental SO splittings for the In, I, and At atoms at a variety of levels of treatment, indicating that the L-S contracted SO-CI approach can be implemented quite effectively with relativistic effective core potentials (RECPs) for both very electronegative atoms and also for lighter electropositive elements up through the fifth row of the periodic table. The thallium atom SO splitting is more difficult to obtain accurately because of greater differences between its valence p1/2 and p3/2 spinors than in the other cases studied, but good results are also possible with the contracted SO-CI approach in this instance, provided proper care is given to the inclusion of key singly excited L-S states. The relationship between all-electron two-component SO-CI treatments and those employing RECPs is also analyzed, and it is concluded that triply excited configurations relative to the 2Po ground state are far less important than previously reported.

  8. PSOLV: a code for calculating the potentials and densities in MFTF-B

    SciTech Connect

    Colborn, J.A.

    1983-08-17

    Code PSOLV solves for potential and densities at the cardinal points of MFTF-B. The code is equipped to handle both the throttle-coil and the axicell geometries. For the throttle-coil case, the potential at point MXO is input, while the potentials and densities at points MAI, b, and A are calculated. For the axicell case, the code must additionally solve for the potentials and densities at points X and MXO. PSOLV is intended primarily for use as a subroutine in TREQ, a code being developed by Rensink that calculates the densities and potentials at the cardinal points of MFTF-B as a function of time. TREQ is to be used for modeling start-up behavior.

  9. 17 CFR 230.459 - Calculation of effective date.

    Code of Federal Regulations, 2010 CFR

    2010-04-01

    ... 17 Commodity and Securities Exchanges 2 2010-04-01 2010-04-01 false Calculation of effective date. 230.459 Section 230.459 Commodity and Securities Exchanges SECURITIES AND EXCHANGE COMMISSION GENERAL RULES AND REGULATIONS, SECURITIES ACT OF 1933 Filings; Fees; Effective Date § 230.459 Calculation...

  10. Effect of grid system on finite element calculation

    NASA Technical Reports Server (NTRS)

    Lee, K. D.; Yen, S. M.

    1980-01-01

    Detailed parametric studies of the effect of grid system on finite element calculation for potential flows were made. These studies led to the formulation of a design criteria for optimum mesh system and the development of two methods to generate the optimum mesh system. The guidelines for optimum mesh system are: (1) the mesh structure should be regular; (2) the element should be as regular and equilateral as possible; (3) the distribution of size of element should be consistent with that of flow variables to insure maximum uniformity in error distribution; (4) for non-Dirichlet boundary conditions, smaller boundary elements or higher order interpolation functions should be used; and (5) the mesh should accommodate the boundary geometry as accurately as possible. The results of the parametric studies are presented.

  11. Calculation of the polarization potential for e-N2 collisions

    NASA Technical Reports Server (NTRS)

    Onda, K.; Temkin, A.

    1983-01-01

    A polarization potential V(pol) for e-N2 collisions is calculated by the generalization of the static part of the method of polarized orbitals to molecular targets. Partial differential equations (PDE) are derived for polarized orbitals, which are functions of the distance (r) from the molecular center and angle (theta) from the molecular axis. The equations are solved with the use of the noniterative PDE technique. From the polarized orbitals a polarization potential can be constructed whose r and theta dependence is found to be significantly different from the well-known phenomenological one. The two potentials are further compared by carrying out limited hybrid-theory scattering calculations. Only those scattering results based on the calculated V(pol) are in satisfactory accord with experiment.

  12. Potential energy curves of Li+2 from all-electron EA-EOM-CCSD calculations

    NASA Astrophysics Data System (ADS)

    Musiał, Monika; Medrek, Magdalena; Kucharski, Stanisław A.

    2015-10-01

    The electron attachment (EA) equation-of-motion coupled-cluster theory provides description of the states obtained by the attachment of an electron to the reference system. If the reference is assumed to be a doubly ionised cation, then the EA results relate to the singly ionised ion. In the current work, the above scheme is applied to the calculations of the potential energy curves (PECs) of the Li+2 cation adopting the doubly ionised Li2 +2 structure as the reference system. The advantage of such computational strategy relies on the fact that the closed-shell Li2 +2 reference dissociates into closed-shell fragments (Li2 +2 ⇒ Li+ + Li+), hence the RHF (restricted Hartree-Fock) function can be used as the reference in the whole range of interatomic distances. This scheme offers the first principle method without any model or effective potential parameters for the description of the bond-breaking processes. In this study, the PECs and selected spectroscopic constants for 18 electronic states of the Li+2 ion were computed and compared with experimental and other theoretical results. †In honour of Professor Sourav Pal on the occasion of an anniversary in his private and scientific life.

  13. Exact calculation of a microscopic nucleon spin-orbit potential: Reexamination of Brieva-Rook localization

    NASA Astrophysics Data System (ADS)

    Haider, W.; Rafi, Syed; Rook, J. R.; Gambhir, Y. K.

    2016-05-01

    Brieva-Rook approximations for calculating the equivalent local microscopic nucleon-nucleus spin-orbit potential has been examined in the energy region of 65 calculating both the direct and the exchange parts of the spin-orbit potential. It has been shown that the series expansion used earlier overestimates both the real and the imaginary parts by about 50 % . Our results show that the exchange contribution to the total spin-orbit potential at low energies is only 15 % (much smaller than earlier estimates) and becomes negligible around 200-MeV nucleon incident energy. Predictions of the new spin-orbit potential with minor scaling has been found to be in fair agreement with p -208Pb scattering observables from 65 to 500 MeV.

  14. Comparisons of measured and calculated potential magnetic fields. [in solar corona

    NASA Technical Reports Server (NTRS)

    Hagyard, M. J.; Teuber, D.

    1978-01-01

    Photospheric line-of-sight and transverse-magnetic-field data obtained, with a vector magnetograph system for an isolated sunspot are described. A study of the linear polarization patterns and of the calculated transverse field lines indicates that the magnetic field of the region is very nearly potential. The H-alpha fibril structures of this region as seen in high-resolution photographs corroborate this conclusion. Consequently, a potential-field calculation is described using the measured line-of-sight fields together with assumed Neumann boundary conditions; both are necessary and sufficient for a unique solution. The computed transverse fields are then compared with the measured transverse fields to verify the potential-field model and assumed boundary values. The implications of these comparisons for the validity of magnetic-field extrapolations using potential theory are discussed.

  15. A finite volume method for calculating transonic potential flow around wings from the pressure minimum integral

    NASA Technical Reports Server (NTRS)

    Eberle, A.

    1978-01-01

    Analysis of the pressure minimum integral in the calculation of three-dimensional potential flow around wings makes it possible to use non-rectangular mesh networks for distributing the three-dimensional potential into discrete points. The method is comparatively easily expanded to the treatment of realistic airplane configurations. Shock-pressure affected pressure distributions on any wings are determined with accuracy using this method.

  16. Self-consistent calculation of hyperfine fields and adiabatic potential of impurities in iron

    NASA Astrophysics Data System (ADS)

    Kanamori, Junjiro; Akai, Hisazumi; Akai, Masako

    1984-01-01

    Hyperfine fields of impurities of the atomic number Z=1 56 at the substitutional site and those of light impurities of Z=1 9 at the interstitial sites in ferromagnetic iron are calculated by the KKR method adapted to the system containing a single impurity atom. The potential of the impurity atom is determined self-consistently by use of the local spin density functional formalism. The results for nonmagnetic sp valence impurities agree with those of the previous nonself-consistent calculation by Katayama-Yoshida, Terakura and Kanamori except for a few cases, confirming their theory of the systematic variation of hyperfine fields. The calculation for magnetic impurities of transition elements is presented for the first time in this paper. The calculations mentioned so far assume that impurities are situated at the center of each site. For the purpose of discussing the stability of the impurity positions, the change of the adiabatic potential due to displacements from the center is calculated by carrying out similar self-consistent calculations for off-center impurity positions. It is concluded that positive muon and some light impurities including boron will be displaced from the center when trapped in a vacancy.

  17. Calculating solar photovoltaic potential on residential rooftops in Kailua Kona, Hawaii

    NASA Astrophysics Data System (ADS)

    Carl, Caroline

    As carbon based fossil fuels become increasingly scarce, renewable energy sources are coming to the forefront of policy discussions around the globe. As a result, the State of Hawaii has implemented aggressive goals to achieve energy independence by 2030. Renewable electricity generation using solar photovoltaic technologies plays an important role in these efforts. This study utilizes geographic information systems (GIS) and Light Detection and Ranging (LiDAR) data with statistical analysis to identify how much solar photovoltaic potential exists for residential rooftops in the town of Kailua Kona on Hawaii Island. This study helps to quantify the magnitude of possible solar photovoltaic (PV) potential for Solar World SW260 monocrystalline panels on residential rooftops within the study area. Three main areas were addressed in the execution of this research: (1) modeling solar radiation, (2) estimating available rooftop area, and (3) calculating PV potential from incoming solar radiation. High resolution LiDAR data and Esri's solar modeling tools and were utilized to calculate incoming solar radiation on a sample set of digitized rooftops. Photovoltaic potential for the sample set was then calculated with the equations developed by Suri et al. (2005). Sample set rooftops were analyzed using a statistical model to identify the correlation between rooftop area and lot size. Least squares multiple linear regression analysis was performed to identify the influence of slope, elevation, rooftop area, and lot size on the modeled PV potential values. The equations built from these statistical analyses of the sample set were applied to the entire study region to calculate total rooftop area and PV potential. The total study area statistical analysis findings estimate photovoltaic electric energy generation potential for rooftops is approximately 190,000,000 kWh annually. This is approximately 17 percent of the total electricity the utility provided to the entire island in

  18. A web-based calculator for estimating the profit potential of grain segregation by protein concentration

    Technology Transfer Automated Retrieval System (TEKTRAN)

    By ignoring spatial variability in grain quality, conventional harvesting systems may increase the likelihood that growers will not capture price premiums for high quality grain found within fields. The Grain Segregation Profit Calculator was developed to demonstrate the profit potential of segregat...

  19. Calculation of equation of state of QCD at zero temperature and finite chemical potential

    NASA Astrophysics Data System (ADS)

    Jiang, Yu; Li, Ning; Sun, Wei-Min; Zong, Hong-Shi

    2010-09-01

    In this paper we calculate the equation of state (EOS) of QCD at zero temperature and finite chemical potential by using several models of quark propagators including the Dyson-Schwinger equations (DSEs) model, the hard-dense-loop (HDL) approximation and the quasi-particle model. The results are analyzed and compared with the known results in the literature.

  20. Computer programs for calculating two-dimensional potential flow in and about propulsion system inlets

    NASA Technical Reports Server (NTRS)

    Hawk, J. D.; Stockman, N. O.; Farrell, C. A., Jr.

    1978-01-01

    Incompressible potential flow calculations are presented that were corrected for compressibility in two-dimensional inlets at arbitrary operating conditions. Included are a statement of the problem to be solved, a description of each of the computer programs, and sufficient documentation, including a test case, to enable a user to run the program.

  1. Improved computer programs for calculating potential flow in propulsion system inlets

    NASA Technical Reports Server (NTRS)

    Stockman, N. O.; Farrell, C. A., Jr.

    1977-01-01

    Computer programs to calculate the incompressible potential flow corrected for compressibility in axisymmetric inlets at arbitrary operating conditions are presented. Included are a statement of the problem to be solved, a description of each of the programs and sufficient documentation, including a test case, to enable a user to run the programs.

  2. Computer programs for calculating two-dimensional potential flow through deflected nozzles

    NASA Technical Reports Server (NTRS)

    Hawk, J. D.; Stockman, N. O.

    1979-01-01

    Computer programs to calculate the incompressible potential flow, corrected for compressibility, in two-dimensional nozzles at arbitrary operating conditions are presented. A statement of the problem to be solved, a description of each of the computer programs, and sufficient documentation, including a test case, to enable a user to run the program are included.

  3. The Inversion Potential of Ammonia: An Intrinsic Reaction Coordinate Calculation for Student Investigation

    ERIC Educational Resources Information Center

    Halpern, Arthur M.; Ramachandran, B. R.; Glendening, Eric D.

    2007-01-01

    A report is presented to describe how students can be empowered to construct the full, double minimum inversion potential for ammonia by performing intrinsic reaction coordinate calculations. This work can be associated with the third year physical chemistry lecture laboratory or an upper level course in computational chemistry.

  4. A consistent calculation of the chemical potential for dense simple fluids.

    PubMed

    Bomont, Jean-Marc

    2006-05-28

    A general method to calculate the excess chemical potential betamuex, that is based on the Kirkwood coupling parameter's dependence of the correlation functions, is presented. The expression for the one particle bridge function B(1)r is derived for simple fluids with spherical interactions. Only the knowledge of the bridge function B(2)r is required. The accuracy of our approach is illustrated for a dense hard sphere fluid. As far as B(2)r is considered as exact, B(1)r is found to be, at high densities, the normalized bridge function -B(2)rB(2)(r=0). This expression ensures a consistent calculation of the excess chemical potential by satisfying implicitly the Gibbs-Duhem constraint. Only the pressure-consistency condition is necessary to calculate the structural and thermodynamic properties of the fluid. PMID:16774388

  5. Carbon dioxide hydrate phase equilibrium and cage occupancy calculations using ab initio intermolecular potentials.

    PubMed

    Velaga, Srinath C; Anderson, Brian J

    2014-01-16

    Gas hydrate deposits are receiving increased attention as potential locations for CO2 sequestration, with CO2 replacing the methane that is recovered as an energy source. In this scenario, it is very important to correctly characterize the cage occupancies of CO2 to correctly assess the sequestration potential as well as the methane recoverability. In order to predict accurate cage occupancies, the guest–host interaction potential must be represented properly. Earlier, these potential parameters were obtained by fitting to experimental equilibrium data and these fitted parameters do not match with those obtained by second virial coefficient or gas viscosity data. Ab initio quantum mechanical calculations provide an independent means to directly obtain accurate intermolecular potentials. A potential energy surface (PES) between H2O and CO2 was computed at the MP2/aug-cc-pVTZ level and corrected for basis set superposition error (BSSE), an error caused due to the lower basis set, by using the half counterpoise method. Intermolecular potentials were obtained by fitting Exponential-6 and Lennard-Jones 6-12 models to the ab initio PES, correcting for many-body interactions. We denoted this model as the “VAS” model. Reference parameters for structure I carbon dioxide hydrate were calculated using the VAS model (site–site ab initio intermolecular potentials) as Δμ(w)(0) = 1206 ± 2 J/mol and ΔH(w)(0) = 1260 ± 12 J/mol. With these reference parameters and the VAS model, pure CO2 hydrate equilibrium pressure was predicted with an average absolute deviation of less than 3.2% from the experimental data. Predictions of the small cage occupancy ranged from 32 to 51%, and the large cage is more than 98% occupied. The intermolecular potentials were also tested by calculating the pure CO2 density and diffusion of CO2 in water using molecular dynamics simulations. PMID:24328234

  6. The Effect of Calculator Use on College Students' Mathematical Performance

    ERIC Educational Resources Information Center

    Boyle, Robert W.; Farreras, Ingrid G.

    2015-01-01

    This experiment tested the effect that calculator use had on 200 randomly assigned college students' mathematical performance. The purposes of the current experiment were twofold: to measure the level of mathematical preparation of current college students, and to test whether calculators improve mathematical performance in such students as it…

  7. A study of potential numerical pitfalls in GPU-based Monte Carlo dose calculation

    NASA Astrophysics Data System (ADS)

    Magnoux, Vincent; Ozell, Benoît; Bonenfant, Éric; Després, Philippe

    2015-07-01

    The purpose of this study was to evaluate the impact of numerical errors caused by the floating point representation of real numbers in a GPU-based Monte Carlo code used for dose calculation in radiation oncology, and to identify situations where this type of error arises. The program used as a benchmark was bGPUMCD. Three tests were performed on the code, which was divided into three functional components: energy accumulation, particle tracking and physical interactions. First, the impact of single-precision calculations was assessed for each functional component. Second, a GPU-specific compilation option that reduces execution time as well as precision was examined. Third, a specific function used for tracking and potentially more sensitive to precision errors was tested by comparing it to a very high-precision implementation. Numerical errors were found in two components of the program. Because of the energy accumulation process, a few voxels surrounding a radiation source end up with a lower computed dose than they should. The tracking system contained a series of operations that abnormally amplify rounding errors in some situations. This resulted in some rare instances (less than 0.1%) of computed distances that are exceedingly far from what they should have been. Most errors detected had no significant effects on the result of a simulation due to its random nature, either because they cancel each other out or because they only affect a small fraction of particles. The results of this work can be extended to other types of GPU-based programs and be used as guidelines to avoid numerical errors on the GPU computing platform.

  8. Calculation of transient potential rise on the wind turbine struck by lightning.

    PubMed

    Xiaoqing, Zhang

    2014-01-01

    A circuit model is proposed in this paper for calculating the transient potential rise on the wind turbine struck by lightning. The model integrates the blade, sliding contact site, and tower and grounding system of the wind turbine into an equivalent circuit. The lightning current path from the attachment point to the ground can be fully described by the equivalent circuit. The transient potential responses are obtained in the different positions on the wind turbine by solving the circuit equations. In order to check the validity of the model, the laboratory measurement is made with a reduced-scale wind turbine. The measured potential waveform is compared with the calculated one and a better agreement is shown between them. The practical applicability of the model is also examined by a numerical example of a 2 MW Chinese-built wind turbine. PMID:25254231

  9. Calculation of Transient Potential Rise on the Wind Turbine Struck by Lightning

    PubMed Central

    Xiaoqing, Zhang

    2014-01-01

    A circuit model is proposed in this paper for calculating the transient potential rise on the wind turbine struck by lightning. The model integrates the blade, sliding contact site, and tower and grounding system of the wind turbine into an equivalent circuit. The lightning current path from the attachment point to the ground can be fully described by the equivalent circuit. The transient potential responses are obtained in the different positions on the wind turbine by solving the circuit equations. In order to check the validity of the model, the laboratory measurement is made with a reduced-scale wind turbine. The measured potential waveform is compared with the calculated one and a better agreement is shown between them. The practical applicability of the model is also examined by a numerical example of a 2 MW Chinese-built wind turbine. PMID:25254231

  10. Full dimension Rb2He ground triplet potential energy surface and quantum scattering calculations

    NASA Astrophysics Data System (ADS)

    Guillon, Grégoire; Viel, Alexandra; Launay, Jean-Michel

    2012-05-01

    We have developed a three-dimensional potential energy surface for the lowest triplet state of the Rb2He complex. A global analytic fit is provided as in the supplementary material [see supplementary material at http://dx.doi.org/10.1063/1.4709433E-JCPSA6-136-034218 for the corresponding Fortran code]. This surface is used to perform quantum scattering calculations of 4He and 3He colliding with 87Rb2 in the partial wave J = 0 at low and ultralow energies. For the heavier helium isotope, the computed vibrational relaxation probabilities show a broad and strong shape resonance for a collisional energy of 0.15 K and a narrow Feshbach resonance at about 17 K for all initial Rb2 vibrational states studied. The broad resonance corresponds to an efficient relaxation mechanism that does not occur when 3He is the colliding partner. The Feshbach resonance observed at higher collisional energy is robust with respect to the isotopic substitution. However, its effect on the vibrational relaxation mechanism is faint for both isotopes.

  11. An exploratory study of a finite difference method for calculating unsteady transonic potential flow

    NASA Technical Reports Server (NTRS)

    Bennett, R. M.; Bland, S. R.

    1979-01-01

    A method for calculating transonic flow over steady and oscillating airfoils was developed by Isogai. The full potential equation is solved with a semi-implicit, time-marching, finite difference technique. Steady flow solutions are obtained from time asymptotic solutions for a steady airfoil. Corresponding oscillatory solutions are obtained by initiating an oscillation and marching in time for several cycles until a converged periodic solution is achieved. The method is described in general terms and results for the case of an airfoil with an oscillating flap are presented for Mach numbers 0.500 and 0.875. Although satisfactory results are obtained for some reduced frequencies, it is found that the numerical technique generates spurious oscillations in the indicial response functions and in the variation of the aerodynamic coefficients with reduced frequency. These oscillations are examined with a dynamic data reduction method to evaluate their effects and trends with reduced frequency and Mach number. Further development of the numerical method is needed to eliminate these oscillations.

  12. Full dimension Rb2He ground triplet potential energy surface and quantum scattering calculations.

    PubMed

    Guillon, Grégoire; Viel, Alexandra; Launay, Jean-Michel

    2012-05-01

    We have developed a three-dimensional potential energy surface for the lowest triplet state of the Rb(2)He complex. A global analytic fit is provided as in the supplementary material [see supplementary material at http://dx.doi.org/10.1063/1.4709433 for the corresponding Fortran code]. This surface is used to perform quantum scattering calculations of (4)He and (3)He colliding with (87)Rb(2) in the partial wave J = 0 at low and ultralow energies. For the heavier helium isotope, the computed vibrational relaxation probabilities show a broad and strong shape resonance for a collisional energy of 0.15 K and a narrow Feshbach resonance at about 17 K for all initial Rb(2) vibrational states studied. The broad resonance corresponds to an efficient relaxation mechanism that does not occur when (3)He is the colliding partner. The Feshbach resonance observed at higher collisional energy is robust with respect to the isotopic substitution. However, its effect on the vibrational relaxation mechanism is faint for both isotopes. PMID:22583230

  13. Calculation procedures for potential and viscous flow solutions for engine inlets

    NASA Technical Reports Server (NTRS)

    Albers, J. A.; Stockman, N. O.

    1973-01-01

    The method and basic elements of computer solutions for both potential flow and viscous flow calculations for engine inlets are described. The procedure is applicable to subsonic conventional (CTOL), short-haul (STOL), and vertical takeoff (VTOL) aircraft engine nacelles operating in a compressible viscous flow. The calculated results compare well with measured surface pressure distributions for a number of model inlets. The paper discusses the uses of the program in both the design and analysis of engine inlets, with several examples given for VTOL lift fans, acoustic splitters, and for STOL engine nacelles. Several test support applications are also given.

  14. Exploring the MRCI method for calculating interaction energies: application to the HeNe potential curve

    NASA Astrophysics Data System (ADS)

    van de Bovenkamp, J.; van Mourik, T.; van Duijneveldt, F. B.

    A multi-reference configuration interaction (MRCI) method is described, which is devised for the calculation of interaction energies of van der Waals complexes and applied to calculating the HeNe potential energy curve. The MRCI calculations make use of a generalized Poplecorrection in order to account for the lack of size consistency. The orbital space is partitioned into three subspaces: the first active space (AS1), which contains the strongly occupied orbitals; the second active space (AS2), which contains the main intra-correlating orbitals; and the external space (ES). It is shown that, to keep the error below 0.2K in the excitation scheme and the active orbital space it is sufficient to include only sigma-orbitals in AS2 and to use an excitation scheme (labelled Qq-MRCI) that encompasses only up to quadruply excited configurations. The final active orbital space (AS2) turned out to be 2s(He), 2psigma(He), 3s(Ne), 3psigma(Ne) and 3dsigma(Ne). Other MRCI variants, in which most or all quadruply excited configurations were deleted from the CI expansion (Qt- and Tt-MRCI), were found to be inadequate. Using the Qq-MRCI scheme together with a 197-orbital 'interaction optimized' basis set (IO197), the MRCI interaction energy at R = 5.7 a0 was calculated to be-21.12K. The corresponding values at the MP4 and CCSD(T) levels of theory are-20.06K and-20.99K, respectively, indicating that the MP4 method is inappropriate for highly accurate calculations on this system. Fitting the calculated data using a generalized Morse function, including an additional C6/R6 term to account for a correct long-range behaviour of the potential, the MRCI well depth was calculated to be-21.16K at Req = 5.73 a0. The MRCI and CCSD(T) potentials have the same quality and are found to be in good agreement with the HartreeFock dispersion (HFD-B) potential of Keil, M., Danielson, L. J., and Dunlop, P. J., 1991, J. Chem. Phys., 94, 296. It is concluded that, for basis IO197, the CCSD(T) method is

  15. A refined quartic potential energy surface and large scale vibrational calculations for S0 thiophosgene.

    PubMed

    Rashev, Svetoslav; Moule, David C

    2015-04-01

    In this work we present a full 6D quartic potential energy surface (PES) for S0 thiophosgene in curvilinear symmetrized bond-angle coordinates. The PES was refined starting from an ab initio field derived from acc-pVTZ basis set with CCSD(T) corrections for electron correlation. In the present calculations we used our variational method that was recently tested on formaldehyde and some of its isotopomers, along with additional improvements. The lower experimentally known vibrational levels for 35Cl2CS were reproduced quite well in the calculations, which can be regarded as a test for the feasibility of the obtained quartic PES. PMID:25615683

  16. A general method for constructing multidimensional molecular potential energy surfaces from {ital ab} {ital initio} calculations

    SciTech Connect

    Ho, T.; Rabitz, H.

    1996-02-01

    A general interpolation method for constructing smooth molecular potential energy surfaces (PES{close_quote}s) from {ital ab} {ital initio} data are proposed within the framework of the reproducing kernel Hilbert space and the inverse problem theory. The general expression for an {ital a} {ital posteriori} error bound of the constructed PES is derived. It is shown that the method yields globally smooth potential energy surfaces that are continuous and possess derivatives up to second order or higher. Moreover, the method is amenable to correct symmetry properties and asymptotic behavior of the molecular system. Finally, the method is generic and can be easily extended from low dimensional problems involving two and three atoms to high dimensional problems involving four or more atoms. Basic properties of the method are illustrated by the construction of a one-dimensional potential energy curve of the He{endash}He van der Waals dimer using the exact quantum Monte Carlo calculations of Anderson {ital et} {ital al}. [J. Chem. Phys. {bold 99}, 345 (1993)], a two-dimensional potential energy surface of the HeCO van der Waals molecule using recent {ital ab} {ital initio} calculations by Tao {ital et} {ital al}. [J. Chem. Phys. {bold 101}, 8680 (1994)], and a three-dimensional potential energy surface of the H{sup +}{sub 3} molecular ion using highly accurate {ital ab} {ital initio} calculations of R{umlt o}hse {ital et} {ital al}. [J. Chem. Phys. {bold 101}, 2231 (1994)]. In the first two cases the constructed potentials clearly exhibit the correct asymptotic forms, while in the last case the constructed potential energy surface is in excellent agreement with that constructed by R{umlt o}hse {ital et} {ital al}. using a low order polynomial fitting procedure. {copyright} {ital 1996 American Institute of Physics.}

  17. Comparison of molecular dynamics methods and interatomic potentials for calculating the thermal conductivity of silicon

    NASA Astrophysics Data System (ADS)

    Howell, P. C.

    2012-12-01

    We compare the molecular dynamics Green-Kubo and direct methods for calculating thermal conductivity κ, using as a test case crystalline silicon at temperatures T in the range 500-1000 K (classical regime). We pay careful attention to the convergence with respect to simulation size and duration and to the procedures used to fit the simulation data. We show that in the Green-Kubo method the heat current autocorrelation function is characterized by three decay processes, of which the slowest lasts several tens of picoseconds so that convergence requires several tens of nanoseconds of data. Using the Stillinger-Weber potential we find excellent agreement between the two methods. We also use the direct method to calculate κ(T) for the Tersoff potential and find that the magnitude and the temperature-dependence are different for the two potentials and that neither potential agrees with experimental data. We argue that this implies that using the Stillinger-Weber or Tersoff potentials to predict trends in kappa as some system parameter is varied may yield results which are specific to the potential but not intrinsic to Si.

  18. Comparison of molecular dynamics methods and interatomic potentials for calculating the thermal conductivity of silicon.

    PubMed

    Howell, P C

    2012-12-14

    We compare the molecular dynamics Green-Kubo and direct methods for calculating thermal conductivity κ, using as a test case crystalline silicon at temperatures T in the range 500-1000 K (classical regime). We pay careful attention to the convergence with respect to simulation size and duration and to the procedures used to fit the simulation data. We show that in the Green-Kubo method the heat current autocorrelation function is characterized by three decay processes, of which the slowest lasts several tens of picoseconds so that convergence requires several tens of nanoseconds of data. Using the Stillinger-Weber potential we find excellent agreement between the two methods. We also use the direct method to calculate κ(T) for the Tersoff potential and find that the magnitude and the temperature-dependence are different for the two potentials and that neither potential agrees with experimental data. We argue that this implies that using the Stillinger-Weber or Tersoff potentials to predict trends in kappa as some system parameter is varied may yield results which are specific to the potential but not intrinsic to Si. PMID:23248991

  19. A modified W-W interatomic potential based on ab initio calculations

    NASA Astrophysics Data System (ADS)

    Wang, J.; Zhou, Y. L.; Li, M.; Hou, Q.

    2014-01-01

    In this paper we have developed a Finnis-Sinclair-type interatomic potential for W-W interactions that is based on ab initio calculations. The modified potential is able to reproduce the correct formation energies of self-interstitial atom (SIA) defects in tungsten, offering a significant improvement over the Ackland-Thetford tungsten potential. Using the modified potential, the thermal expansion is calculated in a temperature range from 0 to 3500 K. The results are in reasonable agreement with the experimental data, thus overcoming the shortcomings of the negative thermal expansion using the Derlet-Nguyen-Manh-Dudarev tungsten potential. The W-W potential presented here is also applied to study in detail the diffusion of SIAs in tungsten. We reveal that the initial SIA initiates a sequence of tungsten atom displacements and replacements in the <1 1 1> direction. An Arrhenius fit to the diffusion data at temperatures below 550 K indicates a migration energy of 0.022 eV, which is in reasonable agreement with the experimental data.

  20. Complex-scaling of screened Coulomb potentials for resonance calculations utilizing the modified Bessel functions

    NASA Astrophysics Data System (ADS)

    Jiao, Li-Guang; Ho, Yew Kam

    2014-05-01

    The screened Coulomb potential (SCP) has been extensively used in atomic physics, nuclear physics, quantum chemistry and plasma physics. However, an accurate calculation for atomic resonances under SCP is still a challenging task for various methods. Within the complex-scaling computational scheme, we have developed a method utilizing the modified Bessel functions to calculate doubly-excited resonances in two-electron atomic systems with configuration interaction-type basis. To test the validity of our method, we have calculated S- and P-wave resonance states of the helium atom with various screening strengths, and have found good agreement with earlier calculations using different methods. Our present method can be applied to calculate high-lying resonances associated with high excitation thresholds of the He+ ion, and with high-angular-momentum states. The derivation and calculation details of our present investigation together with new results of high-angular-momentum states will be presented at the meeting. Supported by NSC of Taiwan.

  1. Calculating excess chemical potentials using dynamic simulations in the fourth dimension

    SciTech Connect

    Pomes, R.; Eisenmesser, E.; Post, C.B.; Roux, B.

    1999-08-01

    A general method for computing excess chemical potentials is presented. The excess chemical potential of a solute or ligand molecule is estimated from the potential of mean-force (PMF) calculated along a nonphysical fourth spatial dimension, {ital w}, into which the molecule is gradually inserted or from which it is gradually abstracted. According to this {open_quotes}4D-PMF{close_quotes} (four dimensional) scheme, the free energy difference between two limiting states defines the excess chemical potential: At w={plus_minus}{infinity}, the molecule is not interacting with the rest of the system, whereas at w=0, it is fully interacting. Use of a fourth dimension avoids the numerical instability in the equations of motion encountered upon growing or shrinking solute atoms in conventional free energy perturbation simulations performed in three dimensions, while benefiting from the efficient sampling of configurational space afforded by PMF calculations. The applicability and usefulness of the method are illustrated with calculations of the hydration free energy of simple Lennard-Jones (LJ) solutes, a water molecule, and camphor, using molecular dynamics simulations and umbrella sampling. Physical insight into the nature of the PMF profiles is gained from a continuum treatment of short- and long-range interactions. The short-range barrier for dissolution of a LJ solute in the added dimension provides an apparent surface tension of the solute. An approximation to the long-range behavior of the PMF profiles is made in terms of a continuum treatment of LJ dispersion and electrostatic interactions. Such an analysis saves the need for configurational sampling in the long-range limit of the fourth dimension. The 4D-PMF method of calculating excess chemical potentials should be useful for neutral solute and ligand molecules with a wide range of sizes, shapes, and polarities. {copyright} {ital 1999 American Institute of Physics.}

  2. Regularizing the molecular potential in electronic structure calculations. II. Many-body methods

    SciTech Connect

    Bischoff, Florian A.

    2014-11-14

    In Paper I of this series [F. A. Bischoff, “Regularizing the molecular potential in electronic structure calculations. I. SCF methods,” J. Chem. Phys. 141, 184105 (2014)] a regularized molecular Hamilton operator for electronic structure calculations was derived and its properties in SCF calculations were studied. The regularization was achieved using a correlation factor that models the electron-nuclear cusp. In the present study we extend the regularization to correlated methods, in particular the exact solution of the two-electron problem, as well as second-order many body perturbation theory. The nuclear and electronic correlation factors lead to computations with a smaller memory footprint because the singularities are removed from the working equations, which allows coarser grid resolution while maintaining the precision. Numerical examples are given.

  3. Ionization potential of {sup 9}Be calculated including nuclear motion and relativistic corrections

    SciTech Connect

    Stanke, Monika; Kedziera, Dariusz; Bubin, Sergiy; Adamowicz, Ludwik

    2007-05-15

    Variational calculations employing explicitly correlated Gaussian functions have been performed for the ground states of {sup 9}Be and {sup 9}Be{sup +} including the nuclear motion [i.e., without assuming the Born-Oppenheimer (BO) approximation]. An approach based on the analytical energy gradient calculated with respect to the Gaussian exponential parameters was employed, leading to energies of the two systems noticeably improved over those found in the recent paper of Pachucki and Komasa [Phys. Rev. A 73, 052502 (2006)]. The non-BO wave functions were used to calculate the {alpha}{sup 2} relativistic corrections ({alpha}=e{sup 2}/({Dirac_h}/2{pi})c). With those corrections and the {alpha}{sup 3} and {alpha}{sup 4} corrections taken from Pachucki and Komasa, a new value of the ionization potential (IP) of {sup 9}Be was determined. It agrees very well with the most recent experimental IP.

  4. Classical calculation of the equilibrium constants for true bound dimers using complete potential energy surface

    SciTech Connect

    Buryak, Ilya; Vigasin, Andrey A.

    2015-12-21

    The present paper aims at deriving classical expressions which permit calculation of the equilibrium constant for weakly interacting molecular pairs using a complete multidimensional potential energy surface. The latter is often available nowadays as a result of the more and more sophisticated and accurate ab initio calculations. The water dimer formation is considered as an example. It is shown that even in case of a rather strongly bound dimer the suggested expression permits obtaining quite reliable estimate for the equilibrium constant. The reliability of our obtained water dimer equilibrium constant is briefly discussed by comparison with the available data based on experimental observations, quantum calculations, and the use of RRHO approximation, provided the latter is restricted to formation of true bound states only.

  5. Shear-deformation-potential constant of the conduction-band minima of Si: Pseudopotential calculations

    NASA Astrophysics Data System (ADS)

    Li, Ming-Fu; Gu, Zong-Quan; Wang, Jian-Qing

    1990-09-01

    We have calculated the value of the shear-deformation-potential constant Ξu of the conduction-band minima of Si and its temperature coefficient dΞu/dT. The value of Ξu is 9.0 eV for an ab initio pseudopotential calculation and 10.8 eV by the empirical-pseudopotential method (EPM), in good agreement with our experiment. The EPM calculations of the temperature dependence of Ξu yield the values of (dΞu/dT)||DW=-0.04 meV/K due to the Debye-Waller contribution, and (dΞu/dT)||TE=-0.04 meV/K for thermal expansion. We suspect and suggest that the existing experimental value of dΞu/dT~=+3 meV/K is unreliable due to large experimental uncertainty.

  6. Classical calculation of the equilibrium constants for true bound dimers using complete potential energy surface

    NASA Astrophysics Data System (ADS)

    Buryak, Ilya; Vigasin, Andrey A.

    2015-12-01

    The present paper aims at deriving classical expressions which permit calculation of the equilibrium constant for weakly interacting molecular pairs using a complete multidimensional potential energy surface. The latter is often available nowadays as a result of the more and more sophisticated and accurate ab initio calculations. The water dimer formation is considered as an example. It is shown that even in case of a rather strongly bound dimer the suggested expression permits obtaining quite reliable estimate for the equilibrium constant. The reliability of our obtained water dimer equilibrium constant is briefly discussed by comparison with the available data based on experimental observations, quantum calculations, and the use of RRHO approximation, provided the latter is restricted to formation of true bound states only.

  7. Artificial Bee Colony Optimization of Capping Potentials for Hybrid Quantum Mechanical/Molecular Mechanical Calculations.

    PubMed

    Schiffmann, Christoph; Sebastiani, Daniel

    2011-05-10

    We present an algorithmic extension of a numerical optimization scheme for analytic capping potentials for use in mixed quantum-classical (quantum mechanical/molecular mechanical, QM/MM) ab initio calculations. Our goal is to minimize bond-cleavage-induced perturbations in the electronic structure, measured by means of a suitable penalty functional. The optimization algorithm-a variant of the artificial bee colony (ABC) algorithm, which relies on swarm intelligence-couples deterministic (downhill gradient) and stochastic elements to avoid local minimum trapping. The ABC algorithm outperforms the conventional downhill gradient approach, if the penalty hypersurface exhibits wiggles that prevent a straight minimization pathway. We characterize the optimized capping potentials by computing NMR chemical shifts. This approach will increase the accuracy of QM/MM calculations of complex biomolecules. PMID:26610125

  8. The ground-state potential energy curve of the radium dimer from relativistic coupled cluster calculations

    NASA Astrophysics Data System (ADS)

    Teodoro, Tiago Quevedo; Haiduke, Roberto Luiz Andrade; Dammalapati, Umakanth; Knoop, Steven; Visscher, Lucas

    2015-08-01

    The potential energy curve for the ground-state of radium dimer (Ra2) is provided by means of atomic and molecular relativistic coupled cluster calculations. The short-range part of this curve is defined by an equilibrium bond length of 5.324 Å, a dissociation energy of 897 cm-1, and a harmonic vibrational frequency of 20.5 cm-1. The asymptotic behavior at large interatomic distances is characterized by the van der Waals coefficients C6 = 5.090 × 103, C8 = 6.978 × 105, and C10 = 8.786 × 107 atomic units. The two regions are matched in an analytical potential to provide a convenient representation for use in further calculations, for instance, to model cold collisions between radium atoms. This might become relevant in future experiments on ultracold, optically trapped, radioactive radium atoms that are used to search for a permanent electric dipole moment.

  9. Calculation of high-order virial coefficients for the square-well potential

    NASA Astrophysics Data System (ADS)

    Do, Hainam; Feng, Chao; Schultz, Andrew J.; Kofke, David A.; Wheatley, Richard J.

    2016-07-01

    Accurate virial coefficients BN(λ ,ɛ ) (where ɛ is the well depth) for the three-dimensional square-well and square-step potentials are calculated for orders N = 5 - 9 and well widths λ =1.1 -2.0 using a very fast recursive method. The efficiency of the algorithm is enhanced significantly by exploiting permutation symmetry and by storing integrands for reuse during the calculation. For N = 9 the storage requirements become sufficiently large that a parallel algorithm is developed. The methodology is general and is applicable to other discrete potentials. The computed coefficients are precise even near the critical temperature, and thus open up possibilities for analysis of criticality of the system, which is currently not accessible by any other means.

  10. Calculation of high-order virial coefficients for the square-well potential.

    PubMed

    Do, Hainam; Feng, Chao; Schultz, Andrew J; Kofke, David A; Wheatley, Richard J

    2016-07-01

    Accurate virial coefficients B_{N}(λ,ɛ) (where ɛ is the well depth) for the three-dimensional square-well and square-step potentials are calculated for orders N=5-9 and well widths λ=1.1-2.0 using a very fast recursive method. The efficiency of the algorithm is enhanced significantly by exploiting permutation symmetry and by storing integrands for reuse during the calculation. For N=9 the storage requirements become sufficiently large that a parallel algorithm is developed. The methodology is general and is applicable to other discrete potentials. The computed coefficients are precise even near the critical temperature, and thus open up possibilities for analysis of criticality of the system, which is currently not accessible by any other means. PMID:27575230

  11. Microscopic calculation of {alpha}-decay half-lives with a deformed potential

    SciTech Connect

    Ni Dongdong; Ren Zhongzhou

    2009-11-15

    A new version of the generalized density-dependent cluster model is presented to describe an {alpha} particle tunneling through a deformed potential barrier. The microscopic deformed potential is numerically constructed in the double-folding model by the multipole expansion method. The decay width is computed using the coupled-channel Schroedinger equation with outgoing wave boundary conditions. We perform a systematic calculation on {alpha}-decay half-lives of even-even nuclei ranging from Z=52 to Z=104, including 65 well-deformed ones. The calculated {alpha}-decay half-lives are found to be in good agreement with the experimental values. There also exists good agreement with the available experimental branching ratios for well-deformed systems.

  12. Multireference configuration interaction calculations of the first six ionization potentials of the uranium atom

    SciTech Connect

    Bross, David H.; Parmar, Payal; Peterson, Kirk A.

    2015-11-14

    The first 6 ionization potentials (IPs) of the uranium atom have been calculated using multireference configuration interaction (MRCI+Q) with extrapolations to the complete basis set limit using new all-electron correlation consistent basis sets. The latter was carried out with the third-order Douglas-Kroll-Hess Hamiltonian. Correlation down through the 5s5p5d electrons has been taken into account, as well as contributions to the IPs due to the Lamb shift. Spin-orbit coupling contributions calculated at the 4-component Kramers restricted configuration interaction level, as well as the Gaunt term computed at the Dirac-Hartree-Fock level, were added to the best scalar relativistic results. The final ionization potentials are expected to be accurate to at least 5 kcal/mol (0.2 eV) and thus more reliable than the current experimental values of IP{sub 3} through IP{sub 6}.

  13. Phenomenological calculation of nuclear binding energy and density with Yukawa-potentials

    NASA Astrophysics Data System (ADS)

    Scheid, W.

    2016-01-01

    In this paper, we study a phenomenological collective model for the calculation of the nuclear density and ground state binding energy of nuclei. The proton density is assumed proportional to the nuclear density. The total binding energy of the nuclear matter consists of the binding energy of infinite nuclear matter, of two Yukawa-potentials, of the Coulomb-energy and of the symmetry-energy. The parameters of the Yukawa-potential are fitted with the Bethe-Weizsäcker (BW) mass formula. The resulting binding energies and nuclear densities agree quite satisfying with known nuclear values.

  14. Thermal conductivity of silicene calculated using an optimized Stillinger-Weber potential

    NASA Astrophysics Data System (ADS)

    Zhang, Xiaoliang; Xie, Han; Hu, Ming; Bao, Hua; Yue, Shengying; Qin, Guangzhao; Su, Gang

    2014-02-01

    Silicene, the silicon-based counterpart of graphene with a two-dimensional honeycomb lattice, has attracted tremendous interest both theoretically and experimentally due to its significant potential industrial applications. From the aspect of theoretical study, the widely used classical molecular dynamics simulation is an appropriate way to investigate the transport phenomena and mechanisms in nanostructures such as silicene. Unfortunately, no available interatomic potential can precisely characterize the unique features of silicene. Here, we optimized the Stillinger-Weber potential parameters specifically for a single-layer Si sheet, which can accurately reproduce the low buckling structure of silicene and the full phonon dispersion curves obtained from ab initio calculations. By performing equilibrium and nonequilibrium molecular dynamics simulations and anharmonic lattice dynamics calculations with the new potential, we reveal that the three methods consistently yield an extremely low thermal conductivity of silicene and a short phonon mean-free path, suggesting silicene as a potential candidate for high-efficiency thermoelectric materials. Moreover, by qualifying the relative contributions of lattice vibrations in different directions, we found that the longitudinal phonon modes dominate the thermal transport in silicene, which is fundamentally different from graphene, despite the similarity of their two-dimensional honeycomb lattices.

  15. Three-body calculations for the K ‑ pp system within potential models

    NASA Astrophysics Data System (ADS)

    Kezerashvili, R. Ya; Tsiklauri, S. M.; Filikhin, I.; Suslov, V. M.; Vlahovic, B.

    2016-06-01

    We present three-body nonrelativistic calculations within the framework of a potential model for the kaonic cluster K ‑ pp using two methods: the method of hyperspherical harmonics in the momentum representation and the method of Faddeev equations in configuration space. To perform numerical calculations, different NN and antikaon–nucleon interactions are applied. The results of the calculations for the ground-state energy for the K ‑ pp system obtained by both methods are in reasonable agreement. Although the ground-state energy is not sensitive to the pp interaction, it shows very strong dependence on the K ‑ p potential. We show that the dominant clustering of the {K}-{pp} system in the configuration Λ (1405) + p allows us to calculate the binding energy to good accuracy within a simple cluster approach for the differential Faddeev equations. The theoretical discrepancies in the binding energy and width for the K ‑ pp system related to the different pp and K ‑ p interactions are addressed.

  16. Calculations of {sup 8}He+p elastic cross sections using a microscopic optical potential

    SciTech Connect

    Lukyanov, V. K.; Zemlyanaya, E. V.; Lukyanov, K. V.; Kadrev, D. N.; Antonov, A. N.; Gaidarov, M. K.; Massen, S. E.

    2009-08-15

    An approach to calculate microscopic optical potential with the real part obtained by a folding procedure and with the imaginary part inherent in the high-energy approximation is applied to study the {sup 8}He+p elastic-scattering data at energies of tens of MeV/nucleon. The neutron and proton density distributions obtained in different models for {sup 8}He are used in the calculations of the differential cross sections. The role of the spin-orbit potential is studied. Comparison of the calculations with the available experimental data on the elastic-scattering differential cross sections at beam energies of 15.7, 26.25, 32, 66, and 73 MeV/nucleon is performed. The problem of the ambiguities of the depths of each component of the optical potential is considered by means of the imposed physical criterion related to the known behavior of the volume integrals as functions of the incident energy. It is shown also that the role of the surface absorption is rather important, in particular for the lowest incident energies (e.g., 15.7 and 26.25 MeV/nucleon)

  17. Evaluation of uncertainty in gravity wave potential energy calculations through GPS radio occultation measurements

    NASA Astrophysics Data System (ADS)

    Luna, D.; Alexander, P.; de la Torre, A.

    2013-09-01

    The application of the Global Positioning System (GPS) radio occultation (RO) method to the atmosphere enables the determination of height profiles of temperature, among other variables. From these measurements, gravity wave activity is usually quantified by calculating the potential energy through the integration of the ratio of perturbation and background temperatures between two given altitudes in each profile. The uncertainty in the estimation of wave activity depends on the systematic biases and random errors of the measured temperature, but also on additional factors like the selected vertical integration layer and the separation method between background and perturbation temperatures. In this study, the contributions of different parameters and variables to the uncertainty in the calculation of gravity wave potential energy in the lower stratosphere are investigated and quantified. In particular, a Monte Carlo method is used to evaluate the uncertainty that results from different GPS RO temperature error distributions. In addition, our analysis shows that RO data above 30 km height becomes dubious for gravity waves potential energy calculations.

  18. Medication calculation: the potential role of digital game-based learning in nurse education.

    PubMed

    Foss, Brynjar; Mordt Ba, Petter; Oftedal, Bjørg F; Løkken, Atle

    2013-12-01

    Medication dose calculation is one of several medication-related activities that are conducted by nurses daily. However, medication calculation skills appear to be an area of global concern, possibly because of low numeracy skills, test anxiety, low self-confidence, and low self-efficacy among student nurses. Various didactic strategies have been developed for student nurses who still lack basic mathematical competence. However, we suggest that the critical nature of these skills demands the investigation of alternative and/or supplementary didactic approaches to improve medication calculation skills and to reduce failure rates. Digital game-based learning is a possible solution because of the following reasons. First, mathematical drills may improve medication calculation skills. Second, games are known to be useful during nursing education. Finally, mathematical drill games appear to improve the attitudes of students toward mathematics. The aim of this article was to discuss common challenges of medication calculation skills in nurse education, and we highlight the potential role of digital game-based learning in this area. PMID:24107685

  19. Calculation of water drop trajectories to and about arbitrary three-dimensional lifting and nonlifting bodies in potential airflow

    NASA Technical Reports Server (NTRS)

    Norment, H. G.

    1985-01-01

    Subsonic, external flow about nonlifting bodies, lifting bodies or combinations of lifting and nonlifting bodies is calculated by a modified version of the Hess lifting code. Trajectory calculations can be performed for any atmospheric conditions and for all water drop sizes, from the smallest cloud droplet to large raindrops. Experimental water drop drag relations are used in the water drop equations of motion and effects of gravity settling are included. Inlet flow can be accommodated, and high Mach number compressibility effects are corrected for approximately. Seven codes are described: (1) a code used to debug and plot body surface description data; (2) a code that processes the body surface data to yield the potential flow field; (3) a code that computes flow velocities at arrays of points in space; (4) a code that computes water drop trajectories from an array of points in space; (5) a code that computes water drop trajectories and fluxes to arbitrary target points; (6) a code that computes water drop trajectories tangent to the body; and (7) a code that produces stereo pair plots which include both the body and trajectories. Accuracy of the calculations is discussed, and trajectory calculation results are compared with prior calculations and with experimental data.

  20. Bounces and the calculation of quantum tunneling effects

    NASA Astrophysics Data System (ADS)

    Liang, Jiu-Qing; Müller-Kirsten, H. J. W.

    1992-04-01

    The imaginary part of the energy of the metastable ground state for the inverted double-well potential is calculated by using the path-integral method. The tunneling process is dominated by bounces. It is shown that the evaluation of the determinant of the second variation of the action at the bounce can be avoided, and that the imaginary part of the energy results directly from characteristic properties of the bounce itself, namely, the antisymmetry of its first time derivative under time reversal. The imaginary part of the result is in exact agreement with that of the well-known WKB calculation of Bender and Wu.

  1. Electronic Structure Calculations of Inter-Ring Torsional Potentials of Regioregular Poly (3-METHYL Thiophene) Oligomers

    NASA Astrophysics Data System (ADS)

    Bhatta, Ram S.; Perry, David S.

    2010-06-01

    The inter-ring torsional potentials of poly (3-methyl thiophene) (P3MT) oligomers are investigated by means of electronic structure calculations. Single layer and ONIOM calculations were performed at B3LYP level with 6-31++G(d,p) basis on the partially optimized geometries of dimer, tetramer and hexamer of P3MT oligomers. Potential energy surfaces are computed as a function of the multiple inter-ring torsional angles involved. The following conclusions are reached: (i) A mixture of cis and trans geometries can be expected in a disordered polymer. (ii) The cis-trans barrier is low enough to allow cis-trans conversion at room temperature. (iii) In the dimer, the potential energy minima are about 30^0 from the cis and trans planar geometries, but planar geometries are stabilized as the chain length increases. (iv) The extended conjugation causes the torsional potential about one inter-ring bond to be coupled to other torsions along the oligomer chain.

  2. Quantum Calculations on Salt Bridges with Water: Potentials, Structure, and Properties

    SciTech Connect

    Liao, Sing; Green, Michael E.

    2011-01-01

    Salt bridges are electrostatic links between acidic and basic amino acids in a protein; quantum calculations are used here to determine the energetics and other properties of one form of these species, in the presence of water molecules. The acidic groups are carboxylic acids (aspartic and glutamic acids); proteins have two bases with pK above physiological pH: one, arginine, with a guanidinium basic group, the other lysine, which is a primary amine. Only arginine is modeled here, by ethyl guanidinium, while propionic acid is used as a model for either carboxylic acid. The salt bridges are accompanied by 0-12 water molecules; for each of the 13 systems, the energy-bond distance relation, natural bond orbitals (NBO), frequency calculations allowing thermodynamic corrections to room temperature, and dielectric constant dependence, were all calculated. The water molecules were found to arrange themselves in hydrogen bonded rings anchored to the oxygens of the salt bridge components. This was not surprising in itself, but it was found that the rings lead to a periodicity in the energy, and to a 'water addition' rule. The latter shows that the initial rings, with four oxygen atoms, become five member rings when an additional water molecule becomes available, with the additional water filling in at the bond with the lowest Wiberg index, as calculated using NBO. The dielectric constant dependence is the expected hyperbola, and the fit of the energy to the inverse dielectric constant is determined. There is an energy periodicity related to ring formation upon addition of water molecules. When 10 water molecules have been added, all spaces near the salt bridge are filled, completing the first hydration shell, and a second shell starts to form. The potentials associated with salt bridges depend on their hydration, and potentials assigned without regard to local hydration are likely to cause errors as large as or larger than kBT, thus suggesting a serious problem if these

  3. Regularizing the molecular potential in electronic structure calculations. I. SCF methods

    SciTech Connect

    Bischoff, Florian A.

    2014-11-14

    We present a method to remove the singular nuclear potential in a molecule and replace it with a regularized potential that is more amenable to be represented numerically. The singular nuclear potential is canceled by the similarity-transformed kinetic energy operator giving rise to an effective nuclear potential that contains derivative operators acting on the wave function. The method is fully equivalent to the non-similarity-transformed version. We give numerical examples within the framework of multi-resolution analysis for medium-sized molecules.

  4. Effective UV radiation from model calculations and measurements

    NASA Technical Reports Server (NTRS)

    Feister, Uwe; Grewe, Rolf

    1994-01-01

    Model calculations have been made to simulate the effect of atmospheric ozone and geographical as well as meteorological parameters on solar UV radiation reaching the ground. Total ozone values as measured by Dobson spectrophotometer and Brewer spectrometer as well as turbidity were used as input to the model calculation. The performance of the model was tested by spectroradiometric measurements of solar global UV radiation at Potsdam. There are small differences that can be explained by the uncertainty of the measurements, by the uncertainty of input data to the model and by the uncertainty of the radiative transfer algorithms of the model itself. Some effects of solar radiation to the biosphere and to air chemistry are discussed. Model calculations and spectroradiometric measurements can be used to study variations of the effective radiation in space in space time. The comparability of action spectra and their uncertainties are also addressed.

  5. Potential Ergogenic Effects of Saffron.

    PubMed

    Meamarbashi, Abbas; Rajabi, Ali

    2016-01-01

    Crocus sativus, commonly known as saffron, is a rich source of carotenoids with many health benefits. The muscular strength, pulmonary function, and reaction time are vital to the athlete's performance, and this study aimed to investigate an ergogenic effect of saffron. Twenty-eight nonactive and healthy male university students were randomly assigned into the saffron (n = 14) and control (n = 15) groups. The experimental group received dried saffron stigma (300 mg/day for 10 days) and the control group received a placebo. After one session, familiarization with the tests, anthropometric parameters, visual and audio reaction times, and the maximum isometric and isotonic forces on a leg press machine were measured accordingly, 1 day before and after the supplementation period. This study shows that 10 days of supplementation with saffron significantly increased (10.1%) the isometric force (p < .0001; effect size (EF) = 0.432) and increased 6.1% the isotonic force (p < .0001; effect size = 0.662), as well as effecting faster visual (p < .05; EF = 0.217) and audio (p < .05; EF = 0.214) reaction times. The ergogenic effect of saffron (increase in the forces) may contribute to increase in the muscle mitochondrial biogenesis and positive effect on the motor cortex, both of which may explain faster audio and visual reaction times. Saffron supplementation was also possibly responsible for improvement of muscle blood perfusion and facilitation in the oxygen transport. PMID:26811090

  6. The effect of dose calculation accuracy on inverse treatment planning

    NASA Astrophysics Data System (ADS)

    Jeraj, Robert; Keall, Paul J.; Siebers, Jeffrey V.

    2002-02-01

    The effect of dose calculation accuracy during inverse treatment planning for intensity modulated radiotherapy (IMRT) was studied in this work. Three dose calculation methods were compared: Monte Carlo, superposition and pencil beam. These algorithms were used to calculate beamlets, which were subsequently used by a simulated annealing algorithm to determine beamlet weights which comprised the optimal solution to the objective function. Three different cases (lung, prostate and head and neck) were investigated and several different objective functions were tested for their effect on inverse treatment planning. It is shown that the use of inaccurate dose calculation introduces two errors in a treatment plan, a systematic error and a convergence error. The systematic error is present because of the inaccuracy of the dose calculation algorithm. The convergence error appears because the optimal intensity distribution for inaccurate beamlets differs from the optimal solution for the accurate beamlets. While the systematic error for superposition was found to be ~1% of Dmax in the tumour and slightly larger outside, the error for the pencil beam method is typically ~5% of Dmax and is rather insensitive to the given objectives. On the other hand, the convergence error was found to be very sensitive to the objective function, is only slightly correlated to the systematic error and should be determined for each case individually. Our results suggest that because of the large systematic and convergence errors, inverse treatment planning systems based on pencil beam algorithms alone should be upgraded either to superposition or Monte Carlo based dose calculations.

  7. Theoretical calculations and vibrational potential energy surface of 4-silaspiro(3,3)heptane

    SciTech Connect

    Ocola, Esther J.; Medders, Cross; Laane, Jaan; Meinander, Niklas

    2014-04-28

    Theoretical computations have been carried out on 4-silaspiro(3,3)heptane (SSH) in order to calculate its molecular structure and conformational energies. The molecule has two puckered four-membered rings with dihedral angles of 34.2° and a tilt angle of 9.4° between the two rings. Energy calculations were carried out for different conformations of SSH. These results allowed the generation of a two-dimensional ring-puckering potential energy surface (PES) of the form V = a(x{sub 1}{sup 4} + x{sub 2}{sup 4}) – b(x{sub 1}{sup 2} + x{sub 2}{sup 2}) + cx{sub 1}{sup 2}x{sub 2}{sup 2}, where x{sub 1} and x{sub 2} are the ring-puckering coordinates for the two rings. The presence of sufficiently high potential energy barriers prevents the molecule from undergoing pseudorotation. The quantum states, wave functions, and predicted spectra resulting from the PESs were calculated.

  8. Application of a sitting MIRD phantom for effective dose calculations.

    PubMed

    Olsher, Richard H; Van Riper, Kenneth A

    2005-01-01

    In typical realistic scenarios, dose factors due to 60Co contaminated steel, used in consumer products, cannot be approximated by standard exposure geometries. It is then necessary to calculate the effective dose using an appropriate anthropomorphic phantom. MCNP calculations were performed using a MIRD human model in two settings. In the first, a male office worker is sitting in a chair containing contaminated steel, surrounded by contaminated furniture. In the second, a male driver is seated inside an automobile, the steel of which is uniformly contaminated. To accurately calculate the dose to lower body organs, especially the gonads, it was essential to modify the MIRD model to simulate two sitting postures: chair and driving position. The phantom modifications are described, and the results of the calculations are presented. In the case of the automobile scenarios, results are compared to those obtained using an isotropic fluence-to-dose conversion function. PMID:16604666

  9. Electronic coupling calculations with transition charges, dipoles, and quadrupoles derived from electrostatic potential fitting.

    PubMed

    Fujimoto, Kazuhiro J

    2014-12-01

    A transition charge, dipole, and quadrupole from electrostatic potential (TrESP-CDQ) method for electronic coupling calculations is proposed. The TrESP method is based on the classical description of electronic Coulomb interaction between transition densities for individual molecules. In the original TrESP method, only the transition charge interactions were considered as the electronic coupling. In the present study, the TrESP method is extended to include the contributions from the transition dipoles and quadrupoles as well as the transition charges. Hence, the self-consistent transition density is employed in the ESP fitting procedure. To check the accuracy of the present approach, several test calculations are performed to a helium dimer, a methane dimer, and an ethylene dimer. As a result, the TrESP-CDQ method gives a much improved description of the electronic coupling, compared with the original TrESP method. The calculated results also show that the self-consistent treatment to the transition densities contributes significantly to the accuracy of the electronic coupling calculations. Based on the successful description of the electronic coupling, the contributions to the electronic coupling are also analyzed. This analysis clearly shows a negligible contribution of the transition charge interaction to the electronic coupling. Hence, the distribution of the transition density is found to strongly influence the magnitudes of the transition charges, dipoles, and quadrupoles. The present approach is useful for analyzing and understanding the mechanism of excitation-energy transfer. PMID:25481127

  10. Calculation of subsonic and supersonic steady and unsteady aerodynamic forces using velocity potential aerodynamic elements

    NASA Technical Reports Server (NTRS)

    Haviland, J. K.; Yoo, Y. S.

    1976-01-01

    Expressions for calculation of subsonic and supersonic, steady and unsteady aerodynamic forces are derived, using the concept of aerodynamic elements applied to the downwash velocity potential method. Aerodynamic elements can be of arbitrary out of plane polygon shape, although numerical calculations are restricted to rectangular elements, and to the steady state case in the supersonic examples. It is suggested that the use of conforming, in place of rectangular elements, would give better results. Agreement with results for subsonic oscillating T tails is fair, but results do not converge as the number of collocation points is increased. This appears to be due to the form of expression used in the calculations. The methods derived are expected to facilitate automated flutter analysis on the computer. In particular, the aerodynamic element concept is consistent with finite element methods already used for structural analysis. The method is universal for the complete Mach number range, and, finally, the calculations can be arranged so that they do not have to be repeated completely for every reduced frequency.

  11. Electronic coupling calculations with transition charges, dipoles, and quadrupoles derived from electrostatic potential fitting

    SciTech Connect

    Fujimoto, Kazuhiro J.

    2014-12-07

    A transition charge, dipole, and quadrupole from electrostatic potential (TrESP-CDQ) method for electronic coupling calculations is proposed. The TrESP method is based on the classical description of electronic Coulomb interaction between transition densities for individual molecules. In the original TrESP method, only the transition charge interactions were considered as the electronic coupling. In the present study, the TrESP method is extended to include the contributions from the transition dipoles and quadrupoles as well as the transition charges. Hence, the self-consistent transition density is employed in the ESP fitting procedure. To check the accuracy of the present approach, several test calculations are performed to a helium dimer, a methane dimer, and an ethylene dimer. As a result, the TrESP-CDQ method gives a much improved description of the electronic coupling, compared with the original TrESP method. The calculated results also show that the self-consistent treatment to the transition densities contributes significantly to the accuracy of the electronic coupling calculations. Based on the successful description of the electronic coupling, the contributions to the electronic coupling are also analyzed. This analysis clearly shows a negligible contribution of the transition charge interaction to the electronic coupling. Hence, the distribution of the transition density is found to strongly influence the magnitudes of the transition charges, dipoles, and quadrupoles. The present approach is useful for analyzing and understanding the mechanism of excitation-energy transfer.

  12. Mesh size and code option effects of strength calculations

    SciTech Connect

    Kaul, Ann M

    2010-12-10

    Modern Lagrangian hydrodynamics codes include numerical methods which allow calculations to proceed past the point obtainable by a purely Lagrangian scheme. These options can be employed as the user deems necessary to 'complete' a calculation. While one could argue that any calculation is better than none, to truly understand the calculated results and their relationship to physical reality, the user needs to understand how their runtime choices affect the calculated results. One step toward this goal is to understand the effect of each runtime choice on particular pieces of the code physics. This paper will present simulation results for some experiments typically used for strength model validation. Topics to be covered include effect of mesh size, use of various ALE schemes for mesh detangling, and use of anti-hour-glassing schemes. Experiments to be modeled include the lower strain rate ({approx} 10{sup 4} s{sup -1}) gas gun driven Taylor impact experiments and the higher strain rate ({approx} 10{sup 5}-10{sup 6} s{sup -1}) HE products driven perturbed plate experiments. The necessary mesh resolution and the effect of the code runtime options are highly dependent on the amount of localization of strain and stress in each experiment. In turn, this localization is dependent on the geometry of the experimental setup and the drive conditions.

  13. Determination of a silane intermolecular force field potential model from an ab initio calculation

    SciTech Connect

    Li, Arvin Huang-Te; Chao, Sheng D.; Chang, Chien-Cheng

    2010-12-15

    Intermolecular interaction potentials of the silane dimer in 12 orientations have been calculated by using the Hartree-Fock (HF) self-consistent theory and the second-order Moeller-Plesset (MP2) perturbation theory. We employed basis sets from Pople's medium-size basis sets [up to 6-311++G(3df, 3pd)] and Dunning's correlation consistent basis sets (up to the triply augmented correlation-consistent polarized valence quadruple-zeta basis set). We found that the minimum energy orientations were the G and H conformers. We have suggested that the Si-H attractions, the central silicon atom size, and electronegativity play essential roles in weakly binding of a silane dimer. The calculated MP2 potential data were employed to parametrize a five-site force field for molecular simulations. The Si-Si, Si-H, and H-H interaction parameters in a pairwise-additive, site-site potential model for silane molecules were regressed from the ab initio energies.

  14. Calculations of the ionization potentials of the halogens by the relativistic Hartree-Rock-Dirac method taking account of superposition of configurations

    SciTech Connect

    Tupitsyn, I.I.

    1988-03-01

    The ionization potentials of the halogen group have been calculated. The calculations were carried out using the relativistic Hartree-Fock method taking into account correlation effects. Comparison of theoretical results with experimental data for the elements F, Cl, Br, and I allows an estimation of the accuracy and reliability of the method. The theoretical values of the ionization potential of astatine obtained here may be of definite interest for the chemistry of astatine.

  15. Using Electronic Energy Derivative Information in Automated Potential Energy Surface Construction for Vibrational Calculations.

    PubMed

    Sparta, Manuel; Hansen, Mikkel B; Matito, Eduard; Toffoli, Daniele; Christiansen, Ove

    2010-10-12

    The availability of an accurate representation of the potential energy surface (PES) is an essential prerequisite in an anharmonic vibrational calculation. At the same time, the high dimensionality of the fully coupled PES and the adverse scaling properties with respect to the molecular size make the construction of an accurate PES a computationally demanding task. In the past few years, our group tested and developed a series of tools and techniques aimed at defining computationally efficient, black-box protocols for the construction of PESs for use in vibrational calculations. This includes the definition of an adaptive density-guided approach (ADGA) for the construction of PESs from an automatically generated set of evaluation points. Another separate aspect has been the exploration of the use of derivative information through modified Shepard (MS) interpolation/extrapolation procedures. With this article, we present an assembled machinery where these methods are embedded in an efficient way to provide both a general machinery as well as concrete computational protocols. In this framework we introduce and discuss the accuracy and computational efficiency of two methods, called ADGA[2gx3M] and ADGA[2hx3M], where the ADGA recipe is used (with MS interpolation) to automatically define modest sized grids for up to two-mode couplings, while MS extrapolation based on, respectively, gradients only and gradients and Hessians from the ADGA determined points provides access to sufficiently accurate three-mode couplings. The performance of the resulting potentials is investigated in vibrational coupled cluster (VCC) calculations. Three molecular systems serve as benchmarks: a trisubstituted methane (CHFClBr), methanimine (CH2NH), and oxazole (C3H3NO). Furthermore, methanimine and oxazole are addressed in accurate calculations aiming to reproduce experimental results. PMID:26616778

  16. Semi-continuum electrostatic calculations of redox potentials in photosystem I.

    PubMed

    Ptushenko, Vasily V; Cherepanov, Dmitry A; Krishtalik, Lev I; Semenov, Alexey Yu

    2008-07-01

    The midpoint redox potentials (E(m)) of all cofactors in photosystem I from Synechococcus elongatus as well as of the iron-sulfur (Fe(4)S(4)) clusters in two soluble ferredoxins from Azotobacter vinelandii and Clostridium acidiurici were calculated within the framework of a semi-continuum dielectric approach. The widely used treatment of proteins as uniform media with single dielectric permittivity is oversimplified, particularly, because permanent charges are considered both as a source for intraprotein electric field and as a part of dielectric polarizability. Our approach overcomes this inconsistency by using two dielectric constants: optical epsilon(o)=2.5 for permanent charges pre-existing in crystal structure, and static epsilon(s) for newly formed charges. We also take into account a substantial dielectric heterogeneity of photosystem I revealed by photoelectric measurements and a liquid junction potential correction for E(m) values of relevant redox cofactors measured in aprotic solvents. We show that calculations based on a single permittivity have the discrepancy with experimental data larger than 0.7 V, whereas E(m) values calculated within our approach fall in the range of experimental estimates. The electrostatic analysis combined with quantum chemistry calculations shows that (i) the energy decrease upon chlorophyll dimerization is essential for the downhill mode of primary charge separation between the special pair P(700) and the primary acceptor A(0); (ii) the primary donor is apparently P(700) but not a pair of accessory chlorophylls; (iii) the electron transfer from the A branch quinone Q(A) to the iron-sulfur cluster F(X) is most probably downhill, whereas that from the B branch quinone Q(B) to F(X) is essentially downhill. PMID:18483776

  17. Intrinsic Defect Properties in GaN Calculated By Ab Initio and Empirical Potential Methods

    SciTech Connect

    Gao, Fei; Bylaska, Eric J.; Weber, William J.

    2004-12-13

    Density functional theory (DFT) has been used to investigate the formation, properties, and atomic configurations of vacancies, antisite defects and interstitials in GaN, and the DFT results are compared with those calculated by molecular dynamics (MD) simulations using two representative potentials. The DFT calculations reveal that the relaxation of vacancies is generally small, but the relaxation around antisite defects is large, especially for the Ga antisite that is not stable and converts to a N+-N<0001> split interstitial plus a Ga vacancy at the original site. The N interstitials, starting from all possible sites, eventually relax into a N+-N<11-20> split interstitial. In the case of Ga interstitials, the most stable configuration is a Ga octahedral interstitial, but the energy difference among all the interstitials is small. The Ga+-Ga<11-20> split interstitial can bridge the gap between non-bonded Ga atoms, thereby leading to a chain of four Ga atoms along the <11-20> direction in GaN. The formation energies of vacancies and antisite defects obtained using the Stillinger-Weber potential (SW) are in reasonable agreement with those obtained by DFT calculations, whereas the Tersoff-Brenner (TB) potential better describes the behavior of N interstitials. In the case of Ga interstitials, the most stable configuration predicted by the TB-model is a Ga+-N<11-20> split interstitial; while for the SW-model the Ga tetrahedral configuration is more stable, which is in contrast to DFT results.

  18. Effect size estimates: current use, calculations, and interpretation.

    PubMed

    Fritz, Catherine O; Morris, Peter E; Richler, Jennifer J

    2012-02-01

    The Publication Manual of the American Psychological Association (American Psychological Association, 2001, American Psychological Association, 2010) calls for the reporting of effect sizes and their confidence intervals. Estimates of effect size are useful for determining the practical or theoretical importance of an effect, the relative contributions of factors, and the power of an analysis. We surveyed articles published in 2009 and 2010 in the Journal of Experimental Psychology: General, noting the statistical analyses reported and the associated reporting of effect size estimates. Effect sizes were reported for fewer than half of the analyses; no article reported a confidence interval for an effect size. The most often reported analysis was analysis of variance, and almost half of these reports were not accompanied by effect sizes. Partial η2 was the most commonly reported effect size estimate for analysis of variance. For t tests, 2/3 of the articles did not report an associated effect size estimate; Cohen's d was the most often reported. We provide a straightforward guide to understanding, selecting, calculating, and interpreting effect sizes for many types of data and to methods for calculating effect size confidence intervals and power analysis. PMID:21823805

  19. Ranking protein-protein docking results using steered molecular dynamics and potential of mean force calculations.

    PubMed

    Kingsley, Laura J; Esquivel-Rodríguez, Juan; Yang, Ying; Kihara, Daisuke; Lill, Markus A

    2016-07-01

    Crystallization of protein-protein complexes can often be problematic and therefore computational structural models are often relied on. Such models are often generated using protein-protein docking algorithms, where one of the main challenges is selecting which of several thousand potential predictions represents the most near-native complex. We have developed a novel technique that involves the use of steered molecular dynamics (sMD) and umbrella sampling to identify near-native complexes among protein-protein docking predictions. Using this technique, we have found a strong correlation between our predictions and the interface RMSD (iRMSD) in ten diverse test systems. On two of the systems, we investigated if the prediction results could be further improved using potential of mean force calculations. We demonstrated that a near-native (<2.0 Å iRMSD) structure could be identified in the top-1 ranked position for both systems. © 2016 Wiley Periodicals, Inc. PMID:27232548

  20. Ab initio calculation of thermodynamic potentials and entropies for superionic water

    DOE PAGESBeta

    French, Martin; Desjarlais, Michael P.; Redmer, Ronald

    2016-02-25

    We construct thermodynamic potentials for two superionic phases of water [with body-centered cubic (bcc) and face-centered cubic (fcc) oxygen lattice] using a combination of density functional theory (DFT) and molecular dynamics simulations (MD). For this purpose, a generic expression for the free energy of warm dense matter is developed and parametrized with equation of state data from the DFT-MD simulations. A second central aspect is the accurate determination of the entropy, which is done using an approximate two-phase method based on the frequency spectra of the nuclear motion. The boundary between the bcc superionic phase and the ices VII andmore » X calculated with thermodynamic potentials from DFT-MD is consistent with that directly derived from the simulations. As a result, differences in the physical properties of the bcc and fcc superionic phases and their impact on interior modeling of water-rich giant planets are discussed.« less

  1. Ab initio calculation of thermodynamic potentials and entropies for superionic water

    NASA Astrophysics Data System (ADS)

    French, Martin; Desjarlais, Michael P.; Redmer, Ronald

    2016-02-01

    We construct thermodynamic potentials for two superionic phases of water [with body-centered cubic (bcc) and face-centered cubic (fcc) oxygen lattice] using a combination of density functional theory (DFT) and molecular dynamics simulations (MD). For this purpose, a generic expression for the free energy of warm dense matter is developed and parametrized with equation of state data from the DFT-MD simulations. A second central aspect is the accurate determination of the entropy, which is done using an approximate two-phase method based on the frequency spectra of the nuclear motion. The boundary between the bcc superionic phase and the ices VII and X calculated with thermodynamic potentials from DFT-MD is consistent with that directly derived from the simulations. Differences in the physical properties of the bcc and fcc superionic phases and their impact on interior modeling of water-rich giant planets are discussed.

  2. Ab initio calculation of thermodynamic potentials and entropies for superionic water.

    PubMed

    French, Martin; Desjarlais, Michael P; Redmer, Ronald

    2016-02-01

    We construct thermodynamic potentials for two superionic phases of water [with body-centered cubic (bcc) and face-centered cubic (fcc) oxygen lattice] using a combination of density functional theory (DFT) and molecular dynamics simulations (MD). For this purpose, a generic expression for the free energy of warm dense matter is developed and parametrized with equation of state data from the DFT-MD simulations. A second central aspect is the accurate determination of the entropy, which is done using an approximate two-phase method based on the frequency spectra of the nuclear motion. The boundary between the bcc superionic phase and the ices VII and X calculated with thermodynamic potentials from DFT-MD is consistent with that directly derived from the simulations. Differences in the physical properties of the bcc and fcc superionic phases and their impact on interior modeling of water-rich giant planets are discussed. PMID:26986321

  3. Modeling intermolecular interactions of physisorbed organic molecules using pair potential calculations

    SciTech Connect

    Kroeger, Ingo; Stadtmueller, Benjamin; Wagner, Christian; Weiss, Christian; Temirov, Ruslan; Tautz, F. Stefan; Kumpf, Christian

    2011-12-21

    The understanding and control of epitaxial growth of organic thin films is of crucial importance in order to optimize the performance of future electronic devices. In particular, the start of the submonolayer growth plays an important role since it often determines the structure of the first layer and subsequently of the entire molecular film. We have investigated the structure formation of 3,4,9,10-perylene-tetracarboxylic dianhydride and copper-phthalocyanine molecules on Au(111) using pair-potential calculations based on van der Waals and electrostatic intermolecular interactions. The results are compared with the fundamental lateral structures known from experiment and an excellent agreement was found for these weakly interacting systems. Furthermore, the calculations are even suitable for chemisorptive adsorption as demonstrated for copper-phthalocyanine/Cu(111), if the influence of charge transfer between substrate and molecules is known and the corresponding charge redistribution in the molecules can be estimated. The calculations are of general applicability for molecular adsorbate systems which are dominated by electrostatic and van der Waals interaction.

  4. Potential energy surface and second virial coefficient of methane-water from ab initio calculations

    NASA Astrophysics Data System (ADS)

    Akin-Ojo, Omololu; Szalewicz, Krzysztof

    2005-10-01

    Six-dimensional intermolecular potential energy surfaces (PESs) for the interaction of CH4 with H2O are presented, obtained from ab initio calculations using symmetry-adapted perturbation theory (SAPT) at two different levels of intramonomer correlation and the supermolecular approach at three different levels of electron correlation. Both CH4 and H2O are assumed to be rigid molecules with interatomic distances and angles fixed at the average values in the ground-state vibration. A physically motivated analytical expression for each PES has been developed as a sum of site-site functions. The PES of the CH4-H2O dimer has only two symmetry-distinct minima. From the SAPT calculations, the global minimum has an energy of -1.03kcal /mol at a geometry where H2O is the proton donor, HO -H⋯CH4, with the O-H-C angle of 165°, while the secondary minimum, with an energy of -0.72kcal/mol, has CH4 in the role of the proton donor (H3C -H⋯OH2). We estimated the complete basis set limit of the SAPT interaction energy at the global minimum to be -1.06kcal/mol. The classical cross second virial coefficient B12(T) has been calculated for the temperature range 298-653K. Our best results agree well with some experiments, allowing an evaluation of the quality of experimental results.

  5. Potential energy curves via double electron-attachment calculations: Dissociation of alkali metal dimers

    NASA Astrophysics Data System (ADS)

    Musiał, Monika; Kowalska-Szojda, Katarzyna; Lyakh, Dmitry I.; Bartlett, Rodney J.

    2013-05-01

    The recently developed method [M. Musiał, J. Chem. Phys. 136, 134111 (2012), 10.1063/1.3700438] to study double electron attached states has been applied to the description of the ground and excited state potential energy curves of the alkali metal dimers. The method is based on the multireference coupled cluster scheme formulated within the Fock space formalism for the (2,0) sector. Due to the use of the efficient intermediate Hamiltonian formulation, the approach is free from the intruder states problem. The description of the neutral alkali metal dimers is accomplished via attaching two electrons to the corresponding doubly ionized system. This way is particularly advantageous when a closed shell molecule dissociates into open shell subunits while its doubly positive cation generates the closed shell fragments. In the current work, we generate the potential energy curves for the ground and multiple excited states of the Li2 and Na2 molecules. In all cases the potential energy curves are smooth for the entire range of interatomic distances (from the equilibrium point to the dissociation limit). Based on the calculated potential energy curves, we are able to compute spectroscopic parameters of the systems studied.

  6. The ground-state potential energy curve of the radium dimer from relativistic coupled cluster calculations.

    PubMed

    Teodoro, Tiago Quevedo; Haiduke, Roberto Luiz Andrade; Dammalapati, Umakanth; Knoop, Steven; Visscher, Lucas

    2015-08-28

    The potential energy curve for the ground-state of radium dimer (Ra2) is provided by means of atomic and molecular relativistic coupled cluster calculations. The short-range part of this curve is defined by an equilibrium bond length of 5.324 Å, a dissociation energy of 897 cm(-1), and a harmonic vibrational frequency of 20.5 cm(-1). The asymptotic behavior at large interatomic distances is characterized by the van der Waals coefficients C6 = 5.090 × 10(3), C8 = 6.978 × 10(5), and C10 = 8.786 × 10(7) atomic units. The two regions are matched in an analytical potential to provide a convenient representation for use in further calculations, for instance, to model cold collisions between radium atoms. This might become relevant in future experiments on ultracold, optically trapped, radioactive radium atoms that are used to search for a permanent electric dipole moment. PMID:26328843

  7. Multireference configuration interaction calculations of the first six ionization potentials of the uranium atom

    SciTech Connect

    Bross, David H.; Parmar, Payal; Peterson, Kirk A.

    2015-11-12

    The first 6 ionization potentials (IPs) of the uranium atom have been calculated using multireference configuration interaction (MRCI+Q) with extrapolations to the complete basis set (CBS) limit using new all-electron correlation consistent basis sets. The latter were carried out with the third-order Douglas-Kroll-Hess Hamiltonian. Correlation down through the 5s5p5d electrons have been taken into account, as well as contributions to the IPs due to the Lamb shift. Spin-orbit coupling contributions calculated at the 4-component Kramers restricted configuration interaction level, as well as the Gaunt term computed at the Dirac-Hartree-Fock level, were added to the best scalar relativistic results. As a result, the final ionization potentials are expected to be accurate to at least 5 kcal/mol (0.2 eV), and thus more reliable than the current experimental values of IP3 through IP6.

  8. Calculation of free energies and chemical potentials for gas hydrates using Monte Carlo simulations.

    PubMed

    Wierzchowski, S J; Monson, P A

    2007-06-28

    We describe a method for calculating free energies and chemical potentials for molecular models of gas hydrate systems using Monte Carlo simulations. The method has two components: (i) thermodynamic integration to obtain the water and guest molecule chemical potentials as functions of the hydrate occupancy; (ii) calculation of the free energy of the zero-occupancy hydrate system using thermodynamic integration from an Einstein crystal reference state. The approach is applicable to any classical molecular model of a hydrate. We illustrate the methodology with an application to the structure-I methane hydrate using two molecular models. Results from the method are also used to assess approximations in the van der Waals-Platteeuw theory and some of its extensions. It is shown that the success of the van der Waals-Platteeuw theory is in part due to a cancellation of the error arising from the assumption of a fixed configuration of water molecules in the hydrate framework with that arising from the neglect of methane-methane interactions. PMID:17530795

  9. Multireference configuration interaction calculations of the first six ionization potentials of the uranium atom

    DOE PAGESBeta

    Bross, David H.; Parmar, Payal; Peterson, Kirk A.

    2015-11-12

    The first 6 ionization potentials (IPs) of the uranium atom have been calculated using multireference configuration interaction (MRCI+Q) with extrapolations to the complete basis set (CBS) limit using new all-electron correlation consistent basis sets. The latter were carried out with the third-order Douglas-Kroll-Hess Hamiltonian. Correlation down through the 5s5p5d electrons have been taken into account, as well as contributions to the IPs due to the Lamb shift. Spin-orbit coupling contributions calculated at the 4-component Kramers restricted configuration interaction level, as well as the Gaunt term computed at the Dirac-Hartree-Fock level, were added to the best scalar relativistic results. As amore » result, the final ionization potentials are expected to be accurate to at least 5 kcal/mol (0.2 eV), and thus more reliable than the current experimental values of IP3 through IP6.« less

  10. Calculation of turbulence effects in an upward-refracting atmosphere

    NASA Astrophysics Data System (ADS)

    Gilbert, Kenneth E.; di, Xiao; Raspet, Richard

    1990-06-01

    The effect of atmospheric turbulence on sound propagation was investigated for both nonrefractive and refractive atmospheres, using the parabolic equation method of Gilbert and White (1989) in conjunction with a two-dimensional atmospheric turbulence model. The calculations for a nonrefractive atmosphere gave good agreement with experimental data and with Daigle's (1979) theory, while calculations for an upward-refractive atmosphere gave reasonable agreement with the data of Weiner and Keast (1959). It is concluded that, for a receiver deep in a shadow zone and for frequencies greater than a few hundred hertz, the measured sound-pressure level is due almost entirely to the sound scattered into the shadow zone by atmospheric turbulence. Consequently, for upward refraction and frequencies above a few hundred hertz, turbulence must be included in long-range propagation calculations.

  11. Seasonal variation of redox species and redox potentials in shallow groundwater: A comparison of measured and calculated redox potentials

    NASA Astrophysics Data System (ADS)

    Ramesh Kumar, A.; Riyazuddin, P.

    2012-06-01

    SummaryThe seasonal variation of redox potential (Eh) and redox species such as As(V)/As(III), Cr(VI)/Cr(III), Fe(III)/Fe(II), NO3-/NO2-, and Se(VI)/Se(IV) were studied in a shallow groundwater for a period of three years (May, 2004-January, 2007). The study area was Chrompet area of Chennai city, India. Groundwater samples from 65 wells were monitored for pH, electrical conductivity, dissolved oxygen (DO), and major ions during pre-(May) and post-monsoon (January) seasons. The objective of the study was to gain insight into the temporal variation of the redox species due to groundwater recharge and to identify the redox reactions controlling the measured Eh of the groundwater. The study revealed that the shallow groundwater was "oxic" with DO ranging between 0.25 and 5.00 mg L-1, and between 0.38 and 5.05 mg L-1 during pre-(May, 2004) and post-monsoon (January, 2005) seasons, respectively. The measured Eh (with respect to standard hydrogen electrode, SHE) ranged between 65 and 322 mV, and between 110 and 330 mV during pre- and post-monsoon seasons, respectively. During post-monsoon seasons, DO and Eh increased in most of the wells due to groundwater recharge. The calculated Eh using the redox couples As(V)/As(III), NO3-/NO2-, O2/H2O and Se(VI)/Se(IV) neither agreed among themselves nor with the measured Eh during all the seasons. It shows that in the shallow groundwater, the various redox couples are in disequilibrium among themselves and with the Pt electrode. However, 41% (n = 122) of the Eh values calculated from Fe(III)/Fe(II) couple agreed with the measured Eh within ±30 mV, the uncertainty of Pt-electrode measurement. The post-monsoon seasons showed higher values of As(V)/As(III) and Se(VI)/Se(IV) compared to the pre-monsoon seasons, whereas Fe(III)/Fe(II) behaved in the opposite manner. This pattern of variation is consistent with the increased oxidizing nature, as shown by the higher DO and Eh values observed during post-monsoon seasons. The results

  12. Hybrid MP2/MP4 potential surfaces in VSCF calculations of IR spectra: applications for organic molecules.

    PubMed

    Knaanie, Roie; Šebek, Jiří; Kalinowski, Jaroslaw; Benny Gerber, R

    2014-02-01

    This study introduces an improved hybrid MP2/MP4 ab initio potential for vibrational spectroscopy calculations which is very accurate, yet without high computational demands. The method uses harmonic vibrational calculations with the MP4(SDQ) potential to construct an improved MP2 potential by coordinate scaling. This improved MP2 potential is used for the anharmonic VSCF calculation. The method was tested spectroscopically for four molecules: butane, acetone, ethylene and glycine. Very good agreement with experiment was found. For most of the systems, the more accurate harmonic treatment considerably improved the MP2 anharmonic results. PMID:23838574

  13. Effective Connectivity Reveals Strategy Differences in an Expert Calculator

    PubMed Central

    Minati, Ludovico; Sigala, Natasha

    2013-01-01

    Mathematical reasoning is a core component of cognition and the study of experts defines the upper limits of human cognitive abilities, which is why we are fascinated by peak performers, such as chess masters and mental calculators. Here, we investigated the neural bases of calendrical skills, i.e. the ability to rapidly identify the weekday of a particular date, in a gifted mental calculator who does not fall in the autistic spectrum, using functional MRI. Graph-based mapping of effective connectivity, but not univariate analysis, revealed distinct anatomical location of “cortical hubs” supporting the processing of well-practiced close dates and less-practiced remote dates: the former engaged predominantly occipital and medial temporal areas, whereas the latter were associated mainly with prefrontal, orbitofrontal and anterior cingulate connectivity. These results point to the effect of extensive practice on the development of expertise and long term working memory, and demonstrate the role of frontal networks in supporting performance on less practiced calculations, which incur additional processing demands. Through the example of calendrical skills, our results demonstrate that the ability to perform complex calculations is initially supported by extensive attentional and strategic resources, which, as expertise develops, are gradually replaced by access to long term working memory for familiar material. PMID:24086291

  14. Absorbed Dose and Dose Equivalent Calculations for Modeling Effective Dose

    NASA Technical Reports Server (NTRS)

    Welton, Andrew; Lee, Kerry

    2010-01-01

    While in orbit, Astronauts are exposed to a much higher dose of ionizing radiation than when on the ground. It is important to model how shielding designs on spacecraft reduce radiation effective dose pre-flight, and determine whether or not a danger to humans is presented. However, in order to calculate effective dose, dose equivalent calculations are needed. Dose equivalent takes into account an absorbed dose of radiation and the biological effectiveness of ionizing radiation. This is important in preventing long-term, stochastic radiation effects in humans spending time in space. Monte carlo simulations run with the particle transport code FLUKA, give absorbed and equivalent dose data for relevant shielding. The shielding geometry used in the dose calculations is a layered slab design, consisting of aluminum, polyethylene, and water. Water is used to simulate the soft tissues that compose the human body. The results obtained will provide information on how the shielding performs with many thicknesses of each material in the slab. This allows them to be directly applicable to modern spacecraft shielding geometries.

  15. Infrared lens thermal effect: equivalent focal shift and calculating model

    NASA Astrophysics Data System (ADS)

    Zhang, Cheng-shuo; Shi, Zelin; Feng, Bin; Xu, Bao-shu

    2014-11-01

    It's well-know that the focal shift of infrared lens is the major factor in degeneration of imaging quality when temperature change. In order to figure out the connection between temperature change and focal shift, partial differential equations of thermal effect on light path are obtained by raytrace method, to begin with. The approximately solution of the PDEs show that focal shift is proportional to temperature change. And a formula to compute the proportional factor is given. In order to understand infrared lens thermal effect deeply, we use defocus by image plane shift at constant temperature to equivalently represent thermal effect on infrared lens. So equivalent focal shift (EFS) is defined and its calculating model is proposed at last. In order to verify EFS and its calculating model, Physical experimental platform including a motorized linear stage with built-in controller, blackbody, target, collimator, IR detector, computer and other devices is developed. The experimental results indicate that EFS make the image plane shift at constant temperature have the same influence on infrared lens as thermal effect and its calculating model is correct.

  16. Calculation of the Lightning Potential Index and electric field in numerical weather prediction models

    NASA Astrophysics Data System (ADS)

    Dementyeva, S. O.; Ilin, N. V.; Mareev, E. A.

    2015-03-01

    Modern methods for predicting thunderstorms and lightnings with the use of high-resolution numerical models are considered. An analysis of the Lightning Potential Index (LPI) is performed for various microphysics parameterizations with the use of the Weather Research and Forecasting (WRF) model. The maximum index values are shown to depend significantly on the type of parameterization. This makes it impossible to specify a single threshold LPI for various parameterizations as a criterion for the occurrence of lightning flashes. The topographic LPI maps underestimate the sizes of regions of likely thunderstorm-hazard events. Calculating the electric field under the assumption that ice and graupel are the main charge carriers is considered a new algorithm of lightning prediction. The model shows that the potential difference (between the ground and cloud layer at a given altitude) sufficient to generate a discharge is retained in a larger region than is predicted by the LPI. The main features of the spatial distribution of the electric field and potential agree with observed data.

  17. Hybrid Density Functional Calculations of Redox Potentials of Transition Metal Compounds

    NASA Astrophysics Data System (ADS)

    Armiento, Rickard; Chevrier, Vincent; Ong, Shyue Ping; Ceder, Gerbrand

    2010-03-01

    Prior works have shown that density functional theory (DFT) with the DFT+U method resolves the underestimation of redox potentials calculated by conventional functionals for a number of transition metal compounds relevant for battery applications, including the olivine LixMPO4 (M = Fe, Mn, Co, Ni), layered LixMO2 (M = Co, Ni) and spinel-like LixMn2O4. We show that the redox potentials of these compounds are also well reproduced by the hybrid density functional by Heyd-Scuseria-Ernzerhof (HSE06). Hybrid functionals combine a conventional DFT functional with a part of Hartree-Fock (HF) exchange. While the HF part increases the computational expense by at least one order of magnitude, it provides, in contrast to DFT+U, a correction for the self-interaction error that does not rely on special treatment of the occupancies of the orbital states of ions or species-specific parameters. We compare the accuracy of regular DFT, DFT+U and HSE06 for the redox potentials, lattice constants, and other properties. Examples of electron delocalization problems connected to the self-interaction error in the systems are discussed, and shown to be resolved both by the hybrid functional and DFT+U methods. Comments are made on the possibility to approach the delocalization problem with a semi-local functional.

  18. Comparative Study of Defect Properties in GaN: Ab initio and Empirical Potential Calculations

    SciTech Connect

    Gao, Fei; Bylaska, Eric J.; El-Azab, Anter; Weber, William J.; LM Wang, R Fromknecht, LL Snead, DF Downey and H Takahashi

    2004-04-05

    Density functional theory (DFT) is used to study the formation, properties and atomic configurations of monovacancies, antisite defects and possible interstitials in GaN. The relaxation around a vacancy is generally small, but the relaxation around antisite defects is large, particularly for a Ga antisite defect, which is not stable and converts to an N-N<0001> split interstitial. All N interstitials, starting from any possible sites, eventually transfer into the N-N split interstitials, forming N molecules. In the case of Ga interstitials, the most favorable configuration is the Ga octahedral interstitial. However, it is found that the Ga-Ga<> split interstitial can bridge the gap between non-bonded Ga atoms along the <11-2> direction, which leads to the formation of Ga atomic wires in GaN, with bond distance close to those noted in bulk Ga. In addition, two representative potentials, namely Stillinger-Weber and Tersoff-Brenner potentials, have been employed to deter mine the formation of defects using molecular dynamics (MD) method in GaN. The MD results are discussed and compared to DFT calculations. The present DFT and MD results provide guidelines for evaluating the quality and fit of empirical potentials for large-scale simulations of ion-solid interaction and thermal annealing of defects in GaN.

  19. Calculation of Steady-state Evaporation for an Arbitrary Matrix Potential at Ground Surface

    NASA Astrophysics Data System (ADS)

    Liu, X.; Zhan, H.

    2014-12-01

    The water loss from soil by evaporation and the amount of ground water available to plants due to the upward movement of water from a water table is an important topic in many disciplines such as soil science, hydrology, and plant physiology. Although water evaporation in actual field setting is a highly complex process, a nearly steady upward flow from a water table to a bare soil surface may be established if the daily evaporative demand is reasonably uniform for a long period of time. While the maximum potential rate of evaporation from the ground surface depends on atmospheric conditions, the actual flux across the soil surface is limited by the ability of the porous medium for transmitting water from the unsaturated zone.The purpose of this study is to calculate the steady-state evaporation for an arbitrary matrix potential at bare soil surface above a shallow water table, while the unsaturated hydraulic conductivity is a nonlinear function of water content or matrix potential. The Haverkamp function and the Brooks-Corey function for the unsaturated hydraulic conductivity are used, and the study results are contrast among the solution developed from the two retention equation and HYDRUS simulation.

  20. A new Fe–He interatomic potential based on ab initio calculations in α-Fe

    SciTech Connect

    Gao, Fei; Deng, Huiqiu; Heinisch, Howard L.; Kurtz, Richard J.

    2011-11-01

    A new interatomic potential for Fe−He interactions has been developed by fitting to the results obtained from ab initio calculations. Based on the electronic hybridization between Fe d-electrons and He s-electrons, an s-band model, along with a repulsive pair potential, has been developed to describe the Fe−He interaction. The atomic configurations and formation energies of single He defects and small interstitial He clusters are utilized in the fitting process. The binding properties and relative stabilities of the He-vacancy and interstitial He clusters are studied. The present Fe−He potential is also applied to study the emission of self-interstitial atoms from small He clusters in alpha-Fe matrices. It is found that the di-He cluster dissociates when the temperature is higher than 400K, but the larger He clusters can create an interstitial Fe atom. The temperature for kicking out an interstitial Fe atom is found to decrease with increasing size of the He clusters.

  1. Chemically transferable coarse-grained potentials from conditional reversible work calculations.

    PubMed

    Brini, E; van der Vegt, N F A

    2012-10-21

    The representability and transferability of effective pair potentials used in multiscale simulations of soft matter systems is ill understood. In this paper, we study liquid state systems composed of n-alkanes, the coarse-grained (CG) potential of which may be assumed pairwise additive and has been obtained using the conditional reversible work (CRW) method. The CRW method is a free-energy-based coarse-graining procedure, which, by means of performing the coarse graining at pair level, rigorously provides a pair potential that describes the interaction free energy between two mapped atom groups (beads) embedded in their respective chemical environments. The pairwise nature of the interactions combined with their dependence on the chemically bonded environment makes CRW potentials ideally suited in studies of chemical transferability. We report CRW potentials for hexane using a mapping scheme that merges two heavy atoms in one CG bead. It is shown that the model is chemically and thermodynamically transferable to alkanes of different chain lengths in the liquid phase at temperatures between the melting and the boiling point under atmospheric (1 atm) pressure conditions. It is further shown that CRW-CG potentials may be readily obtained from a single simulation of the liquid state using the free energy perturbation method, thereby providing a fast and versatile molecular coarse graining method for aliphatic molecules. PMID:23083154

  2. Coaxial electromagnetic launcher calculations using FE-BE method and hybrid potentials

    NASA Astrophysics Data System (ADS)

    Chang, J. H.; Becker, E. B.; Driga, M. D.

    1993-01-01

    In this paper, a hybrid method combining finite and boundary elements (the FE-BE method) is presented to analyze the transient electromagnetic/mechanical behavior of coaxial induction launchers (i.e., the coilgun). The corresponding initial/boundary value problem is formulated in terms of the hybrid potentials, which mixes the vector and scalar magnetic potential functions. The problem is assumed to be axisymmetric and the forcing currents are circumferential. Thermal effect due to ohmic losses is considered during the launching processes.

  3. Consistent use of the standard model effective potential.

    PubMed

    Andreassen, Anders; Frost, William; Schwartz, Matthew D

    2014-12-12

    The stability of the standard model is determined by the true minimum of the effective Higgs potential. We show that the potential at its minimum when computed by the traditional method is strongly dependent on the gauge parameter. It moreover depends on the scale where the potential is calculated. We provide a consistent method for determining absolute stability independent of both gauge and calculation scale, order by order in perturbation theory. This leads to a revised stability bounds m(h)(pole)>(129.4±2.3)  GeV and m(t)(pole)<(171.2±0.3)  GeV. We also show how to evaluate the effect of new physics on the stability bound without resorting to unphysical field values. PMID:25541764

  4. Reduced neutron spectroscopic factors when using potential geometries constrained by Hartree-Fock calculations

    SciTech Connect

    Lee, Jenny; Brown, B.A.; Delaunay, F.; Lynch, W.G.; Saelim, M. J.; Tsang, M.B.; Tostevin, J.A.

    2006-04-15

    We carry out a systematic analysis of angular distribution measurements for selected ground-state to ground-state (d,p) and (p,d) neutron transfer reactions, including the calcium isotopes. We propose a consistent three-body model reaction methodology in which we constrain the transferred-neutron bound state and nucleon-target optical potential geometries using modern Hartree-Fock calculations. Our deduced neutron spectroscopic factors are found to be suppressed by {approx}30% relative to independent-particle shell-model values, from {sup 40}Ca through {sup 49}Ca. The other nuclei studied, ranging from B to Ti, show similar average suppressions with respect to large-basis shell-model expectations. Our results are consistent with deduced spectroscopic strengths for neutrons and protons from intermediate-energy nucleon knockout reactions and for protons from (e,e{sup '}p) reactions on well-bound nuclei.

  5. Spectroscopic studies, potential energy surface and molecular orbital calculations of pramipexole.

    PubMed

    Muthu, S; Uma Maheswari, J; Srinivasan, S; Isac paulraj, E

    2013-11-01

    A systematic vibrational spectroscopic assignment and analysis of pramipexole [(S)-N(6)-propyl-4,5,6,7-tetrahydro-1,3-benzothiazole-2,6-diamine] has been carried out using FT-IR and FT-Raman spectral data. The vibrational analysis was aided by an electronic structure calculation based on the hybrid density functional method B3LYP using a 6-311G(d, p) and cc-pVTZ basis sets. Molecular equilibrium geometries, electronic energies, IR and Raman intensities, harmonic vibrational frequencies have been computed. The assignments are based on the experimental IR and Raman spectra, and a complete assignment of the observed spectra has been proposed. The UV-visible spectrum of the compound was recorded and the electronic properties, such as HOMO and LUMO energies and the maximum absorption λmax were determined by time-dependent DFT (TD-DFT) method. The geometrical parameters, vibrational frequencies and absorption wavelengths were compared with the experimental data. The complete vibrational assignments are performed on the basis of the potential energy distributions (PEDs) of the vibrational modes in terms of natural internal coordinates. The simulated FT-IR, FT-Raman, and UV spectra of the title compound have been constructed. In addition, the potential energy surface, HOMO and LUMO energies, the molecular electrostatic potential and the first-order hyperpolarizability have been computed. The magnitude of the first-order hyperpolarizability is 5 times larger than that of urea and the title compound may be a potential applicant for the development of NLO materials. PMID:23831980

  6. Full potential calculation of electronics and thermoelectric properties of doped Mg{sub 2}Si

    SciTech Connect

    Poopanya, P.; Yangthaisong, A.

    2013-12-04

    We present the calculations of the electronic structure and transport properties on the anti-fluorite Mg{sub 2}Si using the full potential linearized augmented plane-wave (FP-LAPW) method and the semi-classical Boltzmann theory. The modified Becke-Johnson (mBJ) exchange potentials are used to derive energy gaps and correct band gaps according to experimental values. It is found that Mg{sub 2}Si is an indirect band gap (Γ→X) material with the gap of 0.56 eV which is in good agreement with the experimental observation. Note that the band structure of Mg{sub 2}Si is directly used in combination with the semi-classical Boltzmann theory to obtain the transport coefficients. It is found that the material is the n-type semiconductor with the lowest electron concentration of 3.03×10{sup 14} cm{sup −3} at 300 K. We have also calculated the thermoelectric properties of Mg{sub 2}Si based on the rigid band approximation by varying the p-type and n-type doping levels. At room temperature, the highest power factor for p-type and n-type dopants are obtained at the hole and electron concentration of 1.63×10{sup 20} cm{sup −3} and 1.15×1021 cm{sup −3}, respectively. From the electronic states, we also found that the n-type doping region is dominated by the Mg−2p{sup 6} 3s{sup 2} and Si−3p{sup 2} states, while the Mg−2p{sup 6} and Si−3p{sup 2} states are important in the p-type doped Mg{sub 2}Si.

  7. Potential Health Effects from Groundwater Pollution.

    ERIC Educational Resources Information Center

    Goyer, Robert A.

    1985-01-01

    Discusses the growing awareness of potential toxicological effects of synthetic organic chemicals contaminating groundwater. Problems concerning pesticides, chlorination, epidemiologic studies, cancer, nephrotoxicity, and considerations of risk are addressed. Additional research in this area is advocated. (DH)

  8. Calculation of the transport properties of carbon dioxide. II. Thermal conductivity and thermomagnetic effects.

    PubMed

    Bock, Steffen; Bich, Eckard; Vogel, Eckhard; Dickinson, Alan S; Vesovic, Velisa

    2004-05-01

    The transport properties of pure carbon dioxide have been calculated from the intermolecular potential using the classical trajectory method. Results are reported in the dilute-gas limit for thermal conductivity and thermomagnetic coefficients for temperatures ranging from 200 K to 1000 K. Three recent carbon dioxide potential energy hypersurfaces have been investigated. Since thermal conductivity is influenced by vibrational degrees of freedom, not included in the rigid-rotor classical trajectory calculation, a correction for vibration has also been employed. The calculations indicate that the second-order thermal conductivity corrections due to the angular momentum polarization (< 2%) and velocity polarization (< 1%) are both small. Thermal conductivity values calculated using the potential energy hypersurface by Bukowski et al. (1999) are in good agreement with the available experimental data. They underestimate the best experimental data at room temperature by 1% and in the range up to 470 K by 1%-3%, depending on the data source. Outside this range the calculated values, we believe, may be more reliable than the currently available experimental data. Our results are consistent with measurements of the thermomagnetic effect at 300 K only when the vibrational degrees of freedom are considered fully. This excellent agreement for these properties indicates that particularly the potential surface of Bukowski et al. provides a realistic description of the anisotropy of the surface. PMID:15267716

  9. High speed inlet calculations with real gas effects

    NASA Technical Reports Server (NTRS)

    Coirier, William J.

    1988-01-01

    A 2-D steady-state Navier-Stokes solver has been upgraded to include the effects of frozen and equilibrium air chemistry for applications to high speed flight vehicles. To provide a computationally economical first order approximation to the high temperature physics, variable thermodynamic data is used for the chemically frozen mode to allow for a variation with temperature of the air specific heats and enthalpy. For calculations involving air in chemical equilibrium, a specially modified version of the NASA Lewis Chemical Equilibrium Code, CEC, is used to compute the chemical composition and resultant thermochemical properties. The upgraded solver is demonstrated by comparing results from calorically perfect (C sub p=constant), thermally perfect (frozen) and equilibrium air calculations for a variety of geometries, and flight Mach numbers.

  10. Ab initio molecular orbital calculations of DNA radical ions. 5. Scaling of calculated electron affinities and ionization potentials to experimental values

    SciTech Connect

    Sevilla, M.D.; Colson, A.O. ); Besler, B. )

    1995-01-19

    Ab initio molecular orbital calculations of the electron affinities (EAs) and ionization potentials (IPs) of the DNA bases are presented in this work. Comparisons of calculated and experimental values are made for a series of compounds of size and/or structure similar to the DNA bases. Excellent correlations between calculated and experimental values are found for both Koopmans EAs at the 6-31G[sup *] and D95v levels and calculated vertical EAs of the model compounds. Several basis sets are considered: 6-31G[sup *], 6-31+G(d), and D95v. Calculations at 6-31G[sup *] and 6-31+G(d) using both ROHF and ROMP2 theories show a consistent difference between calculated vertical and adiabatic EAs. This allows for a good estimate of DNA base adiabatic EAs. i.e., -0.7, -0.3, 0.2, 0.3, and 0.4 eV; from the vertical EAs -1.23, -0.74, -0.40, -0.32, and -0.19 eV for G, A, C, T, and U respectively. While EAs must be scaled, we find that Koopmans IPs calculated at the simple 3-21G level predict vertical IPs of the DNA bases with only a 0.15 eV average absolute deviation from the experimentally reported values and calculations at MP2/6-31+G(d)//6-31G[sup *] for the adiabatic ionization potentials of the DNA bases are all within 0.1 eV of experiment. 41 refs., 2 figs., 5 tabs.

  11. Perturbative Calculation of Quasi-Potential in Non-equilibrium Diffusions: A Mean-Field Example

    NASA Astrophysics Data System (ADS)

    Bouchet, Freddy; Gawȩdzki, Krzysztof; Nardini, Cesare

    2016-04-01

    In stochastic systems with weak noise, the logarithm of the stationary distribution becomes proportional to a large deviation rate function called the quasi-potential. The quasi-potential, and its characterization through a variational problem, lies at the core of the Freidlin-Wentzell large deviations theory (Freidlin and Wentzell, Random perturbations of dynamical systems, 2012). In many interacting particle systems, the particle density is described by fluctuating hydrodynamics governed by Macroscopic Fluctuation Theory (Bertini et al., arXiv:1404.6466, 2014), which formally fits within Freidlin-Wentzell's framework with a weak noise proportional to 1/√{N} , where N is the number of particles. The quasi-potential then appears as a natural generalization of the equilibrium free energy to non-equilibrium particle systems. A key physical and practical issue is to actually compute quasi-potentials from their variational characterization for non-equilibrium systems for which detailed balance does not hold. We discuss how to perform such a computation perturbatively in an external parameter λ , starting from a known quasi-potential for λ =0 . In a general setup, explicit iterative formulae for all terms of the power-series expansion of the quasi-potential are given for the first time. The key point is a proof of solvability conditions that assure the existence of the perturbation expansion to all orders. We apply the perturbative approach to diffusive particles interacting through a mean-field potential. For such systems, the variational characterization of the quasi-potential was proven by Dawson and Gartner (Stochastics 20:247-308, 1987; Stochastic differential systems, vol 96, pp 1-10, 1987). Our perturbative analysis provides new explicit results about the quasi-potential and about fluctuations of one-particle observables in a simple example of mean field diffusions: the Shinomoto-Kuramoto model of coupled rotators (Prog Theoret Phys 75:1105-1110, [74]). This

  12. Perturbative Calculation of Quasi-Potential in Non-equilibrium Diffusions: A Mean-Field Example

    NASA Astrophysics Data System (ADS)

    Bouchet, Freddy; Gawȩdzki, Krzysztof; Nardini, Cesare

    2016-06-01

    In stochastic systems with weak noise, the logarithm of the stationary distribution becomes proportional to a large deviation rate function called the quasi-potential. The quasi-potential, and its characterization through a variational problem, lies at the core of the Freidlin-Wentzell large deviations theory (Freidlin and Wentzell, Random perturbations of dynamical systems, 2012). In many interacting particle systems, the particle density is described by fluctuating hydrodynamics governed by Macroscopic Fluctuation Theory (Bertini et al., arXiv:1404.6466 , 2014), which formally fits within Freidlin-Wentzell's framework with a weak noise proportional to 1/√{N}, where N is the number of particles. The quasi-potential then appears as a natural generalization of the equilibrium free energy to non-equilibrium particle systems. A key physical and practical issue is to actually compute quasi-potentials from their variational characterization for non-equilibrium systems for which detailed balance does not hold. We discuss how to perform such a computation perturbatively in an external parameter λ , starting from a known quasi-potential for λ =0. In a general setup, explicit iterative formulae for all terms of the power-series expansion of the quasi-potential are given for the first time. The key point is a proof of solvability conditions that assure the existence of the perturbation expansion to all orders. We apply the perturbative approach to diffusive particles interacting through a mean-field potential. For such systems, the variational characterization of the quasi-potential was proven by Dawson and Gartner (Stochastics 20:247-308, 1987; Stochastic differential systems, vol 96, pp 1-10, 1987). Our perturbative analysis provides new explicit results about the quasi-potential and about fluctuations of one-particle observables in a simple example

  13. Radion effective potential in brane gas cosmology

    SciTech Connect

    Kim, Jin Young

    2008-09-15

    We consider a cosmological solution which can explain anisotropic evolution of spatial dimensions and the stabilization of extra dimensions in brane gas formalism. We evaluate the effective potentials, induced by brane gas, bulk flux and supergravity particles, which govern the sizes of the observed three and the extra dimensions. It is possible that the wrapped internal volume can oscillate between two turning points or sit at the minimum of the potential while the unwrapped three-dimensional volume can expand monotonically. Including the supergravity particles makes the effective potential steeper as the internal volume shrinks.

  14. A comparison of estimated and calculated effective porosity

    NASA Astrophysics Data System (ADS)

    Stephens, Daniel B.; Hsu, Kuo-Chin; Prieksat, Mark A.; Ankeny, Mark D.; Blandford, Neil; Roth, Tracy L.; Kelsey, James A.; Whitworth, Julia R.

    Effective porosity in solute-transport analyses is usually estimated rather than calculated from tracer tests in the field or laboratory. Calculated values of effective porosity in the laboratory on three different textured samples were compared to estimates derived from particle-size distributions and soil-water characteristic curves. The agreement was poor and it seems that no clear relationships exist between effective porosity calculated from laboratory tracer tests and effective porosity estimated from particle-size distributions and soil-water characteristic curves. A field tracer test in a sand-and-gravel aquifer produced a calculated effective porosity of approximately 0.17. By comparison, estimates of effective porosity from textural data, moisture retention, and published values were approximately 50-90% greater than the field calibrated value. Thus, estimation of effective porosity for chemical transport is highly dependent on the chosen transport model and is best obtained by laboratory or field tracer tests. Résumé La porosité effective dans les analyses de transport de soluté est habituellement estimée, plutôt que calculée à partir d'expériences de traçage sur le terrain ou au laboratoire. Les valeurs calculées de la porosité effective au laboratoire sur trois échantillons de textures différentes ont été comparées aux estimations provenant de distributions de taille de particules et de courbes caractéristiques sol-eau. La concordance était plutôt faible et il semble qu'il n'existe aucune relation claire entre la porosité effective calculée à partir des expériences de traçage au laboratoire et la porosité effective estimée à partir des distributions de taille de particules et de courbes caractéristiques sol-eau. Une expérience de traçage de terrain dans un aquifère de sables et de graviers a fourni une porosité effective calculée d'environ 0,17. En comparaison, les estimations de porosité effective de données de

  15. Electron Affinity Calculations for Atoms: Sensitive Probe of Many-Body Effects

    NASA Astrophysics Data System (ADS)

    Felfli, Z.; Msezane, A. Z.

    2016-05-01

    Electron-electron correlations and core-polarization interactions are crucial for the existence and stability of most negative ions. Therefore, they can be used as a sensitive probe of many-body effects in the calculation of the electron affinities (EAs) of atoms. The importance of relativistic effects in the calculation of the EAs of atoms has recently been assessed to be insignificant up to Z of 85. Here we use the complex angular momentum (CAM) methodology wherein is embedded fully the electron-electron correlations, to investigate core-polarization interactions in low-energy electron elastic scattering from the atoms In, Sn, Eu, Au and At through the calculation of their EAs. For the core-polarization interaction we use the rational function approximation of the Thomas-Fermi potential, which can be analytically continued into the complex plane. The EAs are extracted from the large resonance peaks in the CAM calculated low-energy electron-atom scattering total cross sections and compared with those from measurements and sophisticated theoretical methods. It is concluded that when the electron-electron correlations and core polarization interactions (both major many-body effects) are accounted for adequately the importance of relativity on the calculation of the EAs of atoms can be assessed. Even for the high Z (85) At atom relativistic effects are estimated to contribute a maximum of 3.6% to its EA calculation.

  16. Accelerating quantum instanton calculations of the kinetic isotope effects.

    PubMed

    Karandashev, Konstantin; Vaníček, Jiří

    2015-11-21

    Path integral implementation of the quantum instanton approximation currently belongs among the most accurate methods for computing quantum rate constants and kinetic isotope effects, but its use has been limited due to the rather high computational cost. Here, we demonstrate that the efficiency of quantum instanton calculations of the kinetic isotope effects can be increased by orders of magnitude by combining two approaches: The convergence to the quantum limit is accelerated by employing high-order path integral factorizations of the Boltzmann operator, while the statistical convergence is improved by implementing virial estimators for relevant quantities. After deriving several new virial estimators for the high-order factorization and evaluating the resulting increase in efficiency, using ⋅Hα + HβHγ → HαHβ + ⋅ Hγ reaction as an example, we apply the proposed method to obtain several kinetic isotope effects on CH4 + ⋅ H ⇌ ⋅ CH3 + H2 forward and backward reactions. PMID:26590524

  17. Dynamics study of the OH + NH3 hydrogen abstraction reaction using QCT calculations based on an analytical potential energy surface

    NASA Astrophysics Data System (ADS)

    Monge-Palacios, M.; Corchado, J. C.; Espinosa-Garcia, J.

    2013-06-01

    To understand the reactivity and mechanism of the OH + NH3 → H2O + NH2 gas-phase reaction, which evolves through wells in the entrance and exit channels, a detailed dynamics study was carried out using quasi-classical trajectory calculations. The calculations were performed on an analytical potential energy surface (PES) recently developed by our group, PES-2012 [Monge-Palacios et al. J. Chem. Phys. 138, 084305 (2013)], 10.1063/1.4792719. Most of the available energy appeared as H2O product vibrational energy (54%), reproducing the only experimental evidence, while only the 21% of this energy appeared as NH2 co-product vibrational energy. Both products appeared with cold and broad rotational distributions. The excitation function (constant collision energy in the range 1.0-14.0 kcal mol-1) increases smoothly with energy, contrasting with the only theoretical information (reduced-dimensional quantum scattering calculations based on a simplified PES), which presented a peak at low collision energies, related to quantized states. Analysis of the individual reactive trajectories showed that different mechanisms operate depending on the collision energy. Thus, while at high energies (Ecoll ≥ 6 kcal mol-1) all trajectories are direct, at low energies about 20%-30% of trajectories are indirect, i.e., with the mediation of a trapping complex, mainly in the product well. Finally, the effect of the zero-point energy constraint on the dynamics properties was analyzed.

  18. Effect of hemodialysis on intraocular lens power calculation.

    PubMed

    Çalışkan, Sinan; Çelikay, Osman; Biçer, Tolga; Aylı, Mehmet Deniz; Gürdal, Canan

    2016-03-01

    Background To evaluate changes in ocular biometric parameters after hemodialysis (HD) in patients with end-stage renal disease (ESRD). Methods Forty eyes of 40 patients undergoing HD were included in this cross-sectional study. Keratometry (K) readings, white-to-white (WTW) distance, central corneal thickness (CCT), anterior chamber depth (ACD), pupil diameter, lens thickness (LT), axial length (AL), and intraocular lens (IOL) power calculation were measured with Lenstar LS 900 (Haag Streit AG, Koeniz, Switzerland) before and after hemodialysis. Intraocular pressure (IOP) was measured with a non-contact tonometer (Tonopachy NT-530P, Nidek Co., LTD, Tokyo, Japan). Main outcomes were changes in biometric parameters after HD. Reliability of the measurements (intraclass correlation coefficients (ICCs)) and the effect size (Cohen's d) were also calculated. Results Mean difference in AL before and after HD was -0.041 ± 0.022 mm with ICCs > 0.90 (p < 0.001 and Cohen's d = 0.06). Pupil diameter was also significantly different before (4.28 ± 0.81 mm) and after (4.44 ± 0.79 mm) HD with ICCs > 0.90 (p = 0.041 and Cohen's d = 0.20). Hemodialysis had no significant effect on K readings, WTW distance, CCT, ACD, LT, or IOP. Conclusion Axial length and pupil diameter increase after HD with small effect size, while HD does not significantly affect IOL power calculations. PMID:26707804

  19. Calculating effective resistances on underlying networks of association schemes

    NASA Astrophysics Data System (ADS)

    Jafarizadeh, M. A.; Sufiani, R.; Jafarizadeh, S.

    2008-07-01

    Recently, in the work of Jafarizadeh et al. [J. Phys, A: Math. Theor. 40, 4949 (2007); e-print arXiv:0705.2480], calculation of effective resistances on distance-regular networks was investigated, where in the first paper, the calculation was based on stratification and Stieltjes functions associated with the network, whereas in the latter one a recursive formula for effective resistances was given based on the Christoffel-Darboux identity. In this paper, evaluation of effective resistances on more general networks that are underlying networks of association schemes is considered, where by using the algebraic combinatoric structures of association schemes such as stratification and Bose-Mesner algebras, an explicit formula for effective resistances on these networks is given in terms of the parameters of the corresponding association schemes. Moreover, we show that for particular underlying networks of association schemes with diameter d such that the adjacency matrix A possesses d +1 distinct eigenvalues, all of the other adjacency matrices Ai, i ≠0, 1 can be written as polynomials of A, i.e., Ai=Pi(A), where Pi is not necessarily of degree i. Then, we use this property for these particular networks and assume that all of the conductances except for one of them, say, c ≡c1=1, are zero to give a procedure for evfor a galuating effective resistances on these networks. The preference of this procedure is that one can evaluate effective resistances by using the structure of their Bose-Mesner algebra without any need to know the spectrum of the adjacency matrices.

  20. SUSY Spectrum Calculator Effects on Dark Matter Observables

    NASA Astrophysics Data System (ADS)

    Bergeron, Paul

    2016-03-01

    Due to questions like the dark matter and hierarchy problem, searches for physics beyond the Standard Model (SM) have not only moved beyond the simplest extensions but the LHC has offered the first tantalizing hints of a new particle at 750 GeV. But as we move into new regimes beyond the minimal extensions to the SM, we need to ask ourselves how accurate are our existing simulations of theory therein? In my talk I will explore the effects of different supersymmetric, spectrum calculators (e.g. ISAJET and SPheno) on predictions for dark matter and Higgs sector observables.

  1. Full potential calculations on the electron bandstructures of Sphalerite, Pyrite and Chalcopyrite

    NASA Astrophysics Data System (ADS)

    Edelbro, R.; Sandström, Å.; Paul, J.

    2003-02-01

    The bulk electronic structures of Sphalerite, Pyrite and Chalcopyrite have been calculated within an ab initio, full potential, density functional approach. The exchange term was approximated with the Dirac exchange functional, the Vosko-Wilk-Nusair parameterization of the Cepler-Alder free electron gas was used for correlation and linear combinations of Gaussian type orbitals were used as basis functions. The Sphalerite (zinc blende) band gap was calculated to be direct with a width of 2.23 eV. The Sphalerite valence band was 5.2 eV wide and composed of a mixture of sulfur and zinc orbitals. The band below the valence band located around -6.2 eV was mainly composed of Zn 3d orbitals. The S 3s orbitals gave rise to a band located around -12.3 eV. Pyrite was calculated to be a semiconductor with an indirect band gap of 0.51 eV, and a direct gap of 0.55 eV. The valence band was 1.25 eV wide and mainly composed of non-bonding Fe 3d orbitals. The band below the valence band was 4.9 eV wide and composed of a mixture of sulfur and iron orbitals. Due to the short inter-atomic distance between the sulfur dumbbells, the S 3s orbitals in Pyrite were split into a bonding and an anti-bonding range. Chalcopyrite was predicted to be a conductor, with no band-crossings at the Fermi level. The bands at -13.2 eV originate from the sulfur 3s orbitals and were quite similar to the sulfur 3s bands in Sphalerite, though somewhat shifted to lower energy. The top of the valence band consisted of a mixture of orbitals from all the atoms. The lower part of the same band showed metal character. Computational modeling as a tool for illuminating the flotation and leaching processes of Pyrite and Chalcopyrite, in connection with surface science experiments, is discussed.

  2. Dechanneling of Positrons in Disordered Lattices Effect of Anharmonic Potential

    NASA Astrophysics Data System (ADS)

    Abu-Assy, M. K.; El-Ashry, M. Y.; Mohamed, A. A.

    2005-01-01

    Dechanneling of positrons due to lattice disorder has been investigated for two stable configurations of the disordered face-centered cubic(fcc) lattices, Dumb-bell configuration (DBC) and Body-centered interstitial (BCI) for channeled positrons with incident energy (10 200) MeV in Cu single crystal in the planar direction (100). The effects of anharmonic terms in the channeling potential have been considered in the calculations. The calculations covered the transition-channeling probability, dechanneling probability, transmission and dechanneling coefficients. It has been found that the transition-channeling probability from the normal into the disordered region occurs only for the transitions n (normal) → n (disordered). Also the dependence of the transmission and dechanneling coefficients on the incident beam position has been studied by using a planar potential function based on shell structure model and compared with the results of a planar potential based on Lindhard's model.

  3. Calculations of the Triton Binding Energy with a Lorentz Boosted Nucleon-Nucleon Potential

    NASA Astrophysics Data System (ADS)

    Kamada, H.; Glöckle, W.; Witała, H.; Golak, J.; Skibiński, R.; Polyzou, W. N.; Elster, Ch.

    2010-04-01

    We study the binding energy of the three-nucleon system in relativistic models that use two different relativistic treatments of the potential that are phase equivalent to realistic NN interactions. One is based on a unitary scale transformation that relates the non-relativistic center-of-mass Hamiltonian to the relativistic mass (rest energy) operator and the other uses a non-linear equation that relates the interaction in the relativistic mass operator to the non-relativistic interaction. In both cases Lorentz-boosted interactions are used in the relativistic Faddeev equation to solve for the three-nucleon binding energy. Using the same realistic NN potentials as input, the solution of the relativistic three-nucleon Faddeev equation for 3H shows slightly less binding energy than the corresponding nonrelativistic result. The effect of the Wigner spin rotation on the binding is very small.

  4. Hybrid density functional calculations of redox potentials and formation energies of transition metal compounds

    NASA Astrophysics Data System (ADS)

    Chevrier, V. L.; Ong, S. P.; Armiento, R.; Chan, M. K. Y.; Ceder, G.

    2010-08-01

    We compare the accuracy of conventional semilocal density functional theory (DFT), the DFT+U method, and the Heyd-Scuseria-Ernzerhof (HSE06) hybrid functional for structural parameters, redox reaction energies, and formation energies of transition metal compounds. Conventional DFT functionals significantly underestimate redox potentials for these compounds. Zhou [Phys. Rev. B 70, 235121 (2004)]10.1103/PhysRevB.70.235121 addressed this issue with DFT+U and a linear-response scheme for calculating U values. We show that the Li intercalation potentials of prominent Li-ion intercalation battery materials, such as the layered LixMO2 ( M=Co and Ni), LixTiS2 ; olivine LixMPO4 ( M=Mn , Fe, Co, and Ni); and spinel-like LixMn2O4 , LixTi2O4 , are also well reproduced by HSE06, due to the self-interaction error correction from the partial inclusion of Hartree-Fock exchange. For formation energies, HSE06 performs well for transition metal compounds, which typically are not well reproduced by conventional DFT functionals but does not significantly improve the results of nontransition metal oxides. Hence, we find that hybrid functionals provide a good alternative to DFT+U for transition metal applications when the large extra computational effort is compensated by the benefits of (i) avoiding species-specific adjustable parameters and (ii) a more universal treatment of the self-interaction error that is not exclusive to specific atomic orbital projections on selected ions.

  5. Calculation of effective doses for broad parallel photon beams.

    PubMed

    Kim, C H; Reece, W D; Poston, J W

    1999-02-01

    Values of effective dose (E) were calculated for the entire range of incident directions of broad parallel photon beams for selected photon energies using the Monte Carlo N-Particle (MCNP) transport code with a hermaphroditic phantom. The calculated results are presented in terms of conversion coefficients transforming air kerma to effective dose. This study also compared the numerical values of E and H(E) over the entire range of incident beam directions. E was always less than H(E) considering all beam directions and photon energies, but the differences were not significant except when a photon beam approaches some specific directions (overhead and underfoot). This result suggests that the current H(E) values can be directly interpreted as E or, at least, as a conservative value of E without knowing the details of irradiation geometries. Finally, based on the distributions of H(E) and E over the beam directions, this study proposes ideal angular response factors for personal dosimeters that can be used to improve the angular response properties of personal dosimeters for off-normal incident photons. PMID:9929126

  6. Molecular electrostatic potential at the atomic sites in the effective core potential approximation.

    PubMed

    Lesiuk, Michał; Zachara, Janusz

    2013-02-21

    Considering calculations of the molecular electrostatic potential at the atomic sites (MEP@AS) in the presence of effective core potentials (ECP), we found that the consequent use of the definition of MEP@AS based on the energy derivative with respect to nuclear charge leads to a formula that differs by one term from the result of simple application of Coulomb's law. We have developed a general method to analytically treat derivatives of ECP with respect to nuclear charge. Benchmarking calculations performed on a set of simple molecules show that our formula leads to a systematic decrease in the error connected with the introduction of ECP when compared to all-electron results. Because of a straightforward implementation and relatively low costs of the developed procedure we suggest to use it by default. PMID:23444997

  7. Molecular electrostatic potential at the atomic sites in the effective core potential approximation

    NASA Astrophysics Data System (ADS)

    Lesiuk, Michał; Zachara, Janusz

    2013-02-01

    Considering calculations of the molecular electrostatic potential at the atomic sites (MEP@AS) in the presence of effective core potentials (ECP), we found that the consequent use of the definition of MEP@AS based on the energy derivative with respect to nuclear charge leads to a formula that differs by one term from the result of simple application of Coulomb's law. We have developed a general method to analytically treat derivatives of ECP with respect to nuclear charge. Benchmarking calculations performed on a set of simple molecules show that our formula leads to a systematic decrease in the error connected with the introduction of ECP when compared to all-electron results. Because of a straightforward implementation and relatively low costs of the developed procedure we suggest to use it by default.

  8. Analysis of water flow paths: methodology and example calculations for a potential geological repository in Sweden.

    PubMed

    Werner, Kent; Bosson, Emma; Berglund, Sten

    2006-12-01

    Safety assessment related to the siting of a geological repository for spent nuclear fuel deep in the bedrock requires identification of potential flow paths and the associated travel times for radionuclides originating at repository depth. Using the Laxemar candidate site in Sweden as a case study, this paper describes modeling methodology, data integration, and the resulting water flow models, focusing on the Quaternary deposits and the upper 150 m of the bedrock. Example simulations identify flow paths to groundwater discharge areas and flow paths in the surface system. The majority of the simulated groundwater flow paths end up in the main surface waters and along the coastline, even though the particles used to trace the flow paths are introduced with a uniform spatial distribution at a relatively shallow depth. The calculated groundwater travel time, determining the time available for decay and retention of radionuclides, is on average longer to the coastal bays than to other biosphere objects at the site. Further, it is demonstrated how GIS-based modeling can be used to limit the number of surface flow paths that need to be characterized for safety assessment. Based on the results, the paper discusses an approach for coupling the present models to a model for groundwater flow in the deep bedrock. PMID:17334050

  9. Numerical calculation of steady inviscid full potential compressible flow about wind turbine blades

    NASA Technical Reports Server (NTRS)

    Dulikravich, D. S.

    1980-01-01

    The air flow through a propeller-type wind turbine rotor is characterized by three-dimensional rotating cascade effects about the inner portions of the rotor blades and compressibility effects about the tip regions of the blades. In the case of large rotor diameter and/or increased rotor angular speed, the existence of small supersonic zones terminated by weak shocks is possible. An exact nonlinear mathematical model (called a steady Full Potential Equation - FPE) that accounts for the above phenomena has been rederived. An artificially time dependent version of FPE was iteratively solved by a finite volume technique involving an artificial viscosity and a three-level consecutive mesh refinement. The exact boundary conditions were applied by generating a boundary conforming periodic computation mesh.

  10. VORSTAB: A computer program for calculating lateral-directional stability derivatives with vortex flow effect

    NASA Technical Reports Server (NTRS)

    Lan, C. Edward

    1985-01-01

    A computer program based on the Quasi-Vortex-Lattice Method of Lan is presented for calculating longitudinal and lateral-directional aerodynamic characteristics of nonplanar wing-body combination. The method is based on the assumption of inviscid subsonic flow. Both attached and vortex-separated flows are treated. For the vortex-separated flow, the calculation is based on the method of suction analogy. The effect of vortex breakdown is accounted for by an empirical method. A summary of the theoretical method, program capabilities, input format, output variables and program job control set-up are described. Three test cases are presented as guides for potential users of the code.

  11. Effects of potential functions on stochastic resonance

    NASA Astrophysics Data System (ADS)

    Li, Jian-Long; Zeng, Ling-Zao

    2011-01-01

    In this paper, the effects of a bistable potential function U(x) = —ax2/2 + b|x|2γ/(2γ) on stochastic resonance (SR) is discussed. We investigate the effects of index γ on the performance of the SR system with fixed parameters a and b, and with fixed potential barriers, respectively. To measure the performance of the SR system in the presence of an aperiodic input, the bit error rate is employed, as is commonly used in binary communications. The numerical simulations strongly support the theoretical results. The goal of this investigation is to explore the effects of the shape of potential functions on SR and give a guidance of nonlinear systems in the application of information processing.

  12. A new hydrocarbon empirical potential in angle bending calculation for the molecular dynamics simulation

    SciTech Connect

    Ping, Tan Ai; Hoe, Yeak Su

    2014-07-10

    Typically, short range potential only depends on neighbouring atoms and its parameters function can be categorized into bond stretching, angle bending and bond rotation potential. In this paper, we present our work called Angle Bending (AB) potential, whereas AB potential is the extension of our previous work namely Bond Stretching (BS) potential. Basically, potential will tend to zero after truncated region, potential in specific region can be represented by different piecewise polynomial. We proposed the AB piecewise potential which is possible to solve a system involving three atoms. AB potential able to handle the potential of covalent bonds for three atoms as well as two atoms cases due to its degeneracy properties. Continuity for the piecewise polynomial has been enforced by coupling with penalty methods. There are still plenty of improvement spaces for this AB potential. The improvement for three atoms AB potential will be studied and further modified into torsional potential which are the ongoing current research.

  13. Environmental Perchlorate Exposure: Potential Adverse Thyroid Effects

    PubMed Central

    Leung, Angela M.; Pearce, Elizabeth N.; Braverman, Lewis E.

    2014-01-01

    Purpose of review This review will present a general overview of the sources, human studies, and proposed regulatory action regarding environmental perchlorate exposure. Recent findings Some recent studies have reported significant associations between urinary perchlorate concentrations, thyroid dysfunction, and decreased infant IQ in groups who would be particularly susceptible to perchlorate effects. An update regarding the recent proposed regulatory actions and potential costs surrounding amelioration of perchlorate contamination is provided. Summary The potential adverse thyroidal effects of environmental perchlorate exposure remain controversial, and further research is needed to further define its relationship to human health among pregnant and lactating women and their infants. PMID:25106002

  14. Electron structure of molecules with very heavy atoms using effective core potentials

    SciTech Connect

    Pitzer, K.S.

    1982-01-01

    Topics covered include effective potential, Hamiltonian for valence-electron motion, molecular calculations, spin-spin coupling, L-S coupling, numerical results of molecular calculations, and results of configuration-interaction Slater-orbital calculations in L-S coupling. (GHT)

  15. The Effects of Relativity on First-Principles Calculations

    NASA Astrophysics Data System (ADS)

    Semi, T.; Mattsson, A. E.; Wills, J. M.

    2011-06-01

    The construction of the equation of state for a given material is of central importance to its characterization. Hugoniots can be calculated using Density Functional Theory (DFT), and DFT points compared to available experimental results. By evaluating the accuracy in a relevant phase space, confidence is gained in the DFT method. This bolsters the dependability of DFT data in phase spaces in which experiment may be difficult or impossible to perform, and verifies its usefulness. The equation of state is comprised of the cold curve and thermal electronic and ionic terms. We discuss differences in the cold curve of Ce produced by first principles calculations using the Scalar Dirac equation with variational spin-orbit coupling treatment and that generated by the full relativistic formulation, both with the same DFT functionals. The relativistic effects of f-electrons in systems like Ce are expected to be of a magnitude substantial enough to be consequential to the description of these structures. ``Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.''

  16. Dynamical properties of non-ideal plasma on the basis of effective potentials

    SciTech Connect

    Ramazanov, T. S.; Kodanova, S. K.; Moldabekov, Zh. A.; Issanova, M. K.

    2013-11-15

    In this work, stopping power has been calculated on the basis of the Coulomb logarithm using the effective potentials. Calculations of the Coulomb logarithm and stopping power for different interaction potentials and degrees of ionization are compared. The comparison with the data of other theoretical and experimental works was carried out.

  17. Calculation of the piezoelectric and flexoelectric effects in nanowires using a decoupled finite element analysis method

    NASA Astrophysics Data System (ADS)

    Zhang, Zhiqiang; Geng, Dalong; Wang, Xudong

    2016-04-01

    A simple and effective decoupled finite element analysis method was developed for simulating both the piezoelectric and flexoelectric effects of zinc oxide (ZnO) and barium titanate (BTO) nanowires (NWs). The piezoelectric potential distribution on a ZnO NW was calculated under three deformation conditions (cantilever, three-point, and four-point bending) and compared to the conventional fully coupled method. The discrepancies of the electric potential maximums from these two methods were found very small, validating the accuracy and effectiveness of the decoupled method. Both ZnO and BTO NWs yielded very similar potential distributions. Comparing the potential distributions induced by the piezoelectric and flexoelectric effects, we identified that the middle segment of a four-point bending NW beam is the ideal place for measuring the flexoelectric coefficient, because the uniform parallel plate capacitor-like potential distribution in this region is exclusively induced by the flexoelectric effect. This decoupled method could provide a valuable guideline for experimental measurements of the piezoelectric effects and flexoelectric effects in the nanometer scale.

  18. Full Dimensional Vibrational Calculations for Methane Using AN Accurate New AB Initio Based Potential Energy Surface

    NASA Astrophysics Data System (ADS)

    Majumder, Moumita; Dawes, Richard; Wang, Xiao-Gang; Carrington, Tucker; Li, Jun; Guo, Hua; Manzhos, Sergei

    2014-06-01

    New potential energy surfaces for methane were constructed, represented as analytic fits to about 100,000 individual high-level ab initio data. Explicitly-correlated multireference data (MRCI-F12(AE)/CVQZ-F12) were computed using Molpro [1] and fit using multiple strategies. Fits with small to negligible errors were obtained using adaptations of the permutation-invariant-polynomials (PIP) approach [2,3] based on neural-networks (PIP-NN) [4,5] and the interpolative moving least squares (IMLS) fitting method [6] (PIP-IMLS). The PESs were used in full-dimensional vibrational calculations with an exact kinetic energy operator by representing the Hamiltonian in a basis of products of contracted bend and stretch functions and using a symmetry adapted Lanczos method to obtain eigenvalues and eigenvectors. Very close agreement with experiment was produced from the purely ab initio PESs. References 1- H.-J. Werner, P. J. Knowles, G. Knizia, 2012.1 ed. 2012, MOLPRO, a package of ab initio programs. see http://www.molpro.net. 2- Z. Xie and J. M. Bowman, J. Chem. Theory Comput 6, 26, 2010. 3- B. J. Braams and J. M. Bowman, Int. Rev. Phys. Chem. 28, 577, 2009. 4- J. Li, B. Jiang and Hua Guo, J. Chem. Phys. 139, 204103 (2013). 5- S Manzhos, X Wang, R Dawes and T Carrington, JPC A 110, 5295 (2006). 6- R. Dawes, X-G Wang, A.W. Jasper and T. Carrington Jr., J. Chem. Phys. 133, 134304 (2010).

  19. Effective potential in density matrix functional theory.

    PubMed

    Nagy, A; Amovilli, C

    2004-10-01

    In the previous paper it was shown that in the ground state the diagonal of the spin independent second-order density matrix n can be determined by solving a single auxiliary equation of a two-particle problem. Thus the problem of an arbitrary system with even electrons can be reduced to a two-particle problem. The effective potential of the two-particle equation contains a term v(p) of completely kinetic origin. Virial theorem and hierarchy of equations are derived for v(p) and simple approximations are proposed. A relationship between the effective potential u(p) of the shape function equation and the potential v(p) is established. PMID:15473719

  20. Accelerating quantum instanton calculations of the kinetic isotope effects

    SciTech Connect

    Karandashev, Konstantin; Vaníček, Jiří

    2015-11-21

    Path integral implementation of the quantum instanton approximation currently belongs among the most accurate methods for computing quantum rate constants and kinetic isotope effects, but its use has been limited due to the rather high computational cost. Here, we demonstrate that the efficiency of quantum instanton calculations of the kinetic isotope effects can be increased by orders of magnitude by combining two approaches: The convergence to the quantum limit is accelerated by employing high-order path integral factorizations of the Boltzmann operator, while the statistical convergence is improved by implementing virial estimators for relevant quantities. After deriving several new virial estimators for the high-order factorization and evaluating the resulting increase in efficiency, using ⋅H{sub α} + H{sub β}H{sub γ} → H{sub α}H{sub β} + ⋅ H{sub γ} reaction as an example, we apply the proposed method to obtain several kinetic isotope effects on CH{sub 4} + ⋅ H ⇌ ⋅ CH{sub 3} + H{sub 2} forward and backward reactions.

  1. Calculations of effective recombination coefficients for nebular astrophysics

    NASA Astrophysics Data System (ADS)

    Fang, Xuan; Liu, Xiaowei; Storey, Pete J.

    2015-08-01

    In the seemingly well established field of nebular astrophysics, there has been a long-standing dichotomy in plasma diagnostics and heavy elemental abundance determinations using the traditional method based on collisionally excited lines on the one hand, and optical recombination lines/continua, on the other. Several mechanisms have been proposed to explain this fundamental problem. Deep spectroscopy and recombination line analysis of emission line nebulae (planetary nebulae and H II regions) in the past decade have pointed to the existence of another previously unknown component of cold, H-deficient material as the culprit. Better constraints are needed on the physical conditions (electron temperature and density), chemical composition, mass, and spatial distribution of the postulated H-deficient inclusions in order to unravel their astrophysical origins. This requires knowledge of the relevant atomic parameters, most importantly the effective recombination coefficients of abundant heavy element ions, such as C II, O II, N II, and Ne II, appropriate for the physical conditions prevailing in those cold inclusions (e.g., electron temperature Te < 1000 K).In this contribution, I will introduce the creation of new effective recombination coefficients for the heavy element optical recombination lines, and review the recent progress in nebular astrophysics since the availability of new and high-quality atomic data. I will also present our new calculations of the effective recombination coefficients for the Ne II recombination line spectrum.

  2. Coupled generator and combustor performance calculations for potential early commercial MHD power plants

    NASA Technical Reports Server (NTRS)

    Dellinger, T. C.; Hnat, J. G.; Marston, C. H.

    1979-01-01

    A parametric study of the performance of the MHD generator and combustor components of potential early commercial open-cycle MHD/steam power plants is presented. Consideration is given to the effects of air heater system concept, MHD combustor type, coal type, thermal input power, oxygen enrichment of the combustion, subsonic and supersonic generator flow and magnetic field strength on coupled generator and combustor performance. The best performance is found to be attained with a 3000 F, indirectly fired air heater, no oxygen enrichment, Illinois no. 6 coal, a two-stage cyclone combustor with 85% slag rejection, a subsonic generator, and a magnetic field configuration yielding a constant transverse electric field of 4 kV/m. Results indicate that optimum net MHD generator power is generally compressor-power-limited rather than electric-stress-limited, with optimum net power a relatively weak function of operating pressure.

  3. WIND: Computer program for calculation of three dimensional potential compressible flow about wind turbine rotor blades

    NASA Technical Reports Server (NTRS)

    Dulikravich, D. S.

    1980-01-01

    A computer program is presented which numerically solves an exact, full potential equation (FPE) for three dimensional, steady, inviscid flow through an isolated wind turbine rotor. The program automatically generates a three dimensional, boundary conforming grid and iteratively solves the FPE while fully accounting for both the rotating cascade and Coriolis effects. The numerical techniques incorporated involve rotated, type dependent finite differencing, a finite volume method, artificial viscosity in conservative form, and a successive line overrelaxation combined with the sequential grid refinement procedure to accelerate the iterative convergence rate. Consequently, the WIND program is capable of accurately analyzing incompressible and compressible flows, including those that are locally transonic and terminated by weak shocks. The program can also be used to analyze the flow around isolated aircraft propellers and helicopter rotors in hover as long as the total relative Mach number of the oncoming flow is subsonic.

  4. Electrode contamination effects of retarding potential analyzer.

    PubMed

    Fang, H K; Oyama, K-I; Cheng, C Z

    2014-01-01

    The electrode contamination in electrostatic analyzers such as Langmuir probes and retarding potential analyzers (RPA) is a serious problem for space measurements. The contamination layer acts as extra capacitance and resistance and leads to distortion in the measured I-V curve, which leads to erroneous measurement results. There are two main effects of the contamination layer: one is the impedance effect and the other is the charge attachment and accumulation due to the capacitance. The impedance effect can be reduced or eliminated by choosing the proper sweeping frequency. However, for RPA the charge accumulation effect becomes serious because the capacitance of the contamination layer is much larger than that of the Langmuir probe of similar dimension. The charge accumulation on the retarding potential grid causes the effective potential, that ions experience, to be changed from the applied voltage. Then, the number of ions that can pass through the retarding potential grid to reach the collector and, thus, the measured ion current are changed. This effect causes the measured ion drift velocity and ion temperature to be changed from the actual values. The error caused by the RPA electrode contamination is expected to be significant for sounding rocket measurements with low rocket velocity (1-2 km/s) and low ion temperature of 200-300 K in the height range of 100-300 km. In this paper we discuss the effects associated with the RPA contaminated electrodes based on theoretical analysis and experiments performed in a space plasma operation chamber. Finally, the development of a contamination-free RPA for sounding rocket missions is presented. PMID:24517809

  5. Symmetry preserving optimised effective potential theory (application to atoms)

    NASA Astrophysics Data System (ADS)

    Theophilou, Andreas K.; Papaconstantinou, Petros G.; Glushkov, Vitaly

    2006-02-01

    In this paper, general symmetry properties of physical systems are used in order to produce a mapping of the external potential of a many electron system to its optimized effective potential (OEP). The so derived effective potential is used to calculate the OEP ground state energies and spin orbitals. Applications are made to atoms and ions and the results are compared to those of the exact Hartree-Fock approximation. The relative deviations from the exact HF theory (triangleE/E) are of the order of 10-4. The same holds for molecules. One of the features of the present theory is that the many electron wave functions derived by the present method, transform according to the irreducible representations of the exact states.

  6. Potential effects on health of global warming

    SciTech Connect

    Haines, A. . Whittington Hospital); Parry, M. . Environmental Change Unit)

    1993-12-01

    Prediction of the impacts of global climate change on health is complicated by a number of factors. These include: the difficulty in predicting regional changes in climate, the capacity for adaptation to climate change, the interactions between the effects of global climate change and a number of other key determinants of health, including population growth and poverty, and the availability of adequate preventive and curative facilities for diseases that may be effected by climate change. Nevertheless, it is of importance to consider the potential health impacts of global climate change for a number of reasons. It is also important to monitor diseases which could be effected by climate change in order to detect changes in incidence as early as possible and study possible interactions with other factors. It seems likely that the possible impacts on health of climate change will be a major determinant of the degree to which policies aimed at reducing global warming are followed, as perceptions of the effect of climate change to human health and well-being are particularly likely to influence public opinion. The potential health impacts of climate change can be divided into direct (primary) and indirect (secondary and tertiary) effects. Primary effects are those related to the effect of temperature on human well-being and disease. Secondary effects include the impacts on health of changes in food production, availability of water and of sea level rise. A tertiary level of impacts can also be hypothesized.

  7. An effective algorithm for calculating the Chandrasekhar function

    NASA Astrophysics Data System (ADS)

    Jablonski, A.

    2012-08-01

    Numerical values of the Chandrasekhar function are needed with high accuracy in evaluations of theoretical models describing electron transport in condensed matter. An algorithm for such calculations should be possibly fast and also accurate, e.g. an accuracy of 10 decimal digits is needed for some applications. Two of the integral representations of the Chandrasekhar function are prospective for constructing such an algorithm, but suitable transformations are needed to obtain a rapidly converging quadrature. A mixed algorithm is proposed in which the Chandrasekhar function is calculated from two algorithms, depending on the value of one of the arguments. Catalogue identifier: AEMC_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEMC_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 567 No. of bytes in distributed program, including test data, etc.: 4444 Distribution format: tar.gz Programming language: Fortran 90 Computer: Any computer with a FORTRAN 90 compiler Operating system: Linux, Windows 7, Windows XP RAM: 0.6 Mb Classification: 2.4, 7.2 Nature of problem: An attempt has been made to develop a subroutine that calculates the Chandrasekhar function with high accuracy, of at least 10 decimal places. Simultaneously, this subroutine should be very fast. Both requirements stem from the theory of electron transport in condensed matter. Solution method: Two algorithms were developed, each based on a different integral representation of the Chandrasekhar function. The final algorithm is edited by mixing these two algorithms and by selecting ranges of the argument ω in which performance is the fastest. Restrictions: Two input parameters for the Chandrasekhar function, x and ω (notation used in the code), are restricted to the range: 0⩽x⩽1 and 0⩽ω⩽1

  8. Effects of Fast Simple Numerical Calculation Training on Neural Systems

    PubMed Central

    Takeuchi, Hikaru; Nagase, Tomomi; Taki, Yasuyuki; Sassa, Yuko; Hashizume, Hiroshi; Nouchi, Rui; Kawashima, Ryuta

    2016-01-01

    Cognitive training, including fast simple numerical calculation (FSNC), has been shown to improve performance on untrained processing speed and executive function tasks in the elderly. However, the effects of FSNC training on cognitive functions in the young and on neural mechanisms remain unknown. We investigated the effects of 1-week intensive FSNC training on cognitive function, regional gray matter volume (rGMV), and regional cerebral blood flow at rest (resting rCBF) in healthy young adults. FSNC training was associated with improvements in performance on simple processing speed, speeded executive functioning, and simple and complex arithmetic tasks. FSNC training was associated with a reduction in rGMV and an increase in resting rCBF in the frontopolar areas and a weak but widespread increase in resting rCBF in an anatomical cluster in the posterior region. These results provide direct evidence that FSNC training alone can improve performance on processing speed and executive function tasks as well as plasticity of brain structures and perfusion. Our results also indicate that changes in neural systems in the frontopolar areas may underlie these cognitive improvements. PMID:26881117

  9. Effect of composition on antiphase boundary energy in Ni3Al based alloys: Ab initio calculations

    NASA Astrophysics Data System (ADS)

    Gorbatov, O. I.; Lomaev, I. L.; Gornostyrev, Yu. N.; Ruban, A. V.; Furrer, D.; Venkatesh, V.; Novikov, D. L.; Burlatsky, S. F.

    2016-06-01

    The effect of composition on the antiphase boundary (APB) energy of Ni-based L 12-ordered alloys is investigated by ab initio calculations employing the coherent potential approximation. The calculated APB energies for the {111} and {001} planes reproduce experimental values of the APB energy. The APB energies for the nonstoichiometric γ' phase increase with Al concentration and are in line with the experiment. The magnitude of the alloying effect on the APB energy correlates with the variation of the ordering energy of the alloy according to the alloying element's position in the 3 d row. The elements from the left side of the 3 d row increase the APB energy of the Ni-based L 12-ordered alloys, while the elements from the right side slightly affect it except Ni. The way to predict the effect of an addition on the {111} APB energy in a multicomponent alloy is discussed.

  10. On Effective Potential in Tortoise Coordinate

    NASA Astrophysics Data System (ADS)

    Ganjali, M. A.

    2012-08-01

    In this paper, we study the field dynamics in Tortoise coordinate where the equation of motion of a scalar can be written as Schrodinger-like form. We obtain a general form for effective potential by finding the Schrodinger equation for scalar and spinor fields and study its global behavior in some black hole backgrounds in three dimension such as BTZ black holes, new type black holes and black holes with no horizon. Especially, we study the asymptotic behavior of potential at infinity, horizons and origin and find that its asymptotic in BTZ and new type solution is completely different from that of vanishing horizon solution. In fact, potential for vanishing horizon goes to a fixed quantity at infinity, while in BTZ and new type black hole we have an infinite barrier.