Calculating activation energies for temperature compensation in circadian rhythms
NASA Astrophysics Data System (ADS)
Bodenstein, C.; Heiland, I.; Schuster, S.
2011-10-01
Many biological species possess a circadian clock, which helps them anticipate daily variations in the environment. In the absence of external stimuli, the rhythm persists autonomously with a period of approximately 24 h. However, single pulses of light, nutrients, chemicals or temperature can shift the clock phase. In the case of light- and temperature-cycles, this allows entrainment of the clock to cycles of exactly 24 h. Circadian clocks have the remarkable property of temperature compensation, that is, the period of the circadian rhythm remains relatively constant within a physiological range of temperatures. For several organisms, temperature-regulated processes within the circadian clock have been identified in recent years. However, how these processes contribute to temperature compensation is not fully understood. Here, we theoretically investigate temperature compensation in general oscillatory systems. It is known that every oscillator can be locally temperature compensated around a reference temperature, if reactions are appropriately balanced. A balancing is always possible if the control coefficient with respect to the oscillation period of at least one reaction in the oscillator network is positive. However, for global temperature compensation, the whole physiological temperature range is relevant. Here, we use an approach which leads to an optimization problem subject to the local balancing principle. We use this approach to analyse different circadian clock models proposed in the literature and calculate activation energies that lead to temperature compensation.
Augustine, R.L.; Lahanas, K.M.; Cole, F.
1992-11-01
An angular overlap calculation has been used to determine the s, p, and d orbital energy levels of the different types of surface sites present on dispersed metal catalysts. These data can permit a Frontier Molecular Orbital treatment of specific site activities as long as the surface orbital availability for overlap with adsorbed substrates is considered along with its energy value and symmetry.
Augustine, R.L.; Lahanas, K.M.; Cole, F.
1992-01-01
An angular overlap calculation has been used to determine the s, p, and d orbital energy levels of the different types of surface sites present on dispersed metal catalysts. These data can permit a Frontier Molecular Orbital treatment of specific site activities as long as the surface orbital availability for overlap with adsorbed substrates is considered along with its energy value and symmetry.
NASA Technical Reports Server (NTRS)
Brown, R. L.; Laufer, A. H.
1981-01-01
Activation energies are calculated by the bond-energy-bond-order (BEBO) and the bond-strength-bond-length (BSBL) methods for the reactions of C2H radicals with H2, CH4, and C2H6 and for the reactions of CN radicals with H2 and CH4. The BSBL technique accurately predicts the activation energies for these reactions while the BEBO method yields energies averaging 9 kcal higher than those observed. A possible reason for the disagreement is considered.
NASA Astrophysics Data System (ADS)
Modine, N. A.; Wright, A. F.; Lee, S. R.
The rate of defect-induced carrier recombination is determined by both defect levels and carrier capture cross-sections. Density functional theory (DFT) has been widely and successfully used to predict defect levels, but only recently has work begun to focus on using DFT to determine carrier capture cross-sections. Lang and Henry developed the theory of carrier-capture by multiphonon emission in the 1970s and showed that carrier-capture cross-sections differ between defects primarily due to differences in their carrier capture activation energies. We present an approach to using DFT to calculate carrier capture activation energies that does not depend on an assumed configuration coordinate and that fully accounts for anharmonic effects, which can substantially modify carrier activation energies. We demonstrate our approach for intrinisic defects in GaAs and GaN and discuss how our results depend on the choice of exchange-correlation functional and the treatment of spin polarization. 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.
NASA Astrophysics Data System (ADS)
Sargsyan, G. N.; Shakhrokh, B.; Harutyunyan, A. B.
2015-02-01
A semi-empirical method is proposed for calculating the activation energy of the unimolecular decomposition of complex compounds using the example of vinyl (ethyl, propyl, and butyl) ethers. The method is based on the concept of the formation of intramolecular hydrogen bonds and the possibility of calculating the energy of deformation of ether molecules upon activation, resulting in the potential surface of the transition state undergoing distortion and the transfer of a hydrogen atom from an alkyl group to a vinyl group. The energy of deformation is calculated using the Mathcad 2001i and MM2 computer programs.
ERIC Educational Resources Information Center
Phillips-Bey, Carol K.
2004-01-01
This article describes TI-73 calculator activities appropriate for middle school students. It was found that the use of the calculator allowed for higher-level thinking and a richer exploration of mathematical ideas by students. [Included with this article are "Dice Roll Worksheet" and "Transforming Tree Worksheet".] (Contains 9 figures.)
Calculating Transition Energy Barriers and Characterizing Activation States for Steps of Fusion.
Ryham, Rolf J; Klotz, Thomas S; Yao, Lihan; Cohen, Fredric S
2016-03-01
We use continuum mechanics to calculate an entire least energy pathway of membrane fusion, from stalk formation, to pore creation, and through fusion pore enlargement. The model assumes that each structure in the pathway is axially symmetric. The static continuum stalk structure agrees quantitatively with experimental stalk architecture. Calculations show that in a stalk, the distal monolayer is stretched and the stored stretching energy is significantly less than the tilt energy of an unstretched distal monolayer. The string method is used to determine the energy of the transition barriers that separate intermediate states and the dynamics of two bilayers as they pass through them. Hemifusion requires a small amount of energy independently of lipid composition, while direct transition from a stalk to a fusion pore without a hemifusion intermediate is highly improbable. Hemifusion diaphragm expansion is spontaneous for distal monolayers containing at least two lipid components, given sufficiently negative diaphragm spontaneous curvature. Conversely, diaphragms formed from single-component distal monolayers do not expand without the continual injection of energy. We identify a diaphragm radius, below which central pore expansion is spontaneous. For larger diaphragms, prior studies have shown that pore expansion is not axisymmetric, and here our calculations supply an upper bound for the energy of the barrier against pore formation. The major energy-requiring deformations in the steps of fusion are: widening of a hydrophobic fissure in bilayers for stalk formation, splay within the expanding hemifusion diaphragm, and fissure widening initiating pore formation in a hemifusion diaphragm. PMID:26958888
Tarr, Lucas; Longcope, Dana; Millhouse, Margaret
2013-06-10
It is well known that photospheric flux emergence is an important process for stressing coronal fields and storing magnetic free energy, which may then be released during a flare. The Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) captured the entire emergence of NOAA AR 11158. This region emerged as two distinct bipoles, possibly connected underneath the photosphere, yet characterized by different photospheric field evolutions and fluxes. The combined active region complex produced 15 GOES C-class, two M-class, and the X2.2 Valentine's Day Flare during the four days after initial emergence on 2011 February 12. The M and X class flares are of particular interest because they are nonhomologous, involving different subregions of the active region. We use a Magnetic Charge Topology together with the Minimum Current Corona model of the coronal field to model field evolution of the complex. Combining this with observations of flare ribbons in the 1600 A channel of the Atmospheric Imaging Assembly on board SDO, we propose a minimization algorithm for estimating the amount of reconnected flux and resulting drop in magnetic free energy during a flare. For the M6.6, M2.2, and X2.2 flares, we find a flux exchange of 4.2 Multiplication-Sign 10{sup 20} Mx, 2.0 Multiplication-Sign 10{sup 20} Mx, and 21.0 Multiplication-Sign 10{sup 20} Mx, respectively, resulting in free energy drops of 3.89 Multiplication-Sign 10{sup 30} erg, 2.62 Multiplication-Sign 10{sup 30} erg, and 1.68 Multiplication-Sign 10{sup 32} erg.
Activation Energy Calculations for the Portevin-Le Chatelier Effect in Nimonic 263 Superalloy
NASA Astrophysics Data System (ADS)
Han, G. M.; Tian, C. G.; Chu, Z. K.; Cui, C. Y.; Hu, Z. Q.; Sun, X. F.
2015-10-01
The Portevin-Le Chatelier (PLC) effect in the Nimonic 263 superalloy was investigated by tensile test in a wide temperature range, from room temperature to 1033 K (760 °C), and at strain rates between 0.1 and 4 × 10-4 s-1. Types A, B, and C serrations were observed depending upon the test temperatures and strain rates. The activation energy ( Q) for serrated flow was determined by employing various methodologies for T < 723 K (450 °C), where a normal PLC effect was observed. The value of Q was found to be independent of the method employed. The average Q value of 70 KJ/mol was found to be in agreement with that for diffusion of substitutional solutes such as Cr and Mo in a Ni matrix by pipe diffusion. At temperatures ranging from 723 K to 923 K (450 °C to 650 °C), type C serrations and an inverse PLC effect were noticed, which may be caused by unlocking process.
Dushanov, E; Kholmurodov, Kh; Yasuoka, K
2013-01-01
Formamide contains the four elements (C, H, O, and N) most required for life and it is attractive as a potential prebiotic starting material for nucleobase synthesis. In the presence of catalysts (for example, TiO2) and with moderate heating, formamide can pass surface energy barriers, yielding a complete set of nucleic bases and acyclonucleosides, and favoring both phosphorylations and transphosphorylations necessary for life. In the reaction mechanism, interaction with water seems to be an essential factor for the formamide molecule to function. In this paper, a formamide–water solution on a TiO$_2$ (anatase) surface is simulated using the molecular dynamics method, and activation energy calculations are performed for the temperature range of T = 250 K to T = 400 K. A correlation is established between the diffusion and density profiles for the formamide and water molecules on an anatase surface. Also, the calculated activation energies of the formamide–water–anatase and formamide–water–platinum systems are compared. A comparative analysis is performed of the behavior of formamide–water and ethanol–water interaction on the same (anatase and platinum) surfaces. PMID:23802018
Biogas - the calculable energy
NASA Astrophysics Data System (ADS)
Kith, Károly; Nagy, Orsolya; Balla, Zoltán; Tamás, András
2015-04-01
EU actions against climate change are rising energy prices, both have emphasized the use of renewable energy,increase investments and energy efficiency. A number of objectives formulated in the EC decree no. 29/2009 by 2020. This document is based on the share of renewable energies in energy consumption should be increased to 20% (EC, 2009). The EU average is 20% but the share of renewables vary from one member state to another. In Hungary in 2020, 14.65% renewable energy share is planned to be achieved. According to the latest Eurostat data, the share of renewable energy in energy consumption of the EU average was 14.1%, while in Hungary, this share was 9.6% in 2012. (EUROSTAT, 2014). The use of renewable energy plant level is influenced by several factors. The most important of these is the cost savings and efficiency gains. Hungarian investments in renewable energy production usually have high associated costs and the payback period is substantially more than five years, depending on the support rate. For example, the payback period is also influenced by the green electricity generated feed prices, which is one of the lowest in Hungary compared the Member States of the European Union. Consequently, it is important to increase the production of green energy. Nowadays, predictable biogas energy is an outstanding type of decentralized energy production. It follows directly that agricultural by-products can be used to produce energy and they also create jobs by the construction of a biogas plant. It is important to dispose of and destroy hazardous and noxious substances in energy production. It follows from this that the construction of biogas plants have a positive impact, in addition to green energy which is prepared to reduce the load on the environment. The production of biogas and green electricity is one of the most environment friendly forms of energy production. Biogas production also has other important ecological effects, such as the substitution of
Activities: More Calculator Capers.
ERIC Educational Resources Information Center
Schmalz, Rosemary
1983-01-01
Provided is an activity designed to give grades 7-12 students opportunities to discover numerical patterns and to derive general conclusions from observing data. The activity focuses attention on patterns in products such as 33x34, 333x334, and 3333x3334. Three worksheets and answers are included. (JN)
Density-functional calculations of prefactors and activation energies for H diffusion in BaZrO3
NASA Astrophysics Data System (ADS)
Sundell, Per G.; Björketun, Mårten E.; Wahnström, Göran
2007-09-01
Density-functional calculations are used to investigate hydrogen diffusion in the solid-state proton conductor BaZrO3 . Activation energies and prefactors for the rate of proton transfer and reorientation are evaluated for a defect-free region of this simple cubic perovskite-structured oxide. Both semiclassical over-barrier jumps and phonon-assisted tunneling transitions between sites are considered. It is found that the classical barriers for the elementary transfer and reorientation steps are both of the order of 0.2eV . The quantum-mechanical zero-point motion effects are found to be sizable, to effectively reduce the barrier heights, and to make the prefactors similar for the transfer and reorientation steps. The Flynn-Stoneham model [Phys. Rev. B 1, 3966 (1970)] of phonon-assisted tunneling yields an activation energy of around 0.2eV and a very small prefactor for proton transfer, whereas the corresponding adiabatic model gives a similar activation energy but a much larger prefactor. It is suggested that the effect of other defects such as dopants has to be included for a proper description of hydrogen diffusion in this material.
Latkowski, J.F.; Sanz, J.; Vujic, J.L.
1996-06-26
Inertial fusion energy (IFE) and magnetic fusion energy (MFE) power plants will probably operate in a pulsed mode. The two different schemes, however, will have quite different time periods. Typical repetition rates for IFE power plants will be 1-5 Hz. MFE power plants will ramp up in current for about 1 hour, shut down for several minutes, and repeat the process. Traditionally, activation calculations for IFE and MFE power plants have assumed continuous operation and used either the ``steady state`` (SS) or ``equivalent steady state`` (ESS) approximations. It has been suggested recently that the SS and ESS methods may not yield accurate results for all radionuclides of interest. The present work expands that of Sisolak, et al. by applying their formulae to conditions which might be experienced in typical IFE and MFE power plants. In addition, complicated, multi-step reaction/decay chains are analyzed using an upgraded version of the ACAB radionuclide generation/depletion code. Our results indicate that the SS method is suitable for application to MFE power plant conditions. We also find that the ESS method generates acceptable results for radionuclides with half-lives more than a factor of three greater than the time between pulses. For components that are subject to 0.05 Hz (or more frequent) irradiation (such as coolant), use of the ESS method is recommended. For components or materials that are subject to less frequent irradiation (such as high-Z target materials), pulsed irradiation calculations should be used.
Maurer, Manuela; de Beer, Stephanie B A; Oostenbrink, Chris
2016-01-01
The periplasmic oligopeptide binding protein A (OppA) represents a well-known example of water-mediated protein-ligand interactions. Here, we perform free-energy calculations for three different ligands binding to OppA, using a thermodynamic integration approach. The tripeptide ligands share a high structural similarity (all have the sequence KXK), but their experimentally-determined binding free energies differ remarkably. Thermodynamic cycles were constructed for the ligands, and simulations conducted in the bound and (freely solvated) unbound states. In the unbound state, it was observed that the difference in conformational freedom between alanine and glycine leads to a surprisingly slow convergence, despite their chemical similarity. This could be overcome by increasing the softness parameter during alchemical transformations. Discrepancies remained in the bound state however, when comparing independent simulations of the three ligands. These difficulties could be traced to a slow relaxation of the water network within the active site. Fluctuations in the number of water molecules residing in the binding cavity occur mostly on a timescale larger than the simulation time along the alchemical path. After extensive simulations, relative binding free energies that were converged to within thermal noise could be obtained, which agree well with available experimental data. PMID:27092480
Neutron and deuteron activation calculations for IFMIF
NASA Astrophysics Data System (ADS)
Forrest, R. A.; Loughlin, M. J.
2007-08-01
The materials for future fusion devices such as DEMO require testing to high neutron fluence. Such testing is planned to be carried out in IFMIF, an accelerator based facility where the neutrons will have maximum energy of about 55 MeV, but with a broad peak near 14 MeV. In order that activation calculations for IFMIF can be carried out, the nuclear data must contain cross sections covering a similar energy range. A description of the EASY-2005 system is given and it is noted that a new library has been added to EASY to cover another significant source of activation from deuteron-induced reactions. Calculations of the neutron activation of materials in many regions of IFMIF have been carried out. These calculations are reported, and the contribution of neutrons above 20 MeV to the activation is discussed. Preliminary calculations using the deuteron library have been made and the activation from deuterons is discussed.
O'Regan, Brian C; Durrant, James R
2006-05-01
Transient photovoltage and photocurrent measurements have been employed to determine the recombination and transport kinetics in operating dye-sensitized photovoltaic cells as a function of potential and temperature. Photocurrent transients have been taken at the open circuit potential, as opposed to the standard measurement at short circuit. Kinetic results have been used to calculate the activation energy as function of the Fermi level position in the TiO(2). In the calculation of activation energies, we have explicitly taken into account the temperature dependence of the offset between the electrolyte redox potential and the conduction band edge. This new method gives activation energies that decrease linearly as the Fermi level position moves toward the conduction band edge, as expected, but not found in previous studies. The results are consistent with the presence of a distribution of traps below the TiO(2) conduction band, the detrapping from which limits both the transport and the recombination of electrons. PMID:16640403
Light Pipe Energy Savings Calculator
NASA Astrophysics Data System (ADS)
Owens, Erin; Behringer, Ernest R.
2009-04-01
Dependence on fossil fuels is unsustainable and therefore a shift to renewable energy sources such as sunlight is required. Light pipes provide a way to utilize sunlight for interior lighting, and can reduce the need for fossil fuel-generated electrical energy. Because consumers considering light pipe installation may be more strongly motivated by cost considerations than by sustainability arguments, an easy means to examine the corresponding costs and benefits is needed to facilitate informed decision-making. The purpose of this American Physical Society Physics and Society Fellowship project is to create a Web-based calculator to allow users to quantify the possible cost savings for their specific light pipe application. Initial calculations show that the illumination provided by light pipes can replace electric light use during the day, and in many cases can supply greater illumination levels than those typically given by electric lighting. While the installation cost of a light pipe is significantly greater than the avoided cost of electricity over the lifetime of the light pipe at current prices, savings may be realized if electricity prices increase.
ACTIVE: a program to calculate and plot reaction rates from ANISN calculated fluxes
Judd, J.L.
1981-12-01
The ACTIVE code calculates spatial heating rates, tritium production rates, neutron reaction rates, and energy spectra from particle fluxes calculated by ANISN. ACTIVE has a variety of input options including the capability to plot all calculated spatial distributions. The code was primarily designed for use with fusion first wall/blanket systems, but could be applied to any one-dimensional problem.
Iwatani, K; Hoshi, M; Shizuma, K; Hiraoka, M; Hayakawa, N; Oka, T; Hasai, H
1994-10-01
A benchmark test of the Monte Carlo neutron and photon transport code system (MCNP) was performed using a bare- and energy-moderated 252Cf fission neutron source which was obtained by transmission through 10-cm-thick iron. An iron plate was used to simulate the effect of the Hiroshima atomic bomb casing. This test includes the activation of indium and nickel for fast neutrons and gold, europium, and cobalt for thermal and epithermal neutrons, which were inserted in the moderators. The latter two activations are also to validate 152Eu and 60Co activity data obtained from the atomic bomb-exposed specimens collected at Hiroshima and Nagasaki, Japan. The neutron moderators used were Lucite and Nylon 6 and the total thickness of each moderator was 60 cm or 65 cm. Measured activity data (reaction yield) of the neutron-irradiated detectors in these moderators decreased to about 1/1,000th or 1/10,000th, which corresponds to about 1,500 m ground distance from the hypocenter in Hiroshima. For all of the indium, nickel, and gold activity data, the measured and calculated values agreed within 25%, and the corresponding values for europium and cobalt were within 40%. From this study, the MCNP code was found to be accurate enough for the bare- and energy-moderated 252Cf neutron activation calculations of these elements using moderators containing hydrogen, carbon, nitrogen, and oxygen. PMID:8083048
Iwatani, Kazuo; Shizuma, Kiyoshi; Hasai, Hiromi; Hoshi, Masaharu; Hiraoka, Masayuki; Hayakawa, Norihiko; Oka, Takamitsu
1994-10-01
A benchmark test of the Monte Carlo neutron and photon transport code system (MCNP) was performed using a bare- and energy-moderated {sup 252}Cf fission neutron source which was obtained by transmission through 10-cm-thick iron. An iron plate was used to simulate the effect of the Hiroshima atomic bomb casing. This test includes the activation of indium and nickel for fast neutrons and gold, europium, and cobalt for thermal and epithermal neutrons, which were inserted in the moderators. The latter two activations are also to validate {sup 152}Eu and {sup 60}Co activity data obtained from the atomic bomb-exposed specimens collected at Hiroshima and Nagasaki, Japan. The neutron moderators used were Lucite and Nylon 6 and the total thickness of each moderator was 60 cm or 65 cm. Measured activity data (reaction yield) of the neutron-irradiated detectors in these moderators decreased to about 1/1,000th or 1/10,000th, which corresponds to about 1,500 m ground distance from the hypocenter in Hiroshima. For all of the indium, nickel, and gold activity data, the measured and calculated values agreed within 25%, and the corresponding values for europium and cobalt were within 40%. From this study, the MCNP code was found to be accurate enough for the bare- and energy-moderated {sup 252}Cf neutron activation calculations of these elements using moderators containing hydrogen, carbon, nitrogen, and oxygen. 18 refs., 10 figs., 4 tabs.
Calculation of water activation for the LHC
NASA Astrophysics Data System (ADS)
Vollaire, Joachim; Brugger, Markus; Forkel-Wirth, Doris; Roesler, Stefan; Vojtyla, Pavol
2006-06-01
The management of activated water in the Large Hadron Collider (LHC) at CERN is a key concern for radiation protection. For this reason, the induced radioactivity of the different water circuits is calculated using the Monte-Carlo (MC) code FLUKA. The results lead to the definition of procedures to be taken into account during the repair and maintenance of the machine, as well as to measures being necessary for a release of water into the environment. In order to assess the validity of the applied methods, a benchmark experiment was carried out at the CERN-EU High Energy Reference Field (CERF) facility, where a hadron beam (120 GeV) is impinging on a copper target. Four samples of water, as used at the LHC, and different in their chemical compositions, were irradiated near the copper target. In addition to the tritium activity measured with a liquid scintillation counter, the samples were also analyzed using gamma spectroscopy in order to determine the activity of the gamma emitting isotopes such as Be7 and Na24. While for the latter an excellent agreement between simulation and measurement was found, for the calculation of tritium a correction factor is derived to be applied for future LHC calculations in which the activity is calculated by direct scoring of produced nuclei. A simplified geometry representing the LHC Arc sections is then used to evaluate the different calculation methods with FLUKA. By comparing these methods and by taking into account the benchmark results, a strategy for the environmental calculations can be defined.
Yao, Yuan; Liu, Junjun; Zhan, Chang-Guo
2012-11-01
Human butyrylcholinesterase (BChE) is recognized as the most promising bioscavenger for organophosphorus (OP) warfare nerve agents. The G117H mutant of human BChE has been identified as a potential catalytic bioscavenger with a remarkably improved activity against OP nerve agents such as sarin, but it still does not satisfy the clinical use. For further design of the higher-activity mutants against OP nerve agents, it is essential to understand how the G117H mutation improves the activity. The reaction mechanisms and the free energy profiles for spontaneous reactivation of wild-type BChE and its G117H mutant phosphorylated by sarin have been explored, in this study, by performing first-principles quantum mechanical/molecular mechanical free energy calculations, and the remarkable role of the G117H mutation on the activity has been elucidated. For both the wild-type and G117H mutant enzymes, H438 acts as a general base to initiate the spontaneous reactivation that consists of two reaction steps: the nucleophilic attack at the phosphorus by a water molecule and decomposition of the pentacoordinated phosphorus intermediate. The calculated overall free energy barriers, i.e., 30.2 and 23.9 kcal/mol for the wild type and G117H mutant, respectively, are in good agreement with available experimental kinetic data. On the basis of the calculated results, the mutated residue (H117 in the G117H mutant) cannot initiate the spontaneous reactivation as a general base. Instead, it skews the oxyanion hole and makes the phosphorus more open to the nucleophilic water molecule, resulting in a remarkable change in the rate-determining step and significantly improved catalytic activity of human BChE. PMID:23092211
Mathematical Creative Activity and the Graphic Calculator
ERIC Educational Resources Information Center
Duda, Janina
2011-01-01
Teaching mathematics using graphic calculators has been an issue of didactic discussions for years. Finding ways in which graphic calculators can enrich the development process of creative activity in mathematically gifted students between the ages of 16-17 is the focus of this article. Research was conducted using graphic calculators with…
Niu, Yuzhen; Pan, Dabo; Shi, Danfeng; Bai, Qifeng; Liu, Huanxiang; Yao, Xiaojun
2015-01-01
As a promising target for the treatment of lung cancer, the MutT Homolog 1 (MTH1) protein can be inhibited by crizotinib. A recent work shows that the inhibitory potency of (S)-crizotinib against MTH1 is about 20 times over that of (R)-crizotinib. But the detailed molecular mechanism remains unclear. In this study, molecular dynamics (MD) simulations and free energy calculations were used to elucidate the mechanism about the effect of chirality of crizotinib on the inhibitory activity against MTH1. The binding free energy of (S)-crizotinib predicted by the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and Adaptive biasing force (ABF) methodologies is much lower than that of (R)-crizotinib, which is consistent with the experimental data. The analysis of the individual energy terms suggests that the van der Waals interactions are important for distinguishing the binding of (S)-crizotinib and (R)-crizotinib. The binding free energy decomposition analysis illustrated that residues Tyr7, Phe27, Phe72 and Trp117 were important for the selective binding of (S)-crizotinib to MTH1. The adaptive biasing force (ABF) method was further employed to elucidate the unbinding process of (S)-crizotinib and (R)-crizotinib from the binding pocket of MTH1. ABF simulation results suggest that the reaction coordinates of the (S)-crizotinib from the binding pocket is different from (R)-crizotinib. The results from our study can reveal the details about the effect of chirality on the inhibition activity of crizotinib to MTH1 and provide valuable information for the design of more potent inhibitors. PMID:26677850
NASA Astrophysics Data System (ADS)
Vinograd, Victor L.; Burton, Benjamin P.; Gale, Julian D.; Allan, Neil L.; Winkler, Björn
2007-02-01
Thermodynamic mixing properties and subsolidus phase relations of the rhombohedral carbonate system, (1 - x) · CaCO 3 - x · MgCO 3, were modelled in the temperature range of 623-2023 K with static structure energy calculations based on well-parameterised empirical interatomic potentials. Relaxed static structure energies of a large set of randomly varied structures in a 4 × 4 × 1 supercell of R3¯c calcite ( a = 19.952 Å, c = 17.061 Å) were calculated with the General Utility Lattice Program (GULP). These energies were cluster expanded in a basis set of 12 pair-wise effective interactions. Temperature-dependent enthalpies of mixing were calculated by the Monte Carlo method. Free energies of mixing were obtained by thermodynamic integration of the Monte Carlo results. The calculated phase diagram is in good agreement with experimental phase boundaries.
Hladilkova, Jana; Prokop, Zbynek; Chaloupkova, Radka; Damborsky, Jiri; Jungwirth, Pavel
2013-11-21
Release of halide ions is an essential step of the catalytic cycle of haloalkane dehalogenases. Here we describe experimentally and computationally the process of release of a halide anion from the buried active site of the haloalkane dehalogenase LinB. Using stopped-flow fluorescence analysis and umbrella sampling free energy calculations, we show that the anion binding is ion-specific and follows the ordering I(-) > Br(-) > Cl(-). We also address the issue of the protonation state of the catalytic His272 residue and its effect on the process of halide release. While deprotonation of His272 increases binding of anions in the access tunnel, we show that the anionic ordering does not change with the switch of the protonation state. We also demonstrate that a sodium cation could relatively easily enter the active site, provided the His272 residue is singly protonated, and replace thus the missing proton. In contrast, Na(+) is strongly repelled from the active site containing the doubly protonated His272 residue. Our study contributes toward understanding of the reaction mechanism of haloalkane dehalogenase enzyme family. Determination of the protonation state of the catalytic histidine throughout the catalytic cycle remains a challenge for future studies. PMID:24151979
Good Practices in Free-energy Calculations
NASA Technical Reports Server (NTRS)
Pohorille, Andrew; Jarzynski, Christopher; Chipot, Christopher
2013-01-01
As access to computational resources continues to increase, free-energy calculations have emerged as a powerful tool that can play a predictive role in drug design. Yet, in a number of instances, the reliability of these calculations can be improved significantly if a number of precepts, or good practices are followed. For the most part, the theory upon which these good practices rely has been known for many years, but often overlooked, or simply ignored. In other cases, the theoretical developments are too recent for their potential to be fully grasped and merged into popular platforms for the computation of free-energy differences. The current best practices for carrying out free-energy calculations will be reviewed demonstrating that, at little to no additional cost, free-energy estimates could be markedly improved and bounded by meaningful error estimates. In energy perturbation and nonequilibrium work methods, monitoring the probability distributions that underlie the transformation between the states of interest, performing the calculation bidirectionally, stratifying the reaction pathway and choosing the most appropriate paradigms and algorithms for transforming between states offer significant gains in both accuracy and precision. In thermodynamic integration and probability distribution (histogramming) methods, properly designed adaptive techniques yield nearly uniform sampling of the relevant degrees of freedom and, by doing so, could markedly improve efficiency and accuracy of free energy calculations without incurring any additional computational expense.
Calculating Free Energies Using Average Force
NASA Technical Reports Server (NTRS)
Darve, Eric; Pohorille, Andrew; DeVincenzi, Donald L. (Technical Monitor)
2001-01-01
A new, general formula that connects the derivatives of the free energy along the selected, generalized coordinates of the system with the instantaneous force acting on these coordinates is derived. The instantaneous force is defined as the force acting on the coordinate of interest so that when it is subtracted from the equations of motion the acceleration along this coordinate is zero. The formula applies to simulations in which the selected coordinates are either unconstrained or constrained to fixed values. It is shown that in the latter case the formula reduces to the expression previously derived by den Otter and Briels. If simulations are carried out without constraining the coordinates of interest, the formula leads to a new method for calculating the free energy changes along these coordinates. This method is tested in two examples - rotation around the C-C bond of 1,2-dichloroethane immersed in water and transfer of fluoromethane across the water-hexane interface. The calculated free energies are compared with those obtained by two commonly used methods. One of them relies on determining the probability density function of finding the system at different values of the selected coordinate and the other requires calculating the average force at discrete locations along this coordinate in a series of constrained simulations. The free energies calculated by these three methods are in excellent agreement. The relative advantages of each method are discussed.
NASA Technical Reports Server (NTRS)
Gadeken, Owen
2002-01-01
Teaming is so common in today's project management environment that most of us assume it comes naturally. We further assume that when presented with meaningful and challenging work, project teams will naturally engage in productive activity to complete their tasks. This assumption is expressed in the simple (but false) equation: Team + Work = Teamwork. Although this equation appears simple and straightforward, it is far from true for most project organizations whose reality is a complex web of institutional norms based on individual achievement and rewards. This is illustrated by the very first successful team experience from my early Air Force career. As a young lieutenant, I was sent to Squadron Officer School, which was the first in the series of Air Force professional military education courses I was required to complete during my career. We were immediately formed into teams of twelve officers. Much of the course featured competition between these teams. As the most junior member of my team, I quickly observed the tremendous pressure to show individual leadership capability. At one point early in the course, almost everyone in our group was vying to become the team leader. This conflict was so intense that it caused us to fail miserably in our first outdoor team building exercise. We spent so much time fighting over leadership that we were unable to complete any of the events on the outdoor obstacle course. This complete lack of success was so disheartening to me that I gave our team little hope for future success. What followed was a very intense period of bickering, conflict, and even shouting matches as our dysfunctional team tried to cope with our early failures and find some way to succeed. British physician and researcher Wilfred Bion (Experiences in Groups, 1961) discovered that there are powerful psychological forces inherent in all groups that divert from accomplishing their primary tasks. To overcome these restraining forces and use the potential
Predicting proteinase specificities from free energy calculations.
Mekonnen, Seble Merid; Olufsen, Magne; Smalås, Arne O; Brandsdal, Bjørn O
2006-10-01
The role of the primary binding residue (P1) in complexes between three different subtilases (subtilisin Carlsberg, thermitase and proteinase K) and their canonical protein inhibitor eglin c have been studied by free energy calculations. Based on the crystal structures of eglin c in complex with subtilisin Carlsberg and thermitase, and a homology model of the eglin c-proteinase K complex, a total of 57 mutants have been constructed and docked into their host proteins. The binding free energy was then calculated using molecular dynamics (MD) simulations combined with the linear interaction energy (LIE) method for all complexes differing only in the nature of the amino acid at the P1 position. LIE calculations for 19 different complexes for each subtilase were thus carried out excluding proline. The effects of substitutions at the P1 position on the binding free energies are found to be very large, and positively charged residues (Arg, Lys and His) are particularly deleterious for all three enzymes. The charged variants of the acidic side chains are found to bind more favorably as compared to their protonated states in all three subtilases. Furthermore, hydrophobic amino acids are accommodated most favorably at the S1-site in all three enzymes. Comparison of the three series of binding free energies shows only minor differences in the 19 computed relative binding free energies among these subtilases. This is further reflected in the correlation coefficient between the 23 relative binding free energies obtained, including the possible protonation states of ionizable side chains, but excluding the P1 Pro, for subtilisin Carlsberg versus thermitase (0.95), subtilisin versus proteinase K (0.94) and thermitase versus proteinase K (0.96). PMID:16386933
Hermes III endpoint energy calculation from photonuclear activation of ^{197}Au and ^{58}Ni foils
Parzyck, Christopher Thomas
2014-09-01
A new process has been developed to characterize the endpoint energy of HERMES III on a shot-to-shot basis using standard dosimetry tools from the Sandia Radiation Measurements Laboratory. Photonuclear activation readings from nickel and gold foils are used in conjunction with calcium fluoride thermoluminescent dosimeters to derive estimated electron endpoint energies for a series of HERMES shots. The results are reasonably consistent with the expected endpoint voltages on those shots.
Bond-Energy and Surface-Energy Calculations in Metals
ERIC Educational Resources Information Center
Eberhart, James G.; Horner, Steve
2010-01-01
A simple technique appropriate for introductory materials science courses is outlined for the calculation of bond energies in metals from lattice energies. The approach is applied to body-centered cubic (bcc), face-centered cubic (fcc), and hexagonal-closest-packed (hcp) metals. The strength of these bonds is tabulated for a variety metals and is…
Zero energy scattering calculation in Euclidean space
NASA Astrophysics Data System (ADS)
Carbonell, J.; Karmanov, V. A.
2016-03-01
We show that the Bethe-Salpeter equation for the scattering amplitude in the limit of zero incident energy can be transformed into a purely Euclidean form, as it is the case for the bound states. The decoupling between Euclidean and Minkowski amplitudes is only possible for zero energy scattering observables and allows determining the scattering length from the Euclidean Bethe-Salpeter amplitude. Such a possibility strongly simplifies the numerical solution of the Bethe-Salpeter equation and suggests an alternative way to compute the scattering length in Lattice Euclidean calculations without using the Luscher formalism. The derivations contained in this work were performed for scalar particles and one-boson exchange kernel. They can be generalized to the fermion case and more involved interactions.
Free-Energy Calculations. A Mathematical Perspective
NASA Technical Reports Server (NTRS)
Pohorille, Andrzej
2015-01-01
conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.
Ultrasonic energy in liposome production: process modelling and size calculation.
Barba, A A; Bochicchio, S; Lamberti, G; Dalmoro, A
2014-04-21
The use of liposomes in several fields of biotechnology, as well as in pharmaceutical and food sciences is continuously increasing. Liposomes can be used as carriers for drugs and other active molecules. Among other characteristics, one of the main features relevant to their target applications is the liposome size. The size of liposomes, which is determined during the production process, decreases due to the addition of energy. The energy is used to break the lipid bilayer into smaller pieces, then these pieces close themselves in spherical structures. In this work, the mechanisms of rupture of the lipid bilayer and the formation of spheres were modelled, accounting for how the energy, supplied by ultrasonic radiation, is stored within the layers, as the elastic energy due to the curvature and as the tension energy due to the edge, and to account for the kinetics of the bending phenomenon. An algorithm to solve the model equations was designed and the relative calculation code was written. A dedicated preparation protocol, which involves active periods during which the energy is supplied and passive periods during which the energy supply is set to zero, was defined and applied. The model predictions compare well with the experimental results, by using the energy supply rate and the time constant as fitting parameters. Working with liposomes of different sizes as the starting point of the experiments, the key parameter is the ratio between the energy supply rate and the initial surface area. PMID:24647821
18 CFR 11.13 - Energy gains calculations.
Code of Federal Regulations, 2010 CFR
2010-04-01
... calculations. 11.13 Section 11.13 Conservation of Power and Water Resources FEDERAL ENERGY REGULATORY... FEDERAL POWER ACT Charges for Headwater Benefits § 11.13 Energy gains calculations. (a) Energy gains at a... not complex or in which headwater benefits are expected to be small, calculations will be made...
A Programmable Calculator Activity, x = 1/x + 1.
ERIC Educational Resources Information Center
Snover, Stephen L.; Spikell, Mark A.
An activity for secondary schools is presented and discussed which may be explored with a programmable calculator. The activity is non-standard and could not be easily explored without the use of a programmable calculator. Related activities are also discussed. Flow charts and programs for different programmable calculators are presented. (MP)
Energy Band Calculations for Maximally Even Superlattices
NASA Astrophysics Data System (ADS)
Krantz, Richard; Byrd, Jason
2007-03-01
Superlattices are multiple-well, semiconductor heterostructures that can be described by one-dimensional potential wells separated by potential barriers. We refer to a distribution of wells and barriers based on the theory of maximally even sets as a maximally even superlattice. The prototypical example of a maximally even set is the distribution of white and black keys on a piano keyboard. Black keys may represent wells and the white keys represent barriers. As the number of wells and barriers increase, efficient and stable methods of calculation are necessary to study these structures. We have implemented a finite-element method using the discrete variable representation (FE-DVR) to calculate E versus k for these superlattices. Use of the FE-DVR method greatly reduces the amount of calculation necessary for the eigenvalue problem.
Groundwater activation calculations for E872
Freeman, W.S.; E872 Collaboration
1995-08-01
The E872 beam dump geometry has been modeled in CASIM and calculations have been done to determine the annual limits for protons n target. Results are presented using both the single resident well model (SRWM) and the newly-approved concentration model (CM). The conclusion is that the target/dump design is adequate for the maximum number of protons on target requested by the experiment, which is >1 {times} 10{sup 18} protons per year at 800 GeV.
Low-energy calculations for nuclear photodisintegration
NASA Astrophysics Data System (ADS)
Deflorian, S.; Efros, V. D.; Leidemann, W.
2016-03-01
In the Standard Solar Model a central role in the nucleosynthesis is played by reactions of the kind {}{Z_1}{A_1}{X_1} + {}{Z_2}{A_2}{X_2} to {}{Z_1 + {Z_2}}{A_1 + {A_2}}Y + γ , which enter the proton-proton chains. These reactions can also be studied through the inverse photodisintegration reaction. One option is to use the Lorentz Integral Transform approach, which transforms the continuum problem into a bound state-like one. A way to check the reliability of such methods is a direct calculation, for example using the Kohn Variational Principle to obtain the scattering wave function and then directly calculate the response function of the reaction.
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.
Calculation of energy deposition distributions for simple geometries
NASA Technical Reports Server (NTRS)
Watts, J. W., Jr.
1973-01-01
When high-energy charged particles pass through a thin detector, the ionization energy loss in that detector is subject to fluctuations or straggling which must be considered in interpreting the data. Under many conditions, which depend upon the charge and energy of the incident particle and the detector geometry, the ionization energy lost by the particle is significantly different from the energy deposited in the detector. This problem divides naturally into a calculation of the energy loss that results in excitation and low-energy secondary electrons which do not travel far from their production points, and a calculation of energy loss that results in high-energy secondary electrons which can escape from the detector. The first calculation is performed using a modification of the Vavilov energy loss distribution. A cutoff energy is introduced above which all electrons are ignored and energy transferred to low energy particles is assumed to be equivalent to the energy deposited by them. For the second calculation, the trajectory of the primary particle is considered as a source of secondary high-energy electrons. The electrons from this source are transported using Monte Carlo techniques and multiple scattering theory, and the energy deposited by them in the detector is calculated. The results of the two calculations are then combined to predict the energy deposition distribution. The results of these calculations are used to predict the charge resolution of parallel-plate pulse ionization chambers that are being designed to measure the charge spectrum of heavy nuclei in the galactic cosmic-ray flux.
Energy assessment: physical activity
Technology Transfer Automated Retrieval System (TEKTRAN)
Physical activity is an important component of total energy expenditure, contributing to energy intake needs; it also provides certain health benefits. This review chapter provides state-of-the-art information to researchers and clinicians who are interested in developing research studies or interv...
Quantum Monte Carlo calculation of the binding energy of the beryllium dimer
NASA Astrophysics Data System (ADS)
Deible, Michael J.; Kessler, Melody; Gasperich, Kevin E.; Jordan, Kenneth D.
2015-08-01
The accurate calculation of the binding energy of the beryllium dimer is a challenging theoretical problem. In this study, the binding energy of Be2 is calculated using the diffusion Monte Carlo (DMC) method, using single Slater determinant and multiconfigurational trial functions. DMC calculations using single-determinant trial wave functions of orbitals obtained from density functional theory calculations overestimate the binding energy, while DMC calculations using Hartree-Fock or CAS(4,8), complete active space trial functions significantly underestimate the binding energy. In order to obtain an accurate value of the binding energy of Be2 from DMC calculations, it is necessary to employ trial functions that include excitations outside the valence space. Our best estimate DMC result for the binding energy of Be2, obtained by using configuration interaction trial functions and extrapolating in the threshold for the configurations retained in the trial function, is 908 cm-1, only slightly below the 935 cm-1 value derived from experiment.
Guidelines for the analysis of free energy calculations.
Klimovich, Pavel V; Shirts, Michael R; Mobley, David L
2015-05-01
Free energy calculations based on molecular dynamics simulations show considerable promise for applications ranging from drug discovery to prediction of physical properties and structure-function studies. But these calculations are still difficult and tedious to analyze, and best practices for analysis are not well defined or propagated. Essentially, each group analyzing these calculations needs to decide how to conduct the analysis and, usually, develop its own analysis tools. Here, we review and recommend best practices for analysis yielding reliable free energies from molecular simulations. Additionally, we provide a Python tool, alchemical-analysis.py, freely available on GitHub as part of the pymbar package (located at http://github.com/choderalab/pymbar), that implements the analysis practices reviewed here for several reference simulation packages, which can be adapted to handle data from other packages. Both this review and the tool covers analysis of alchemical calculations generally, including free energy estimates via both thermodynamic integration and free energy perturbation-based estimators. Our Python tool also handles output from multiple types of free energy calculations, including expanded ensemble and Hamiltonian replica exchange, as well as standard fixed ensemble calculations. We also survey a range of statistical and graphical ways of assessing the quality of the data and free energy estimates, and provide prototypes of these in our tool. We hope this tool and discussion will serve as a foundation for more standardization of and agreement on best practices for analysis of free energy calculations. PMID:25808134
Guidelines for the analysis of free energy calculations
Klimovich, Pavel V.; Shirts, Michael R.; Mobley, David L.
2015-01-01
Free energy calculations based on molecular dynamics (MD) simulations show considerable promise for applications ranging from drug discovery to prediction of physical properties and structure-function studies. But these calculations are still difficult and tedious to analyze, and best practices for analysis are not well defined or propagated. Essentially, each group analyzing these calculations needs to decide how to conduct the analysis and, usually, develop its own analysis tools. Here, we review and recommend best practices for analysis yielding reliable free energies from molecular simulations. Additionally, we provide a Python tool, alchemical–analysis.py, freely available on GitHub at https://github.com/choderalab/pymbar–examples, that implements the analysis practices reviewed here for several reference simulation packages, which can be adapted to handle data from other packages. Both this review and the tool covers analysis of alchemical calculations generally, including free energy estimates via both thermodynamic integration and free energy perturbation-based estimators. Our Python tool also handles output from multiple types of free energy calculations, including expanded ensemble and Hamiltonian replica exchange, as well as standard fixed ensemble calculations. We also survey a range of statistical and graphical ways of assessing the quality of the data and free energy estimates, and provide prototypes of these in our tool. We hope these tools and discussion will serve as a foundation for more standardization of and agreement on best practices for analysis of free energy calculations. PMID:25808134
Spectroscopically Accurate Calculations of the Rovibrational Energies of Diatomic Hydrogen
NASA Astrophysics Data System (ADS)
Perry, Jason
2005-05-01
The Born-Oppenheimer approximation has been used to calculate the rotational and vibrational states of diatomic hydrogen. Because it is an approximation, our group now wants to use a Born-Oppenheimer potential to calculate the electronic energy that has been corrected to match closely with spectroscopic results. We are using a code that has corrections for adiabatic, relativistic, radiative, and non-adiabatic effects. The rovibrational energies have now been calculated for both bound and quasi-bound states. We also want to compute quadrupole transition probabilities for diatomic hydrogen. These calculations aspire to investigate diatomic hydrogen in astrophysical environments.
18 CFR 11.13 - Energy gains calculations.
Code of Federal Regulations, 2011 CFR
2011-04-01
... 18 Conservation of Power and Water Resources 1 2011-04-01 2011-04-01 false Energy gains calculations. 11.13 Section 11.13 Conservation of Power and Water Resources FEDERAL ENERGY REGULATORY COMMISSION, DEPARTMENT OF ENERGY REGULATIONS UNDER THE FEDERAL POWER ACT ANNUAL CHARGES UNDER PART I OF...
18 CFR 11.13 - Energy gains calculations.
Code of Federal Regulations, 2014 CFR
2014-04-01
... 18 Conservation of Power and Water Resources 1 2014-04-01 2014-04-01 false Energy gains calculations. 11.13 Section 11.13 Conservation of Power and Water Resources FEDERAL ENERGY REGULATORY COMMISSION, DEPARTMENT OF ENERGY REGULATIONS UNDER THE FEDERAL POWER ACT ANNUAL CHARGES UNDER PART I OF...
18 CFR 11.13 - Energy gains calculations.
Code of Federal Regulations, 2013 CFR
2013-04-01
... 18 Conservation of Power and Water Resources 1 2013-04-01 2013-04-01 false Energy gains calculations. 11.13 Section 11.13 Conservation of Power and Water Resources FEDERAL ENERGY REGULATORY COMMISSION, DEPARTMENT OF ENERGY REGULATIONS UNDER THE FEDERAL POWER ACT ANNUAL CHARGES UNDER PART I OF...
18 CFR 11.13 - Energy gains calculations.
Code of Federal Regulations, 2012 CFR
2012-04-01
... 18 Conservation of Power and Water Resources 1 2012-04-01 2012-04-01 false Energy gains calculations. 11.13 Section 11.13 Conservation of Power and Water Resources FEDERAL ENERGY REGULATORY COMMISSION, DEPARTMENT OF ENERGY REGULATIONS UNDER THE FEDERAL POWER ACT ANNUAL CHARGES UNDER PART I OF...
WAPA Daily Energy Accounting Activities
Energy Science and Technology Software Center (ESTSC)
1990-10-01
ISA (Interchange, Scheduling, & Accounting) is the interchange scheduling system used by the DOE Western Area Power Administration to perform energy accounting functions associated with the daily activities of the Watertown Operations Office (WOO). The system's primary role is to provide accounting functions for scheduled energy which is exchanged with other power companies and power operating organizations. The system has a secondary role of providing a historical record of all scheduled interchange transactions. The followingmore » major functions are performed by ISA: scheduled energy accounting for received and delivered energy; generation scheduling accounting for both fossil and hydro-electric power plants; metered energy accounting for received and delivered totals; energy accounting for Direct Current (D.C.) Ties; regulation accounting; automatic generation control set calculations; accounting summaries for Basin, Heartland Consumers Power District, and the Missouri Basin Municipal Power Agency; calculation of estimated generation for the Laramie River Station plant; daily and monthly reports; and dual control areas.« less
Science Activities in Energy: Chemical Energy.
ERIC Educational Resources Information Center
Oak Ridge Associated Universities, TN.
Presented is a science activities in energy package which includes 15 activities relating to chemical energy. Activities are simple, concrete experiments for fourth, fifth and sixth grades which illustrate principles and problems relating to energy. Each activity is outlined on a single card which is introduced by a question. A teacher's…
Science Activities in Energy: Electrical Energy.
ERIC Educational Resources Information Center
Oak Ridge Associated Universities, TN.
Presented is a science activities in energy package which includes 16 activities relating to electrical energy. Activities are simple, concrete experiments for fourth, fifth and sixth grades which illustrate principles and problems relating to energy. Each activity is outlined in a single card which is introduced by a question. A teacher's…
Science Activities in Energy: Solar Energy.
ERIC Educational Resources Information Center
Oak Ridge Associated Universities, TN.
Presented is a science activities in energy package which includes 12 activities relating to solar energy. Activities are simple, concrete experiments for fourth, fifth, and sixth grades, which illustrate principles and problems relating to energy. Each activity is outlined on a single card which is introduced by a question. A teacher's supplement…
New approach to calculating the potential energy of colliding nuclei
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.
Calculation of Rydberg energy levels for the francium atom
NASA Astrophysics Data System (ADS)
Huang, Shi-Zhong; Chu, Jin-Min
2010-06-01
Based on the weakest bound electron potential model theory, the Rydberg energy levels and quantum defects of the np2Po1/2 (n = 7-50) and np2Po3/2 (n = 7-50) spectrum series for the francium atom are calculated. The calculated results are in excellent agreement with the 48 measured levels, and 40 energy levels for highly excited states are predicted.
Protein thermostability calculations using alchemical free energy simulations.
Seeliger, Daniel; de Groot, Bert L
2010-05-19
Thermal stability of proteins is crucial for both biotechnological and therapeutic applications. Rational protein engineering therefore frequently aims at increasing thermal stability by introducing stabilizing mutations. The accurate prediction of the thermodynamic consequences caused by mutations, however, is highly challenging as thermal stability changes are caused by alterations in the free energy of folding. Growing computational power, however, increasingly allows us to use alchemical free energy simulations, such as free energy perturbation or thermodynamic integration, to calculate free energy differences with relatively high accuracy. In this article, we present an automated protocol for setting up alchemical free energy calculations for mutations of naturally occurring amino acids (except for proline) that allows an unprecedented, automated screening of large mutant libraries. To validate the developed protocol, we calculated thermodynamic stability differences for 109 mutations in the microbial Ribonuclease Barnase. The obtained quantitative agreement with experimental data illustrates the potential of the approach in protein engineering and design. PMID:20483340
Changing Conceptions of Activation Energy.
ERIC Educational Resources Information Center
Pacey, Philip D.
1981-01-01
Provides background material which relates to the concept of activation energy, fundamental in the study of chemical kinetics. Compares the related concepts of the Arrhenius activation energy, the activation energy at absolute zero, the enthalpy of activation, and the threshold energy. (CS)
Is ring breaking feasible in relative binding free energy calculations?
Liu, Shuai; Wang, Lingle; Mobley, David L
2015-04-27
Our interest is relative binding free energy (RBFE) calculations based on molecular simulations. These are promising tools for lead optimization in drug discovery, computing changes in binding free energy due to modifications of a lead compound. However, in the "alchemical" framework for RBFE calculations, some types of mutations have the potential to introduce error into the computed binding free energies. Here we explore the magnitude of this error in several different model binding calculations. We find that some of the calculations involving ring breaking have significant errors, and these errors are especially large in bridged ring systems. Since the error is a function of ligand strain, which is unpredictable in advance, we believe that ring breaking should be avoided when possible. PMID:25835054
Free energy calculations for molecular solids using GROMACS
NASA Astrophysics Data System (ADS)
Aragones, J. L.; Noya, E. G.; Valeriani, C.; Vega, C.
2013-07-01
In this work, we describe a procedure to evaluate the free energy of molecular solids with the GROMACS molecular dynamics package. The free energy is calculated using the Einstein molecule method that can be regarded as a small modification of the Einstein crystal method. Here, the position and orientation of the molecules is fixed by using an Einstein field that binds with harmonic springs at least three non-collinear atoms (or points of the molecule) to their reference positions. The validity of the Einstein field is tested by performing free-energy calculations of methanol, water (ice), and patchy colloids molecular solids. The free energies calculated with GROMACS show a very good agreement with those obtained using Monte Carlo and with previously published results.
3DRISM Multigrid Algorithm for Fast Solvation Free Energy Calculations.
Sergiievskyi, Volodymyr P; Fedorov, Maxim V
2012-06-12
In this paper we present a fast and accurate method for modeling solvation properties of organic molecules in water with a main focus on predicting solvation (hydration) free energies of small organic compounds. The method is based on a combination of (i) a molecular theory, three-dimensional reference interaction sites model (3DRISM); (ii) a fast multigrid algorithm for solving the high-dimensional 3DRISM integral equations; and (iii) a recently introduced universal correction (UC) for the 3DRISM solvation free energies by properly scaled molecular partial volume (3DRISM-UC, Palmer et al., J. Phys.: Condens. Matter2010, 22, 492101). A fast multigrid algorithm is the core of the method because it helps to reduce the high computational costs associated with solving the 3DRISM equations. To facilitate future applications of the method, we performed benchmarking of the algorithm on a set of several model solutes in order to find optimal grid parameters and to test the performance and accuracy of the algorithm. We have shown that the proposed new multigrid algorithm is on average 24 times faster than the simple Picard method and at least 3.5 times faster than the MDIIS method which is currently actively used by the 3DRISM community (e.g., the MDIIS method has been recently implemented in a new 3DRISM implicit solvent routine in the recent release of the AmberTools 1.4 molecular modeling package (Luchko et al. J. Chem. Theory Comput. 2010, 6, 607-624). Then we have benchmarked the multigrid algorithm with chosen optimal parameters on a set of 99 organic compounds. We show that average computational time required for one 3DRISM calculation is 3.5 min per a small organic molecule (10-20 atoms) on a standard personal computer. We also benchmarked predicted solvation free energy values for all of the compounds in the set against the corresponding experimental data. We show that by using the proposed multigrid algorithm and the 3DRISM-UC model, it is possible to obtain good
S-matrix calculations of energy levels of alkalilike ions
NASA Astrophysics Data System (ADS)
Sapirstein, Jonathan; Cheng, K. T.
2013-05-01
A recent S-matrix based QED calculation of energy levels of the lithium isoelectronic sequence is extended to the general case of a valence electron outside an arbitrary filled core. Formulas are presented that allow calculation of the energy levels of valence ns , np1 / 2 , np3 / 2 , nd3 / 2 , and nd5 / 2 states. Emphasis is placed on modifications of the lithiumlike formulas required because more than one core state is present, and a discussion of an unusual feature of the two-photon exchange contribution involving autoiononizing states is given. The method is illustrated with a calculation of energy levels of the sodium isoelectronic sequence, with results for 3s1 / 2 , 3p1 / 2 , and 3p3 / 2 energies tabulated for the range Z = 20 - 100 . A detailed breakdown of the calculation is given for Z = 74 . Comparison with experiment and other calculations is given, and prospects for extension of the method to ions with more complex electronic structure discussed. The work of JS was supported in part by NSF Grant No. PHY-1068065. The work of KTC was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.
Performance calculation and simulation system of high energy laser weapon
NASA Astrophysics Data System (ADS)
Wang, Pei; Liu, Min; Su, Yu; Zhang, Ke
2014-12-01
High energy laser weapons are ready for some of today's most challenging military applications. Based on the analysis of the main tactical/technical index and combating process of high energy laser weapon, a performance calculation and simulation system of high energy laser weapon was established. Firstly, the index decomposition and workflow of high energy laser weapon was proposed. The entire system was composed of six parts, including classical target, platform of laser weapon, detect sensor, tracking and pointing control, laser atmosphere propagation and damage assessment module. Then, the index calculation modules were designed. Finally, anti-missile interception simulation was performed. The system can provide reference and basis for the analysis and evaluation of high energy laser weapon efficiency.
Calculating fusion neutron energy spectra from arbitrary reactant distributions
NASA Astrophysics Data System (ADS)
Eriksson, J.; Conroy, S.; Andersson Sundén, E.; Hellesen, C.
2016-02-01
The Directional Relativistic Spectrum Simulator (DRESS) code can perform Monte-Carlo calculations of reaction product spectra from arbitrary reactant distributions, using fully relativistic kinematics. The code is set up to calculate energy spectra from neutrons and alpha particles produced in the D(d, n)3He and T(d, n)4He fusion reactions, but any two-body reaction can be simulated by including the corresponding cross section. The code has been thoroughly tested. The kinematics calculations have been benchmarked against the kinematics module of the ROOT Data Analysis Framework. Calculated neutron energy spectra have been validated against tabulated fusion reactivities and against an exact analytical expression for the thermonuclear fusion neutron spectrum, with good agreement. The DRESS code will be used as the core of a detailed synthetic diagnostic framework for neutron measurements at the JET and MAST tokamaks.
Energy Adventure Center. Activity Book.
ERIC Educational Resources Information Center
Carlton, Linda L.
Energy activities are provided in this student activity book. They include: (1) an energy walk; (2) forms of energy in the home; (3) energy conversion; (4) constructing a solar hot dog cooker (with instructions for drawing a parabola); (5) interviewing senior citizens to learn about energy use in the past; (6) packaging materials; (7) insulation;…
Activities Handbook for Energy Education.
ERIC Educational Resources Information Center
DeVito, Alfred; Krockover, Gerald H.
The purpose of this handbook is to present information about energy and to translate this information into learning activities for children. Chapter 1, "Energy: A Delicate Dilemma," presents activities intended to provide an introduction to energy and energy usage. Chapter 2, "What are the Sources of Energy?" provides background information and…
Science Activities in Energy: Conservation.
ERIC Educational Resources Information Center
Oak Ridge Associated Universities, TN.
Presented is a science activities in energy package which includes 14 activities relating to energy conservation. Activities are simple, concrete experiments for fourth, fifth and sixth grades, which illustrate principles and problems relating to energy. Each activity is outlined on a simple card which is introduced by a question. A teacher's…
Lead Optimization Mapper: Automating free energy calculations for lead optimization
Liu, Shuai; Wu, Yujie; Lin, Teng; Abel, Robert; Redmann, Jonathan P.; Summa, Christopher M.; Jaber, Vivian R.; Lim, Nathan M.; Mobley, David L.
2013-01-01
Alchemical free energy calculations hold increasing promise as an aid to drug discovery efforts. However, applications of these techniques in discovery projects have been relatively few, partly because of the difficulty of planning and setting up calculations. Here, we introduce Lead Optimization Mapper, LOMAP, an automated algorithm to plan efficient relative free energy calculations between potential ligands within a substantial library of perhaps hundreds of compounds. In this approach, ligands are first grouped by structural similarity primarily based on the size of a (loosely defined) maximal common substructure, and then calculations are planned within and between sets of structurally related compounds. An emphasis is placed on ensuring that relative free energies can be obtained between any pair of compounds without combining the results of too many different relative free energy calculations (to avoid accumulation of error) and by providing some redundancy to allow for the possibility of error and consistency checking and provide some insight into when results can be expected to be unreliable. The algorithm is discussed in detail and a Python implementation, based on both Schrödinger's and OpenEye's APIs, has been made available freely under the BSD license. PMID:24072356
Calculated viscosity-distance dependence for some actively flowing lavas
NASA Technical Reports Server (NTRS)
Pieri, David
1987-01-01
The importance of viscosity as a gauge of the various energy and momentum dissipation regimes of lava flows has been realized for a long time. Nevertheless, despite its central role in lava dynamics and kinematics, it remains among the most difficult of flow physical properties to measure in situ during an eruption. Attempts at reconstructing the actual emplacement viscosities of lava flows from their solidified topographic form are difficult. Where data are available on the position of an advancing flow front as a function of time, it is possible to calculate the effective viscosity of the front as a function of distance from the vent, under the assumptions of a steady state regime. As an application and test of an equation given, relevant parameters from five recent flows on Mauna Loa and Kilauea were utilized to infer the dynamic structure of their aggregate flow front viscosity as they advanced, up to cessation. The observed form of the viscosity-distance relation for the five active Hawaiian flows examined appears to be exponential, with a rapid increase just before the flows stopped as one would expect.
Total-energy and pressure calculations for random substitutional alloys
Johnson, D.D. ); Nicholson, D.M. ); Pinski, F.J. ); Gyoerffy, B.L. ); Stocks, G.M. )
1990-05-15
We present the details and the derivation of density-functional-based expressions for the total energy and pressure for random substitutional alloys (RSA) using the Korringa-Kohn-Rostoker Green's-function approach in combination with the coherent-potential approximation (CPA) to treat the configurational averaging. This includes algebraic cancellation of various electronic core contributions to the total energy and pressure, as in ordered-solid muffin-tin-potential calculations. Thus, within the CPA, total-energy and pressure calculations for RSA have the same foundation and have been found to have the same accuracy as those obtained in similar calculations for ordered solids. Results of our calculations for the impurity formation energy, and for the bulk moduli, the lattice parameters, and the energy of mixing as a function of concentration in fcc Cu{sub {ital c}}Zn{sub 1{minus}{ital c}} alloys show that this generalized density-functional theory will be useful in studying alloy phase stability.
Magnetic field calculation and measurement of active magnetic bearings
NASA Astrophysics Data System (ADS)
Ding, Guoping; Zhou, Zude; Hu, Yefa
2006-11-01
Magnetic Bearings are typical devices in which electric energy and mechanical energy convert mutually. Magnetic Field indicates the relationship between 2 of the most important parameters in a magnetic bearing - current and force. This paper presents calculation and measurement of the magnetic field distribution of a self-designed magnetic bearing. Firstly, the static Maxwell's equations of the magnetic bearing are presented and a Finite Element Analysis (FEA) is found to solve the equations and get post-process results by means of ANSYS software. Secondly, to confirm the calculation results a Lakeshore460 3-channel Gaussmeter is used to measure the magnetic flux density of the magnetic bearing in X, Y, Z directions accurately. According to the measurement data the author constructs a 3D magnetic field distribution digital model by means of MATLAB software. Thirdly, the calculation results and the measurement data are compared and analyzed; the comparing result indicates that the calculation results are consistent with the measurement data in allowable dimension variation, which means that the FEA calculation method of the magnetic bearing has high precision. Finally, it is concluded that the magnetic field calculation and measurement can accurately reflect the real magnetic distribution in the magnetic bearing and the result can guide the design and analysis of the magnetic bearing effectively.
Ion beam energy spectrum calculation via dosimetry data deconvolution.
Harper-Slaboszewicz, Victor Jozef; Sharp, Andrew Clinton
2010-10-01
The energy spectrum of a H{sup +} beam generated within the HERMES III accelerator is calculated from dosimetry data to refine future experiments. Multiple layers of radiochromic film are exposed to the beam. A graphic user interface was written in MATLAB to align the film images and calculate the beam's dose depth profile. Singular value regularization is used to stabilize the unfolding and provide the H{sup +} beam's energy spectrum. The beam was found to have major contributions from 1 MeV and 8.5 MeV protons. The HERMES III accelerator is typically used as a pulsed photon source to experimentally obtain photon impulse response of systems due to high energy photons. A series of experiments were performed to explore the use of Hermes III to generate an intense pulsed proton beam. Knowing the beam energy spectrum allows for greater precision in experiment predictions and beam model verification.
Calculation of exchange energies using algebraic perturbation theory
Burrows, B. L.; Dalgarno, A.; Cohen, M.
2010-04-15
An algebraic perturbation theory is presented for efficient calculations of localized states and hence of exchange energies, which are the differences between low-lying states of the valence electron of a molecule, formed by the collision of an ion Y{sup +} with an atom X. For the case of a homonuclear molecule these are the gerade and ungerade states and the exchange energy is an exponentially decreasing function of the internuclear distance. For such homonuclear systems the theory is used in conjunction with the Herring-Holstein technique to give accurate exchange energies for a range of intermolecular separations R. Since the perturbation parameter is essentially 1/R, this method is suitable for large R. In particular, exchange energies are calculated for X{sub 2}{sup +} systems, where X is H, Li, Na, K, Rb, or Cs.
Science Activities in Energy: Wind Energy.
ERIC Educational Resources Information Center
Oak Ridge Associated Universities, TN.
Included in this science activities energy package are 12 activities related to wind energy for elementary students. Each activity is outlined on a single card and is introduced by a question. Topics include: (1) At what time of day is there enough wind to make electricity where you live?; (2) Where is the windiest spot on your schoolground?; and…
On the calculation of classical vibrational energy exchange
NASA Astrophysics Data System (ADS)
Gibbons, John P.; Stettler, John D.
1982-07-01
A three-dimensional, Monte Carlo classical model for the calculation of vibrational energy relaxation and transfer rates for both diatomic—monatomic and diatomic—diatomic systems was developed, analyzed and implemented. Mediation by internal angular momentum changes was demonstrated to be important in these energy transfer processes. This mechanism was incorporated into the model in order to achieve statistically significant results within reasonable computer running times. This made possible the extension of the model calculations to much lower temperatures than had been previously investigated. This calculational procedure was applied to Ar—O 2, to He—O 2 and to the near resonant CO—N 2 process at several temperatures between room temperature and 4000 K with the use of exponential repulsive intermolecular potential. Three different sets of potential parameters obtained from three independent sources were used. The results were compared to experiment.
The Calculation of Accurate Metal-Ligand Bond Energies
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W.; Partridge, Harry, III; Ricca, Alessandra; Arnold, James O. (Technical Monitor)
1997-01-01
The optimization of the geometry and calculation of zero-point energies are carried out at the B3LYP level of theory. The bond energies are determined at this level, as well as at the CCSD(T) level using very large basis sets. The successive OH bond energies to the first row transition metal cations are reported. For most systems there has been an experimental determination of the first OH. In general, the CCSD(T) values are in good agreement with experiment. The bonding changes from mostly covalent for the early metals to mostly electrostatic for the late transition metal systems.
Accurate calculation of diffraction-limited encircled and ensquared energy.
Andersen, Torben B
2015-09-01
Mathematical properties of the encircled and ensquared energy functions for the diffraction-limited point-spread function (PSF) are presented. These include power series and a set of linear differential equations that facilitate the accurate calculation of these functions. Asymptotic expressions are derived that provide very accurate estimates for the relative amount of energy in the diffraction PSF that fall outside a square or rectangular large detector. Tables with accurate values of the encircled and ensquared energy functions are also presented. PMID:26368873
Calculation of Mg(+)-ligand relative binding energies
NASA Technical Reports Server (NTRS)
Partridge, Harry; Bauschlicher, Charles W., Jr.
1992-01-01
The calculated relative binding energies of 16 organic molecules to Mg(+) are compared with experimental results where available. The geometries of the ligands and the Mg(+)-ligand complexes arc optimized at the self-consistent field level using a 6-31G* basis set. The Mg(+) binding energies are evaluated using second-order perturbation theory and basis sets of triple-sigma quality augmented with two sets of polarization functions. This level of theory is calibrated against higher levels of theory for selected systems. The computed binding energies are accurate to about 2 kcal/mol.
Gauthier, J.H.; Zieman, N.B.; Miller, W.B.
1991-10-01
The purpose of the the Code Verification (COVE) 2A benchmarking activity is to assess the numerical accuracy of several computer programs for the Yucca Mountain Site Characterization Project of the Department of Energy. This paper presents a brief description of the computer program TOSPAC and a discussion of the calculational effort and results generated by TOSPAC for the COVE 2A problem set. The calculations were performed twice. The initial calculations provided preliminary results for comparison with the results from other COVE 2A participants. TOSPAC was modified in response to the comparison and the final calculations included a correction and several enhancements to improve efficiency. 8 refs.
Three dimensional calculation of flux of low energy atmospheric neutrinos
NASA Technical Reports Server (NTRS)
Lee, H.; Bludman, S. A.
1985-01-01
Results of three-dimensional Monte Carlo calculation of low energy flux of atmospheric neutrinos are presented and compared with earlier one-dimensional calculations 1,2 valid at higher neutrino energies. These low energy neutrinos are the atmospheric background in searching for neutrinos from astrophysical sources. Primary cosmic rays produce the neutrino flux peaking at near E sub=40 MeV and neutrino intensity peaking near E sub v=100 MeV. Because such neutrinos typically deviate by 20 approximately 30 from the primary cosmic ray direction, three-dimensional effects are important for the search of atmospheric neutrinos. Nevertheless, the background of these atmospheric neutrinos is negligible for the detection of solar and supernova neutrinos.
SNS Sample Activation Calculator Flux Recommendations and Validation
McClanahan, Tucker C.; Gallmeier, Franz X.; Iverson, Erik B.; Lu, Wei
2015-02-01
The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) uses the Sample Activation Calculator (SAC) to calculate the activation of a sample after the sample has been exposed to the neutron beam in one of the SNS beamlines. The SAC webpage takes user inputs (choice of beamline, the mass, composition and area of the sample, irradiation time, decay time, etc.) and calculates the activation for the sample. In recent years, the SAC has been incorporated into the user proposal and sample handling process, and instrument teams and users have noticed discrepancies in the predicted activation of their samples. The Neutronics Analysis Team validated SAC by performing measurements on select beamlines and confirmed the discrepancies seen by the instrument teams and users. The conclusions were that the discrepancies were a result of a combination of faulty neutron flux spectra for the instruments, improper inputs supplied by SAC (1.12), and a mishandling of cross section data in the Sample Activation Program for Easy Use (SAPEU) (1.1.2). This report focuses on the conclusion that the SAPEU (1.1.2) beamline neutron flux spectra have errors and are a significant contributor to the activation discrepancies. The results of the analysis of the SAPEU (1.1.2) flux spectra for all beamlines will be discussed in detail. The recommendations for the implementation of improved neutron flux spectra in SAPEU (1.1.3) are also discussed.
Composite electron propagator methods for calculating ionization energies
NASA Astrophysics Data System (ADS)
Díaz-Tinoco, Manuel; Dolgounitcheva, O.; Zakrzewski, V. G.; Ortiz, J. V.
2016-06-01
Accurate ionization energies of molecules may be determined efficiently with composite electron-propagator (CEP) techniques. These methods estimate the results of a calculation with an advanced correlation method and a large basis set by performing a series of more tractable calculations in which large basis sets are used with simpler approximations and small basis sets are paired with more demanding correlation techniques. The performance of several CEP methods, in which diagonal, second-order electron propagator results with large basis sets are combined with higher-order results obtained with smaller basis sets, has been tested for the ionization energies of closed-shell molecules from the G2 set. Useful compromises of accuracy and computational efficiency employ complete-basis-set extrapolation for second-order results and small basis sets in third-order, partial third-order, renormalized partial-third order, or outer valence Green's function calculations. Analysis of results for vertical as well as adiabatic ionization energies leads to specific recommendations on the best use of regular and composite methods. Results for 22 organic molecules of interest in the design of photovoltaic devices, benzo[a]pyrene, Mg-octaethylporphyrin, and C60 illustrate the capabilities of CEP methods for calculations on large molecules.
Calculating Free Energy Changes in Continuum Solvation Models.
Ho, Junming; Ertem, Mehmed Z
2016-02-25
We recently showed for a large data set of pKas and reduction potentials that free energies calculated directly within the SMD continuum model compares very well with corresponding thermodynamic cycle calculations in both aqueous and organic solvents [ Phys. Chem. Chem. Phys. 2015 , 17 , 2859 ]. In this paper, we significantly expand the scope of our study to examine the suitability of this approach for calculating general solution phase kinetics and thermodynamics, in conjunction with several commonly used solvation models (SMD-M062X, SMD-HF, CPCM-UAKS, and CPCM-UAHF) for a broad range of systems. This includes cluster-continuum schemes for pKa calculations as well as various neutral, radical, and ionic reactions such as enolization, cycloaddition, hydrogen and chlorine atom transfer, and SN2 and E2 reactions. On the basis of this benchmarking study, we conclude that the accuracies of both approaches are generally very similar-the mean errors for Gibbs free energy changes of neutral and ionic reactions are approximately 5 and 25 kJ mol(-1), respectively. In systems where there are significant structural changes due to solvation, as is the case for certain ionic transition states and amino acids, the direct approach generally afford free energy changes that are in better agreement with experiment. PMID:26878566
Composite electron propagator methods for calculating ionization energies.
Díaz-Tinoco, Manuel; Dolgounitcheva, O; Zakrzewski, V G; Ortiz, J V
2016-06-14
Accurate ionization energies of molecules may be determined efficiently with composite electron-propagator (CEP) techniques. These methods estimate the results of a calculation with an advanced correlation method and a large basis set by performing a series of more tractable calculations in which large basis sets are used with simpler approximations and small basis sets are paired with more demanding correlation techniques. The performance of several CEP methods, in which diagonal, second-order electron propagator results with large basis sets are combined with higher-order results obtained with smaller basis sets, has been tested for the ionization energies of closed-shell molecules from the G2 set. Useful compromises of accuracy and computational efficiency employ complete-basis-set extrapolation for second-order results and small basis sets in third-order, partial third-order, renormalized partial-third order, or outer valence Green's function calculations. Analysis of results for vertical as well as adiabatic ionization energies leads to specific recommendations on the best use of regular and composite methods. Results for 22 organic molecules of interest in the design of photovoltaic devices, benzo[a]pyrene, Mg-octaethylporphyrin, and C60 illustrate the capabilities of CEP methods for calculations on large molecules. PMID:27305999
Binding Energy Calculations for Novel Ternary Ionic Lattices
NASA Astrophysics Data System (ADS)
Rodríguez-Mijangos, Ricardo; Vazquez-Polo, Gustavo
2002-03-01
Theoretical calculations for the binding energy between metalic ions and negative ions on a novel ternary ionic lattice is carried out for several solid solutions prepared with different concentrations and characterized recently (1). The ternary lattices that reach a good miscibility are: KCl(x)KBr(y)RbCl(z) in three different concentrations: (x=y=z=0.33), (x=0.5, y=0.25, z=0.25) and (x=0.33, y=0.07, z=0.60). The binding energy for these novel structures is calculated from the lattice constants obtained by X ray diffractometry analysis performed on the samples and the Vegard law (2). For the repulsive force exponent m, an average of the m values was considered. The energy values obtained by the Born´expression are compared with corresponding energy values from the lattice with more complex expressions, such as the Born Mayer, Born-Van der Walls. There is a good aggreement between all these calculations. (1)R. R. Mijangos, A. Cordero-Borboa, E. Alvarez, M. Cervantes, Physics Letters A 282 (2001) 195-200. (2) G. Vazquez-Polo, R. R. Mijangos et al. Revista Mexicana de Fisica, 47, Diciembre 2001. In Press.
Science Activities in Energy: Solar Energy II.
ERIC Educational Resources Information Center
Oak Ridge Associated Universities, TN.
Included in this science activities energy package are 14 activities related to solar energy for secondary students. Each activity is outlined on a single card and is introduced by a question such as: (1) how much solar heat comes from the sun? or (2) how many times do you have to run water through a flat-plate collector to get a 10 degree rise in…
Self-consistent calculations of alpha-decay energies
Tolokonnikov, S. V.; Lutostansky, Yu. S.; Saperstein, E. E.
2013-06-15
On the basis of the self-consistent theory of finite Fermi systems, the energies of alphadecay chains were calculated for several new superheavy nuclei discovered recently in experiments of the Dubna-Livermore Collaboration headed by Yu.Ts. Oganessian. The approach in question is implemented on the basis of the generalized method of the density functional proposed by Fayans and his coauthors. The version used here relies on the functional DF3-a proposed recently for describing a wide array of nuclear data, including data on superheavy nuclei. A detailed comparison of the results obtained on this basis with the predictions of different approaches, including the self-consistent Skyrme-Hartree-Fock method, the micro-macro method in the version developed by Sobiczewski and his coauthors, and the phenomenological method of Liran and his coauthors, is performed. The resulting alpha-decay energies are used to calculate respective lifetimes with the aid of the phenomenological Parkhomenko-Sobiczewski formula.
Improved initial guess for minimum energy path calculations.
Smidstrup, Søren; Pedersen, Andreas; Stokbro, Kurt; Jónsson, Hannes
2014-06-01
A method is presented for generating a good initial guess of a transition path between given initial and final states of a system without evaluation of the energy. An objective function surface is constructed using an interpolation of pairwise distances at each discretization point along the path and the nudged elastic band method then used to find an optimal path on this image dependent pair potential (IDPP) surface. This provides an initial path for the more computationally intensive calculations of a minimum energy path on an energy surface obtained, for example, by ab initio or density functional theory. The optimal path on the IDPP surface is significantly closer to a minimum energy path than a linear interpolation of the Cartesian coordinates and, therefore, reduces the number of iterations needed to reach convergence and averts divergence in the electronic structure calculations when atoms are brought too close to each other in the initial path. The method is illustrated with three examples: (1) rotation of a methyl group in an ethane molecule, (2) an exchange of atoms in an island on a crystal surface, and (3) an exchange of two Si-atoms in amorphous silicon. In all three cases, the computational effort in finding the minimum energy path with DFT was reduced by a factor ranging from 50% to an order of magnitude by using an IDPP path as the initial path. The time required for parallel computations was reduced even more because of load imbalance when linear interpolation of Cartesian coordinates was used. PMID:24907989
Improved initial guess for minimum energy path calculations
Smidstrup, Søren; Pedersen, Andreas; Stokbro, Kurt
2014-06-07
A method is presented for generating a good initial guess of a transition path between given initial and final states of a system without evaluation of the energy. An objective function surface is constructed using an interpolation of pairwise distances at each discretization point along the path and the nudged elastic band method then used to find an optimal path on this image dependent pair potential (IDPP) surface. This provides an initial path for the more computationally intensive calculations of a minimum energy path on an energy surface obtained, for example, by ab initio or density functional theory. The optimal path on the IDPP surface is significantly closer to a minimum energy path than a linear interpolation of the Cartesian coordinates and, therefore, reduces the number of iterations needed to reach convergence and averts divergence in the electronic structure calculations when atoms are brought too close to each other in the initial path. The method is illustrated with three examples: (1) rotation of a methyl group in an ethane molecule, (2) an exchange of atoms in an island on a crystal surface, and (3) an exchange of two Si-atoms in amorphous silicon. In all three cases, the computational effort in finding the minimum energy path with DFT was reduced by a factor ranging from 50% to an order of magnitude by using an IDPP path as the initial path. The time required for parallel computations was reduced even more because of load imbalance when linear interpolation of Cartesian coordinates was used.
Calculating Free Energies Using Scaled-Force Molecular Dynamics Algorithm
NASA Technical Reports Server (NTRS)
Darve, Eric; Wilson, Micahel A.; Pohorille, Andrew
2000-01-01
One common objective of molecular simulations in chemistry and biology is to calculate the free energy difference between different states of the system of interest. Examples of problems that have such an objective are calculations of receptor-ligand or protein-drug interactions, associations of molecules in response to hydrophobic, and electrostatic interactions or partition of molecules between immiscible liquids. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), 'native' state. Perhaps the best example of such a problem is folding of proteins or short RNA molecules. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to 'quasi non-ergodicity', whereby a part of phase space is inaccessible on timescales of the simulation. A host of strategies has been developed to improve efficiency of sampling the phase space. For example, some Monte Carlo techniques involve large steps which move the system between low-energy regions in phase space without the need for sampling the configurations corresponding to energy barriers (J-walking). Most strategies, however, rely on modifying probabilities of sampling low and high-energy regions in phase space such that transitions between states of interest are encouraged. Perhaps the simplest implementation of this strategy is to increase the temperature of the system. This approach was successfully used to identify denaturation pathways in several proteins, but it is clearly not applicable to protein folding. It is also not a successful method for determining free energy differences. Finally, the approach is likely to fail for systems with co-existing phases, such as water-membrane systems, because it may lead to spontaneous
Quantum Monte Carlo calculation of the binding energy of the beryllium dimer
Deible, Michael J.; Kessler, Melody; Gasperich, Kevin E.; Jordan, Kenneth D.
2015-08-28
The accurate calculation of the binding energy of the beryllium dimer is a challenging theoretical problem. In this study, the binding energy of Be{sub 2} is calculated using the diffusion Monte Carlo (DMC) method, using single Slater determinant and multiconfigurational trial functions. DMC calculations using single-determinant trial wave functions of orbitals obtained from density functional theory calculations overestimate the binding energy, while DMC calculations using Hartree-Fock or CAS(4,8), complete active space trial functions significantly underestimate the binding energy. In order to obtain an accurate value of the binding energy of Be{sub 2} from DMC calculations, it is necessary to employ trial functions that include excitations outside the valence space. Our best estimate DMC result for the binding energy of Be{sub 2}, obtained by using configuration interaction trial functions and extrapolating in the threshold for the configurations retained in the trial function, is 908 cm{sup −1}, only slightly below the 935 cm{sup −1} value derived from experiment.
Quantum Monte Carlo calculation of the binding energy of the beryllium dimer.
Deible, Michael J; Kessler, Melody; Gasperich, Kevin E; Jordan, Kenneth D
2015-08-28
The accurate calculation of the binding energy of the beryllium dimer is a challenging theoretical problem. In this study, the binding energy of Be2 is calculated using the diffusion Monte Carlo (DMC) method, using single Slater determinant and multiconfigurational trial functions. DMC calculations using single-determinant trial wave functions of orbitals obtained from density functional theory calculations overestimate the binding energy, while DMC calculations using Hartree-Fock or CAS(4,8), complete active space trial functions significantly underestimate the binding energy. In order to obtain an accurate value of the binding energy of Be2 from DMC calculations, it is necessary to employ trial functions that include excitations outside the valence space. Our best estimate DMC result for the binding energy of Be2, obtained by using configuration interaction trial functions and extrapolating in the threshold for the configurations retained in the trial function, is 908 cm(-1), only slightly below the 935 cm(-1) value derived from experiment. PMID:26328827
Non-Equilibrium Properties from Equilibrium Free Energy Calculations
NASA Technical Reports Server (NTRS)
Pohorille, Andrew; Wilson, Michael A.
2012-01-01
Calculating free energy in computer simulations is of central importance in statistical mechanics of condensed media and its applications to chemistry and biology not only because it is the most comprehensive and informative quantity that characterizes the eqUilibrium state, but also because it often provides an efficient route to access dynamic and kinetic properties of a system. Most of applications of equilibrium free energy calculations to non-equilibrium processes rely on a description in which a molecule or an ion diffuses in the potential of mean force. In general case this description is a simplification, but it might be satisfactorily accurate in many instances of practical interest. This hypothesis has been tested in the example of the electrodiffusion equation . Conductance of model ion channels has been calculated directly through counting the number of ion crossing events observed during long molecular dynamics simulations and has been compared with the conductance obtained from solving the generalized Nernst-Plank equation. It has been shown that under relatively modest conditions the agreement between these two approaches is excellent, thus demonstrating the assumptions underlying the diffusion equation are fulfilled. Under these conditions the electrodiffusion equation provides an efficient approach to calculating the full voltage-current dependence routinely measured in electrophysiological experiments.
Free energy calculations for a flexible water model.
Habershon, Scott; Manolopoulos, David E
2011-11-28
In this work, we consider the problem of calculating the classical free energies of liquids and solids for molecular models with intramolecular flexibility. We show that thermodynamic integration from the fully-interacting solid of interest to a Debye crystal reference state, with anisotropic harmonic interactions derived from the Hessian of the original crystal, provides a straightforward route to calculating the Gibbs free energy of the solid. To calculate the molecular liquid free energy, it is essential to correctly account for contributions from both intermolecular and intramolecular motion; we employ thermodynamic integration to a Lennard-Jones reference fluid, coupled with direct evaluation of the molecular ro-vibrational partition function. These approaches are used to study the low-pressure classical phase diagram of the flexible q-TIP4P/F water model. We find that, while the experimental ice-I/liquid and ice-III/liquid coexistence lines are described reasonably well by this model, the ice-II phase is predicted to be metastable. In light of this finding, we go on to examine how the coupling between intramolecular flexibility and intermolecular interactions influences the computed phase diagram by comparing our results with those of the underlying rigid-body water model. PMID:21887423
Free energy calculation of permeation through aquaporin-5
NASA Astrophysics Data System (ADS)
Bastien, David
The work of this paper continues upon the large area of research being done on aquaporins (AQPs). AQPs are proteins that take on the role of facilitating the transfer of substances, mainly water, across cell membranes. There are many different types of AQPs, with each of these highly selective proteins conducting only certain solutes, along with unique permeability rates. The permeation characteristics of aquaporins rely mostly on the residue hydrophobicity and steric restraints of the aromatic arginine (ar/R) region of the protein channel. The purpose of this paper is to analyze the structures of aquaporin-5 (AQP5) and aquaglycerolporin (Glpf), including a radius profile of the respective protein channels, and to compare them to permeation events using steered molecular dynamics (SMD) pulling simulations. Two in silico experiments are performed in order to achieve the free Energy landscape of a single water molecule permeating through the four channels of both Aqp5 and GlpF. The equilibrium free energy curves are calculated from the non-equilibrium, irreversible work measurements using the fluctuation-dissipation theorem (FDT) of Brownian dynamicis (BD). The free energy profiles are then compared and related to the structural profiles of AQP5 and GlpF. The change in free energy across the ar/R region in AQP5 is found to be reasonably larger than that of GlpF. The free energy profiles of AQP5 and GlpF agree with the diameter profile of the channels respectively. Furthermore, free energy calculations are computed for the permeation of Na+ and Cl- ions through the central pore of Aqp5, which provide some insight into the structural mechanisms of AQP5. The free energy barrier for ion transport through the central pore is found to be very large, peaking at around 11 Kcal/mol for chloride and 20 Kcal/mol for sodium.
ERIC Educational Resources Information Center
Vargas, Francisco M.
2014-01-01
The temperature dependence of the Gibbs energy and important quantities such as Henry's law constants, activity coefficients, and chemical equilibrium constants is usually calculated by using the Gibbs-Helmholtz equation. Although, this is a well-known approach and traditionally covered as part of any physical chemistry course, the required…
Activities for Teaching Solar Energy.
ERIC Educational Resources Information Center
Mason, Jack Lee; Cantrell, Joseph S.
1980-01-01
Plans and activities are suggested for teaching elementary children about solar energy. Directions are included for constructing a flat plate collector and a solar oven. Activities for a solar field day are given. (SA)
Variational calculations of the HT{sup +} rovibrational energies
Bekbaev, A. K.; Korobov, V. I.; Dineykhan, M.
2011-04-15
In this Brief Report, we use the exponential explicitly correlated variational basis set of the type exp(-{alpha}{sub n}R-{beta}{sub n}r{sub 1}-{gamma}{sub n}r{sub 2}) to calculate systematically the nonrelativistic bound-state energies for the hydrogen molecular ion HT{sup +}. We perform calculations for the states of the total orbital angular momentum L=0 and 1 with the complete set of vibrational quantum numbers v= 0-23, as well as for the states of L= 2-5 and v= 0-5. The E1 dipole transition moments, which are of importance for the planning of spectroscopic laser experiments, have been obtained as well.
Prediction of Protein Solubility from Calculation of Transfer Free Energy
Tjong, Harianto; Zhou, Huan-Xiang
2008-01-01
Solubility plays a major role in protein purification, and has serious implications in many diseases. We studied the effects of pH and mutations on protein solubility by calculating the transfer free energy from the condensed phase to the solution phase. The condensed phase was modeled as an implicit solvent, with a dielectric constant lower than that of water. To account for the effects of pH, the protonation states of titratable side chains were sampled by running constant-pH molecular dynamics simulations. Conformations were then selected for calculations of the electrostatic solvation energy: once for the condensed phase, and once for the solution phase. The average transfer free energy from the condensed phase to the solution phase was found to predict reasonably well the variations in solubility of ribonuclease Sa and insulin with pH. This treatment of electrostatic contributions combined with a similar approach for nonelectrostatic contributions led to a quantitative rationalization of the effects of point mutations on the solubility of ribonuclease Sa. This study provides valuable insights into the physical basis of protein solubility. PMID:18515380
SCALE Continuous-Energy Eigenvalue Sensitivity Coefficient Calculations
Perfetti, Christopher M.; Rearden, Bradley T.; Martin, William R.
2016-02-25
Sensitivity coefficients describe the fractional change in a system response that is induced by changes to system parameters and nuclear data. The Tools for Sensitivity and UNcertainty Analysis Methodology Implementation (TSUNAMI) code within the SCALE code system makes use of eigenvalue sensitivity coefficients for an extensive number of criticality safety applications, including quantifying the data-induced uncertainty in the eigenvalue of critical systems, assessing the neutronic similarity between different critical systems, and guiding nuclear data adjustment studies. The need to model geometrically complex systems with improved fidelity and the desire to extend TSUNAMI analysis to advanced applications has motivated the developmentmore » of a methodology for calculating sensitivity coefficients in continuous-energy (CE) Monte Carlo applications. The Contributon-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) and Iterated Fission Probability (IFP) eigenvalue sensitivity methods were recently implemented in the CE-KENO framework of the SCALE code system to enable TSUNAMI-3D to perform eigenvalue sensitivity calculations using continuous-energy Monte Carlo methods. This work provides a detailed description of the theory behind the CLUTCH method and describes in detail its implementation. This work explores the improvements in eigenvalue sensitivity coefficient accuracy that can be gained through the use of continuous-energy sensitivity methods and also compares several sensitivity methods in terms of computational efficiency and memory requirements.« less
Exploiting the Properties of Aquaporin to Calculate Free Energy
NASA Astrophysics Data System (ADS)
Espejel, Hugo; Chen, Liao
2010-03-01
Aquaporins' (AQPs) main purpose is to facilitate the transfer of water molecules through a molecular membrane. We can calculate the free energy of the AQP system when water permeates through it. This is performed using the Visual Molecular Dynamics (VMD) and the Nanoscale Molecular Dynamics (NAMD) programs. In our first set of experiments, AQP is submerged in a body of water, in which case a water molecule near AQP is pulled through the protein. The data is then used to calculate the free energy using two different equations: the Jarzynski equality and the fluctuation-dissipation theorem. The values from both equations are then compared to examine their accuracy. The second set of experiments has the same set up, but now AQP is embedded in a lipid bilayer. We found that both equations give values that are much smaller than kT. This verifies that AQP is a channel for water molecules because the pulling of water gives constant values of free energy. We also found that the water molecules' negative poles were all pointing towards the center of the AQP channel. This means that the process of proton transport in AQP is overwhelmingly difficult.
NASA Astrophysics Data System (ADS)
Hafiz, Hasnain; Barbiellini, B.; Jia, Q.; Tylus, U.; Strickland, K.; Bansil, A.; Mukerjee, S.
2015-03-01
Catalysts based on Fe/N/C clusters can support the oxygen-reduction reaction (ORR) without the use of expensive metals such as platinum. These systems can also prevent some poisonous species to block the active sites from the reactant. We have performed spin-polarized calculations on various Fe/N/C fragments using the Vienna Ab initio Simulation Package (VASP) code. Some results are compared to similar calculations obtained with the Gaussian code. We investigate the partial density of states (PDOS) of the 3d orbitals near the Fermi level and calculate the binding energies of several ligands. Correlations of the binding energies with the 3d electronic PDOS's are used to propose electronic descriptors of the ORR associated with the 3d states of Fe. We also suggest a structural model for the most active site with a ferrous ion (Fe2+) in the high spin state or the so-called Doublet 3 (D3).
Activation calculations for trapped protons below 200 MeV: Appendix
NASA Technical Reports Server (NTRS)
Laird, C. E.
1991-01-01
Tables are given displaying of the results of the activation calculations of metal samples and other material aboard the Long Duration Exposure Facility-1 (LDEF-1) and Spacelab-2 with the computer program, PTRAP4. The computer printouts give the reaction, the reactant product, the proton reaction cross sections as a function of the energy of the incident protons, and the activation as a function of distance into the sample from the exposed surface.
Parquet decomposition calculations of the electronic self-energy
NASA Astrophysics Data System (ADS)
Gunnarsson, O.; Schäfer, T.; LeBlanc, J. P. F.; Merino, J.; Sangiovanni, G.; Rohringer, G.; Toschi, A.
2016-06-01
The parquet decomposition of the self-energy into classes of diagrams, those associated with specific scattering processes, can be exploited for different scopes. In this work, the parquet decomposition is used to unravel the underlying physics of nonperturbative numerical calculations. We show the specific example of dynamical mean field theory and its cluster extensions [dynamical cluster approximation (DCA)] applied to the Hubbard model at half-filling and with hole doping: These techniques allow for a simultaneous determination of two-particle vertex functions and self-energies and, hence, for an essentially "exact" parquet decomposition at the single-site or at the cluster level. Our calculations show that the self-energies in the underdoped regime are dominated by spin-scattering processes, consistent with the conclusions obtained by means of the fluctuation diagnostics approach [O. Gunnarsson et al., Phys. Rev. Lett. 114, 236402 (2015), 10.1103/PhysRevLett.114.236402]. However, differently from the latter approach, the parquet procedure displays important changes with increasing interaction: Even for relatively moderate couplings, well before the Mott transition, singularities appear in different terms, with the notable exception of the predominant spin channel. We explain precisely how these singularities, which partly limit the utility of the parquet decomposition and, more generally, of parquet-based algorithms, are never found in the fluctuation diagnostics procedure. Finally, by a more refined analysis, we link the occurrence of the parquet singularities in our calculations to a progressive suppression of charge fluctuations and the formation of a resonance valence bond state, which are typical hallmarks of a pseudogap state in DCA.
Investigation of Multipole Electrostatics in Hydration Free Energy Calculations
Shi, Yue; Wu, Chuanjie; Ponder, Jay W.; Ren, Pengyu
2010-01-01
Hydration free energy (HFE) is generally used for evaluating molecular solubility, which is an important property for pharmaceutical and chemical engineering processes. Accurately predicting HFE is also recognized as one fundamental capability of molecular mechanics force field. Here we present a systematic investigation on HFE calculations with AMOEBA polarizable force field at various parameterization and simulation conditions. The HFEs of seven small organic molecules have been obtained alchemically using the Bennett Acceptance Ratio (BAR) method. We have compared two approaches to derive the atomic multipoles from quantum mechanical (QM) calculations: one directly from the new distributed multipole analysis (DMA) and the other involving fitting to the electrostatic potential around the molecules. Wave functions solved at the MP2 level with four basis sets (6-311G*, 6-311++G(2d,2p), cc-pVTZ, and aug-cc-pVTZ) are used to derive the atomic multipoles. HFEs from all four basis sets show a reasonable agreement with experimental data (root mean square error 0.63 kcal/mol for aug-ccpVTZ). We conclude that aug-cc-pVTZ gives the best performance when used with AMOEBA, and 6-311++G(2d,2p) is comparable but more efficient for larger systems. The results suggest that the inclusion of diffuse basis functions is important for capturing intermolecular interactions. The effect of long-range correction to van der Waals interaction on the hydration free energies is about 0.1 kcal/mol when the cutoff is 12Å, and increases linearly with the number of atoms in the solute/ligand. In addition, we also discussed the results from a hybrid approach that combines polarizable solute with fixed-charge water in the hydration free energy calculation. PMID:20925089
Empirically corrected HEAT method for calculating atomization energies
Brand, Holmann V
2008-01-01
We describe how to increase the accuracy ofthe most recent variants ofthe HEAT method for calculating atomization energies of molecules by means ofextremely simple empirical corrections that depend on stoichiometry and the number ofunpaired electrons in the molecule. Our corrections reduce the deviation from experiment for all the HEAT variants. In particular, our corrections reduce the average absolute deviation and the root-mean-square deviation ofthe 456-QP variant to 0.18 and 0.23 kJoule/mol (i.e., 0.04 and 0.05 kcallmol), respectively.
Multiple scattering calculations of relativistic electron energy loss spectra
NASA Astrophysics Data System (ADS)
Jorissen, K.; Rehr, J. J.; Verbeeck, J.
2010-04-01
A generalization of the real-space Green’s-function approach is presented for ab initio calculations of relativistic electron energy loss spectra (EELS) which are particularly important in anisotropic materials. The approach incorporates relativistic effects in terms of the transition tensor within the dipole-selection rule. In particular, the method accounts for relativistic corrections to the magic angle in orientation resolved EELS experiments. The approach is validated by a study of the graphite CK edge, for which we present an accurate magic angle measurement consistent with the predicted value.
Active neutron multiplicity analysis and Monte Carlo calculations
NASA Astrophysics Data System (ADS)
Krick, M. S.; Ensslin, N.; Langner, D. G.; Miller, M. C.; Siebelist, R.; Stewart, J. E.; Ceo, R. N.; May, P. K.; Collins, L. L., Jr.
Active neutron multiplicity measurements of high-enrichment uranium metal and oxide samples have been made at Los Alamos and Y-12. The data from the measurements of standards at Los Alamos were analyzed to obtain values for neutron multiplication and source-sample coupling. These results are compared to equivalent results obtained from Monte Carlo calculations. An approximate relationship between coupling and multiplication is derived and used to correct doubles rates for multiplication and coupling. The utility of singles counting for uranium samples is also examined.
Model for analytical calculation of nuclear photoabsorption at intermediate energies
NASA Astrophysics Data System (ADS)
Hütt, M.-Th.; Milstein, A. I.; Schumacher, M.
1997-02-01
The universal curve {σ}/{A} of nuclear photoabsorption is investigated within a Fermi gas model of nuclear matter. An energy range from pion threshold up to 400 MeV is considered. The interactions between nucleon, pion, Δ-isobar and photon are considered in the non-relativistic approximation with corrections of the order {1}/{M} taken into account with respect to proton mass. Analytical expressions are obtained, in which the influence of nuclear correlations and two-nucleon contributions is studied explicitly. The contributions of real and virtual pions are found to be sufficient to obtain agreement with experimental data in this energy range. An extension of the model calculation to nucleon knock-out reactions is discussed.
Calculated photoelectron pitch angle and energy spectra. [in upper atmosphere
NASA Technical Reports Server (NTRS)
Mantas, G. P.; Bowhill, S. A.
1974-01-01
Calculations of the steady-state photoelectron energy and angular distribution in the altitude region between 120 and 1000 km are presented. The distribution is found to be isotropic at all altitudes below 250 km, while above this altitude anisotropies in both pitch angle and energy are found. The isotropy found in the angular distribution below 250 km implies that photoelectron transport below 250 km is insignificant, while the angular anisotropy found above this altitude implies a net photoelectron current in the upward direction. The energy anisotropy above 500 km arises from the selective backscattering of the low energy photoelectron population of the upward flux component by Coulomb collisions with the ambient ions. The total photoelectron flux attains its maximum value between about 40 and 70 km above the altitude at which the photoelectron production rate is maximum. The displacement of the maximum of the equilibrium flux is attributed to an increasing (with altitude) photoelectron lifetime. Photoelectrons at altitudes above that where the flux is maximum are on the average more energetic than those below that altitude.
Tables for simplifying calculations of activities produced by thermal neutrons
Senftle, F.E.; Champion, W.R.
1954-01-01
The method of calculation described is useful for the types of work of which examples are given. It is also useful in making rapid comparison of the activities that might be expected from several different elements. For instance, suppose it is desired to know which of the three elements, cobalt, nickel, or vanadium is, under similar conditions, activated to the greatest extent by thermal neutrons. If reference is made to a cross-section table only, the values may be misleading unless properly interpreted by a suitable comparison of half-lives and abundances. In this table all the variables have been combined and the desired information can be obtained directly from the values of A 3??, the activity produced per gram per second of irradiation, under the stated conditions. Hence, it is easily seen that, under similar circumstances of irradiation, vanadium is most easily activated even though the cross section of one of the cobalt isotopes is nearly five times that of vanadium and the cross section of one of the nickel isotopes is three times that of vanadium. ?? 1954 Societa?? Italiana di Fisica.
High energy activation data library (HEAD-2009)
Mashnik, Stepan G; Korovin, Yury A; Natalenko, Anatoly A; Konobeyev, Alexander Yu; Stankovskiy, A Yu
2010-01-01
A proton activation data library for 682 nuclides from 1 H to 210Po in the energy range from 150 MeV up to 1 GeV was developed. To calculate proton activation data, the MCNPX 2.6.0 and CASCADE/INPE codes were chosen. Different intranuclear cascade, preequilibrium, and equilibrium nuclear reaction models and their combinations were used. The optimum calculation models have been chosen on the basis of statistical correlations for calculated and experimental proton data taken from the EXFOR library of experimental nuclear data. All the data are written in ENDF-6 format. The library is called HEPAD-2008 (High-Energy Proton Activation Data). A revision of IEAF-2005 neutron activation data library has been performed. A set of nuclides for which the cross-section data can be (and were) updated using more modern and improved models is specified, and the corresponding calculations have been made in the present work. The new version of the library is called IEAF-2009. The HEPAD-2008 and IEAF-2009 are merged to the final HEAD-2009 library.
High energy activation data library (HEAD-2009)
NASA Astrophysics Data System (ADS)
Korovin, Yu. A.; Natalenko, A. A.; Stankovskiy, A. Yu.; Mashnik, S. G.; Konobeyev, A. Yu.
2010-12-01
A proton activation data library for 682 nuclides from 1H to 210Po in the energy range from 150 MeV up to 1 GeV was developed. To calculate proton activation data, the MCNPX 2.6.0 and CASCADE/INPE codes were chosen. Different intranuclear cascade, preequilibrium, and equilibrium nuclear reaction models and their combinations were used. The optimum calculation models have been chosen on the basis of statistical correlations for calculated and experimental proton data taken from the EXFOR library of experimental nuclear data. All the data are written in ENDF-6 format. The library is called HEPAD-2008 (High-Energy Proton Activation Data). A revision of IEAF-2005 neutron activation data library has been performed. A set of nuclides for which the cross-section data can be (and were) updated using more modern and improved models is specified, and the corresponding calculations have been made in the present work. The new version of the library is called IEAF-2009. The HEPAD-2008 and IEAF-2009 are merged to the final HEAD-2009 library.
Scattering Theory Calculations of Casimir Energies at High Curvature
NASA Astrophysics Data System (ADS)
Graham, Noah; Emig, Thorsten; Forrow, Aden; Jaffe, Robert; Kardar, Mehran; Maghrebi, Mohammad; Rahi, Jamal; Shpunt, Alex
2013-03-01
Scattering theory provides a powerful tool for capturing the response of an object to electromagnetic charge and field fluctuations. Techniques based on scattering theory have made possible a wide range of new calculations of Casimir energies. In this approach, the Casimir interaction energy for a collection of objects can be expressed in terms of the scattering T-matrices for each object individually, combined with universal translation matrices describing the objects' relative positions and orientations. These translation matrices are derived from an expansion of the free Green's function in an appropriate coordinate system, independent of the details of the objects themselves. This method proves particularly valuable for geometries involving high curvature, such as edges and tips. I will describe this approach in general terms and then give results from several problems to which it has been applied successfully. I will also discuss new developments in scattering theory that have been motivated by these problems. I would like to request that this abstract be part of a session on Casimir physics. Supported by the National Science Foundation, the US Department of Energy, the Defense Advanced Research Projects Agency, and the Deutsche Forschungsgemeinschaft
Similarity criteria in calculations of the energy characteristics of a cw oxygen - iodine laser
NASA Astrophysics Data System (ADS)
Mezhenin, A. V.; Azyazov, V. N.
2012-12-01
The calculated and experimental data on the energy efficiency of a cw oxygen - iodine laser (OIL) are analysed based on two similarity criteria, namely, on the ratio of the residence time of the gas mixture in the resonator to the characteristic time of extraction of the energy stored in singlet oxygen td and on the gain-to-loss ratio Π. It is shown that the simplified two-level laser model satisfactorily predicts the output characteristics of OILs with a stable resonator at τd <= 7. Efficient energy extraction from the OIL active medium is achieved in the case of τd = 5 - 7, Π = 4 - 8.
Similarity criteria in calculations of the energy characteristics of a cw oxygen - iodine laser
Mezhenin, A V; Azyazov, V N
2012-12-31
The calculated and experimental data on the energy efficiency of a cw oxygen - iodine laser (OIL) are analysed based on two similarity criteria, namely, on the ratio of the residence time of the gas mixture in the resonator to the characteristic time of extraction of the energy stored in singlet oxygen td and on the gain-to-loss ratio {Pi}. It is shown that the simplified two-level laser model satisfactorily predicts the output characteristics of OILs with a stable resonator at {tau}{sub d} {<=} 7. Efficient energy extraction from the OIL active medium is achieved in the case of {tau}{sub d} = 5 - 7, {Pi} = 4 - 8. (lasers)
Soil biological activity at European scale - two calculation concepts
NASA Astrophysics Data System (ADS)
Krüger, Janine; Rühlmann, Jörg
2014-05-01
The CATCH-C project aims to identify and improve the farm-compatibility of Soil Management Practices including to promote productivity, climate change mitigation and soil quality. The focus of this work concentrates on turnover conditions for soil organic matter (SOM). SOM is fundamental for the maintenance of quality and functions of soils while SOM storage is attributed a great importance in terms of climate change mitigation. The turnover conditions depend on soil biological activity characterized by climate and soil properties. To assess the turnover conditions two model concepts are applied: (I) Biological active time (BAT) regression approach derived from CANDY model (Franko & Oelschlägel 1995) expresses the variation of air temperature, precipitation and soil texture as a timescale and an indicator of biological activity for soil organic matter (SOM) turnover. (II) Re_clim parameter within the Introductory Carbon Balance Model (Andrén & Kätterer 1997) states the soil temperature and soil water to estimate soil biological activity. The modelling includes two strategies to cover the European scale and conditions. BAT was calculated on a 20x20 km grid basis. The European data sets of precipitation and air temperature (time period 1901-2000, monthly resolution), (Mitchell et al. 2004) were used to derive long-term averages. As we focus on agricultural areas we included CORINE data (2006) to extract arable land. The resulting BATs under co-consideration of the main soil textures (clay, silt, sand and loam) were investigated per environmental zone (ENZs, Metzger et al. 2005) that represents similar conditions for precipitation, temperature and relief to identify BAT ranges and hence turnover conditions for each ENZ. Re_clim was quantified by climatic time series of more than 250 weather stations across Europe presented by Klein Tank et al. (2002). Daily temperature, precipitation and potential evapotranspiration (maximal thermal extent) were used to calculate
Valence calculations of lanthanide anion binding energies: a comprehensive study
NASA Astrophysics Data System (ADS)
O'Malley, Steven M.; Beck, Donald R.
2009-05-01
We have applied a methodology of universal jls restrictions on the 4f^n subgroup of relativistic configuration-interaction calculations of progressively more complex lanthanide anionsootnotetextS. M. O'Malley and D. R. Beck, Phys. Rev. A 77, 012505 (2008).^,ootnotetextS. M. O'Malley and D. R. Beck, Phys. Rev. A 78, 012510 (2008).^,ootnotetextS. M. O'Malley and D. R. Beck, Phys. Rev. A, in press.. Our completed study of the row predicts bound 6p attachments to all lanthanide ground state configurations except Yb, additional 6p attachments to excited opposite parity configurations in Tb and Lu, and 6s attachments to excited open-6s thresholds in La, Ce, Pr, and Gd. In total we predict more than 100 bound states for the lanthanide anions, and we hope this comprehensive study encourages further experimentalootnotetexte.g. V. T. Davis et al., Nucl. Instrum. Methods Phys. Res. B 241, 118 (2005).^,ootnotetexte.g. C. W. Walter et al., Phys. Rev. A 76, 052702 (2007). interest in these anions. Such measurements will be useful in ``fine tuning'' these ab initio binding energies to account for missing core-valence correlation and the approximations that were necessary in these complex calculations.
Free-energy calculation via mean-force dynamics using a logarithmic energy landscape.
Morishita, Tetsuya; Itoh, Satoru G; Okumura, Hisashi; Mikami, Masuhiro
2012-06-01
A method for free-energy calculation based on mean-force dynamics (fictitious dynamics on a potential of mean force) is presented. The method utilizes a logarithmic form of free energy to enhance crossing barriers on a free-energy landscape, which results in efficient sampling of "rare" events. Invoking a conserved quantity in mean-force dynamics, free energy can be estimated on-the-fly without postprocessing. This means that an estimate of the free-energy profile can be locally made in contrast to the other methods based on mean-force dynamics such as metadynamics. The method is benchmarked against conventional methods and its high efficiency is demonstrated in the free-energy calculation for a glycine dipeptide molecule. PMID:23005238
Calculation of optical second-harmonic susceptibilities and optical activity for crystals
Levine, Z.H.
1994-12-31
A new generation of nearly first-principles calculations predicts both the linear and second-harmonic susceptibilities for a variety of insulating crystals, including GaAs, GaP, AlAs, AlP, Se, {alpha}-quartz, and c-urea. The results are typically in agreement with experimental measurements. The calculations have been extended to optical activity, with somewhat less success to date. The theory, based on a simple self-energy correction to the local density approximation, and results are reviewed herein.
Distributed activation energy model parameters of some Turkish coals
Gunes, M.; Gunes, S.K.
2008-07-01
A multi-reaction model based on distributed activation energy has been applied to some Turkish coals. The kinetic parameters of distributed activation energy model were calculated via computer program developed for this purpose. It was observed that the values of mean of activation energy distribution vary between 218 and 248 kJ/mol, and the values of standard deviation of activation energy distribution vary between 32 and 70 kJ/mol. The correlations between kinetic parameters of the distributed activation energy model and certain properties of coal have been investigated.
Jia, Xiangyu; Wang, Meiting; Shao, Yihan; König, Gerhard; Brooks, Bernard R; Zhang, John Z H; Mei, Ye
2016-02-01
In this work, the solvation free energies of 20 organic molecules from the 4th Statistical Assessment of the Modeling of Proteins and Ligands (SAMPL4) have been calculated. The sampling of phase space is carried out at a molecular mechanical level, and the associated free energy changes are estimated using the Bennett Acceptance Ratio (BAR). Then the quantum mechanical (QM) corrections are computed through the indirect Non-Boltzmann Bennett's acceptance ratio (NBB) or the thermodynamics perturbation (TP) method. We show that BAR+TP gives a minimum analytic variance for the calculated solvation free energy at the Gaussian limit and performs slightly better than NBB in practice. Furthermore, the expense of the QM calculations in TP is only half of that in NBB. We also show that defining the biasing potential as the difference of the solute-solvent interaction energy, instead of the total energy, can converge the calculated solvation free energies much faster but possibly to different values. Based on the experimental solvation free energies which have been published before, it is discovered in this study that BLYP yields better results than MP2 and some other later functionals such as B3LYP, M06-2X, and ωB97X-D. PMID:26731197
INDIVIDUALISED CALCULATION OF TISSUE IMPARTED ENERGY IN BREAST TOMOSYNTHESIS.
Geeraert, N; Klausz, R; Muller, S; Bloch, I; Bosmans, H
2016-06-01
The imparted energy to the glandular tissue in the breast (glandular imparted energy, GIE) is proposed for an improved assessment of the individual radiation-induced risk resulting from X-ray breast imaging. GIE is computed from an estimation of the quantity and localisation of glandular tissue in the breast. After a digital breast tomosynthesis (DBT) acquisition, the volumetric glandular content (volumetric breast density, VBD) is computed from the central X-ray projection. The glandular tissue distribution is determined by labelling the DBT voxels to ensure the conservation of the VBD. Finally, the GIE is calculated by Monte Carlo computation on the resulting tissue-labelled DBT volume. For verification, the method was applied to 10 breast-shaped digital phantoms made of different glandular spheres in an adipose background, and to a digital anthropomorphic phantom. Results were compared to direct GIE computations on the phantoms considered as 'ground-truth'. The major limitations in accuracy are those of DBT, in particular the limited z-resolution. However, for most phantoms, the results can be considered as acceptable. PMID:27127209
Calculating Capstone depleted uranium aerosol concentrations from beta activity measurements.
Szrom, Frances; Falo, Gerald A; Parkhurst, Mary Ann; Whicker, Jeffrey J; Alberth, David P
2009-03-01
Beta activity measurements were used as surrogate measurements of uranium mass in aerosol samples collected during the field testing phase of the Capstone Depleted Uranium (DU) Aerosol Study. These aerosol samples generated by the perforation of armored combat vehicles were used to characterize the DU source term for the subsequent Human Health Risk Assessment (HHRA) of Capstone aerosols. Establishing a calibration curve between beta activity measurements and uranium mass measurements is straightforward if the uranium isotopes are in equilibrium with their immediate short-lived, beta-emitting progeny. For DU samples collected during the Capstone study, it was determined that the equilibrium between the uranium isotopes and their immediate short-lived, beta-emitting progeny had been disrupted when penetrators had perforated target vehicles. Adjustments were made to account for the disrupted equilibrium and for wall losses in the aerosol samplers. Values for the equilibrium fraction ranged from 0.16 to 1, and the wall loss correction factors ranged from 1 to 1.92. This paper describes the process used and adjustments necessary to calculate uranium mass from proportional counting measurements. PMID:19204483
Identification of HIV Inhibitors Guided by Free Energy Perturbation Calculations
Acevedo, Orlando; Ambrose, Zandrea; Flaherty, Patrick T.; Aamer, Hadega; Jain, Prashi; Sambasivarao, Somisetti V.
2013-01-01
Free energy perturbation (FEP) theory coupled to molecular dynamics (MD) or Monte Carlo (MC) statistical mechanics offers a theoretically precise method for determining the free energy differences of related biological inhibitors. Traditionally requiring extensive computational resources and expertise, it is only recently that its impact is being felt in drug discovery. A review of computer-aided anti-HIV efforts employing FEP calculations is provided here that describes early and recent successes in the design of human immunodeficiency virus type 1 (HIV-1) protease and non-nucleoside reverse transcriptase inhibitors. In addition, our ongoing work developing and optimizing leads for small molecule inhibitors of cyclophilin A (CypA) is highlighted as an update on the current capabilities of the field. CypA has been shown to aid HIV-1 replication by catalyzing the cis/trans isomerization of a conserved Gly-Pro motif in the N-terminal domain of HIV-1 capsid (CA) protein. In the absence of a functional CypA, e.g., by the addition of an inhibitor such as cyclosporine A (CsA), HIV-1 has reduced infectivity. Our simulations of acylurea-based and 1-indanylketone-based CypA inhibitors have determined that their nanomolar and micromolar binding affinities, respectively, are tied to their ability to stabilize Arg55 and Asn102. A structurally novel 1-(2,6-dichlorobenzamido) indole core was proposed to maximize these interactions. FEP-guided optimization, experimental synthesis, and biological testing of lead compounds for toxicity and inhibition of wild-type HIV-1 and CA mutants have demonstrated a dose-dependent inhibition of HIV-1 infection in two cell lines. While the inhibition is modest compared to CsA, the results are encouraging. PMID:22316150
5 CFR 591.220 - How does OPM calculate energy utility cost indexes?
Code of Federal Regulations, 2011 CFR
2011-01-01
... 5 Administrative Personnel 1 2011-01-01 2011-01-01 false How does OPM calculate energy utility... Areas Cost-Of-Living Allowances § 591.220 How does OPM calculate energy utility cost indexes? (a) OPM calculates energy utility cost indexes based on the relative cost of maintaining a standard size dwelling...
5 CFR 591.220 - How does OPM calculate energy utility cost indexes?
Code of Federal Regulations, 2014 CFR
2014-01-01
... 5 Administrative Personnel 1 2014-01-01 2014-01-01 false How does OPM calculate energy utility... Areas Cost-Of-Living Allowances § 591.220 How does OPM calculate energy utility cost indexes? (a) OPM calculates energy utility cost indexes based on the relative cost of maintaining a standard size dwelling...
5 CFR 591.220 - How does OPM calculate energy utility cost indexes?
Code of Federal Regulations, 2013 CFR
2013-01-01
... 5 Administrative Personnel 1 2013-01-01 2013-01-01 false How does OPM calculate energy utility... Areas Cost-Of-Living Allowances § 591.220 How does OPM calculate energy utility cost indexes? (a) OPM calculates energy utility cost indexes based on the relative cost of maintaining a standard size dwelling...
5 CFR 591.220 - How does OPM calculate energy utility cost indexes?
Code of Federal Regulations, 2012 CFR
2012-01-01
... 5 Administrative Personnel 1 2012-01-01 2012-01-01 false How does OPM calculate energy utility... Areas Cost-Of-Living Allowances § 591.220 How does OPM calculate energy utility cost indexes? (a) OPM calculates energy utility cost indexes based on the relative cost of maintaining a standard size dwelling...
5 CFR 591.220 - How does OPM calculate energy utility cost indexes?
Code of Federal Regulations, 2010 CFR
2010-01-01
... 5 Administrative Personnel 1 2010-01-01 2010-01-01 false How does OPM calculate energy utility... Areas Cost-Of-Living Allowances § 591.220 How does OPM calculate energy utility cost indexes? (a) OPM calculates energy utility cost indexes based on the relative cost of maintaining a standard size dwelling...
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
Theoretical model for calculation of helicity in solar active regions
NASA Astrophysics Data System (ADS)
Chatterjee, P.
We (Choudhuri, Chatterjee and Nandy, 2005) calculate helicities of solar active regions based on the idea of Choudhuri (2003) that poloidal flux lines get wrapped around a toroidal flux tube rising through the convection zone, thereby giving rise to the helicity. Rough estimates based on this idea compare favourably with the observed magnitude of helicity. We use our solar dynamo model based on the Babcock--Leighton α-effect to study how helicity varies with latitude and time. At the time of solar maximum, our theoretical model gives negative helicity in the northern hemisphere and positive helicity in the south, in accordance with observed hemispheric trends. However, we find that, during a short interval at the beginning of a cycle, helicities tend to be opposite of the preferred hemispheric trends. Next we (Chatterjee, Choudhuri and Petrovay 2006) use the above idea along with the sunspot decay model of Petrovay and Moreno-Insertis, (1997) to estimate the distribution of helicity inside a flux tube as it keeps collecting more azimuthal flux during its rise through the convection zone and as turbulent diffusion keeps acting on it. By varying parameters over reasonable ranges in our simple 1-d model, we find that the azimuthal flux penetrates the flux tube to some extent instead of being confined to a narrow sheath outside.
Path-breaking schemes for nonequilibrium free energy calculations
NASA Astrophysics Data System (ADS)
Chelli, Riccardo; Gellini, Cristina; Pietraperzia, Giangaetano; Giovannelli, Edoardo; Cardini, Gianni
2013-06-01
We propose a path-breaking route to the enhancement of unidirectional nonequilibrium simulations for the calculation of free energy differences via Jarzynski's equality [C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997)], 10.1103/PhysRevLett.78.2690. One of the most important limitations of unidirectional nonequilibrium simulations is the amount of realizations necessary to reach suitable convergence of the work exponential average featuring the Jarzynski's relationship. In this respect, a significant improvement of the performances could be obtained by finding a way of stopping trajectories with negligible contribution to the work exponential average, before their normal end. This is achieved using path-breaking schemes which are essentially based on periodic checks of the work dissipated during the pulling trajectories. Such schemes can be based either on breaking trajectories whose dissipated work exceeds a given threshold or on breaking trajectories with a probability increasing with the dissipated work. In both cases, the computer time needed to carry out a series of nonequilibrium trajectories is reduced up to a factor ranging from 2 to more than 10, at least for the processes under consideration in the present study. The efficiency depends on several aspects, such as the type of process, the number of check-points along the pathway and the pulling rate as well. The method is illustrated through radically different processes, i.e., the helix-coil transition of deca-alanine and the pulling of the distance between two methane molecules in water solution.
Okeechobee County Energy Education Activities--Middle School Level.
ERIC Educational Resources Information Center
Allen, Rodney F., Ed.
Over 60 energy education activities related to mathematics, science, social studies, and English comprise this manual for middle school teachers. Included are issues for discussion, puzzles, science investigations, story writing exercises, and energy cost calculation problems. Among the topics covered in these lessons are energy consumption…
Caveat Emptor: Calculating All the Costs of Energy.
ERIC Educational Resources Information Center
Zinberg, Dorothy S.
This paper examines the energy problem. Specific topics discussed include the recent history of oil and gas consumption in the United States, conservation, coal, solar energy, and nuclear energy. While solutions to the energy problem differ, there is an urgent need for broad, public debate. Ultimately, the decisions made regarding energy will be…
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.
The importance of geospatial data to calculate the optimal distribution of renewable energies
NASA Astrophysics Data System (ADS)
Díaz, Paula; Masó, Joan
2013-04-01
Specially during last three years, the renewable energies are revolutionizing the international trade while they are geographically diversifying markets. Renewables are experiencing a rapid growth in power generation. According to REN21 (2012), during last six years, the total renewables capacity installed grew at record rates. In 2011, the EU raised its share of global new renewables capacity till 44%. The BRICS nations (Brazil, Russia, India and China) accounted for about 26% of the total global. Moreover, almost twenty countries in the Middle East, North Africa, and sub-Saharan Africa have currently active markets in renewables. The energy return ratios are commonly used to calculate the efficiency of the traditional energy sources. The Energy Return On Investment (EROI) compares the energy returned for a certain source and the energy used to get it (explore, find, develop, produce, extract, transform, harvest, grow, process, etc.). These energy return ratios have demonstrated a general decrease of efficiency of the fossil fuels and gas. When considering the limitations of the quantity of energy produced by some sources, the energy invested to obtain them and the difficulties of finding optimal locations for the establishment of renewables farms (e.g. due to an ever increasing scarce of appropriate land) the EROI becomes relevant in renewables. A spatialized EROI, which uses variables with spatial distribution, enables the optimal position in terms of both energy production and associated costs. It is important to note that the spatialized EROI can be mathematically formalized and calculated the same way for different locations in a reproducible way. This means that having established a concrete EROI methodology it is possible to generate a continuous map that will highlight the best productive zones for renewable energies in terms of maximum energy return at minimum cost. Relevant variables to calculate the real energy invested are the grid connections between
On the possibility of negative activation energies in bimolecular reactions
NASA Technical Reports Server (NTRS)
Jaffe, R. L.
1978-01-01
The temperature dependence of the rate constants for model reacting systems was studied to understand some recent experimental measurements which imply the existence of negative activation energies. A collision theory model and classical trajectory calculations are used to demonstrate that the reaction probability can vary inversely with collision energy for bimolecular reactions occurring on attractive potential energy surfaces. However, this is not a sufficient condition to ensure that the rate constant has a negative temperature dependence. On the basis of these calculations, it seems unlikely that a true bimolecular reaction between neutral molecules will have a negative activation energy.
Solar Energy Project, Activities: Biology.
ERIC Educational Resources Information Center
Tullock, Bruce, Ed.; And Others
This guide contains lesson plans and outlines of science activities which present concepts of solar energy in the context of biology experiments. Each unit presents an introduction; objectives; skills and knowledge needed; materials; methods; questions; recommendations for further work; and a teacher information sheet. The teacher information…
Activation energy measurements of cheese
Technology Transfer Automated Retrieval System (TEKTRAN)
Temperature sweeps of cheeses using small amplitude oscillatory shear tests produced values for activation energy of flow (Ea) between 30 and 44 deg C. Soft goat cheese and Queso Fresco, which are high-moisture cheeses and do not flow when heated, exhibited Ea values between 30 and 60 kJ/mol. The ...
Activation Energies of Plasmonic Catalysts.
Kim, Youngsoo; Dumett Torres, Daniel; Jain, Prashant K
2016-05-11
The activation energy of a catalytic reaction serves not only as a metric of the efficacy of a catalyst but also as a potential indicator of mechanistic differences between the catalytic and noncatalytic reaction. However, activation energies are quite underutilized in the field of photocatalysis. We characterize in detail the effect of visible light excitation on the activation enthalpy of an electron transfer reaction photocatalyzed by plasmonic Au nanoparticles. We find that in the presence of visible light photoexcitation, the activation enthalpy of the Au nanoparticle-catalyzed electron transfer reaction is significantly reduced. The reduction in the activation enthalpy depends on the excitation wavelength, the incident laser power, and the strength of a hole scavenger. On the basis of these results, we argue that the activation enthalpy reduction is directly related to the photoelectrochemical potential built-up on the Au nanoparticle under steady-state light excitation, analogous to electrochemical activation. Under optimum light excitation conditions, a potential as high as 240 mV is measured. The findings constitute more precise insights into the mechanistic role and energetic contribution of plasmonic excitation to chemical reactions catalyzed by transition metal nanoparticles. PMID:27064549
Meng, X. Y.; Qin, G. W.; Li, S.; Ren, Y. P.; Pei, W. L.; Zuo, L.; Wen, X. H.
2011-03-14
To improve photoelectrochemical (PEC) activity of hematite, the modification of energy band by doping 3d transition metal ions Cu and Ti into {alpha}-Fe{sub 2}O{sub 3} were studied via the first-principles calculations with density function theory (DFT)+U method. The results show that the band gap of hematite is {approx}2.1 eV and n-type dopant Ti improves the electric conductivity, confirmed by recent experiments. The p-type dopant Cu enhances the utilization ratio of solar energy, shifts both valance, and conduction band edges to a higher energy level, satisfying hydrogen production in the visible light driven PEC water splitting without voltage bias.
Trypsin-Ligand Binding Free Energy Calculation with AMOEBA
Shi, Yue; Jiao, Dian; Schnieders, Michael J.; Ren, Pengyu
2010-01-01
The binding free energies of several benzamidine-like inhibitors to trypsin were examined using a polarizable potential. All the computed binding free energies are in good agreement with the experimental data. From free energy decomposition, electrostatic interaction was found to be the driving force for the binding. Structural analysis shows that the ligands form hydrogen bonds with residues and water molecules nearby in a competitive fashion. The dependence of binding free energy on molecular dipole moment and polarizability was also studied. While the binding free energy is independent on the dipole moment, it shows a negative correlation with the polarizability. PMID:19965178
Extending the energy range of materials activation modelling
NASA Astrophysics Data System (ADS)
Forrest, R. A.
2004-08-01
Activation calculations are an essential contribution to understanding the interactions of fusion materials with neutrons. The existing state-of-the-art tools such as EASY-2003 enable calculations to be carried out with neutrons up to 20 MeV. Plans to expose fusion components to high neutron fluxes include the IFMIF materials testing facility. This accelerator-based device will produce neutrons with a high-energy tail up to about 55 MeV. In order to carry out activation calculations on materials exposed to such neutrons it is necessary to extend the energy range of the data libraries. An extension of the European Activation System (EASY) to a new version, EASY-2004, for testing has been completed. The existing reactions have been extended up to 60 MeV and new classes of reactions added using calculated cross sections. Results of preliminary calculations in an IFMIF relevant neutron field are given.
NASA Astrophysics Data System (ADS)
Sanchez-Parcerisa, D.; Cortés-Giraldo, M. A.; Dolney, D.; Kondrla, M.; Fager, M.; Carabe, A.
2016-02-01
In order to integrate radiobiological modelling with clinical treatment planning for proton radiotherapy, we extended our in-house treatment planning system FoCa with a 3D analytical algorithm to calculate linear energy transfer (LET) in voxelized patient geometries. Both active scanning and passive scattering delivery modalities are supported. The analytical calculation is much faster than the Monte-Carlo (MC) method and it can be implemented in the inverse treatment planning optimization suite, allowing us to create LET-based objectives in inverse planning. The LET was calculated by combining a 1D analytical approach including a novel correction for secondary protons with pencil-beam type LET-kernels. Then, these LET kernels were inserted into the proton-convolution-superposition algorithm in FoCa. The analytical LET distributions were benchmarked against MC simulations carried out in Geant4. A cohort of simple phantom and patient plans representing a wide variety of sites (prostate, lung, brain, head and neck) was selected. The calculation algorithm was able to reproduce the MC LET to within 6% (1 standard deviation) for low-LET areas (under 1.7 keV μm-1) and within 22% for the high-LET areas above that threshold. The dose and LET distributions can be further extended, using radiobiological models, to include radiobiological effectiveness (RBE) calculations in the treatment planning system. This implementation also allows for radiobiological optimization of treatments by including RBE-weighted dose constraints in the inverse treatment planning process.
NASA Astrophysics Data System (ADS)
Ghazzali, Mohamed; Khattab, Sherine A. N.; Elnakady, Yasser A.; Al-Mekhlafi, Fahd A.; Al-Farhan, Khalid; El-Faham, Ayman
2013-08-01
A series of naphthyl and tolyl sulfonate ester were synthesized and characterized by H NMR. X-ray single crystal diffraction experiments established the molecular structure of three new sulfonate esters derivatives, and spectral data agree with these in solution. The observed hydrogen bonding is discussed on the basis of crystal structural analyses and DFT/MP2 geometry optimization quantum calculations. Antimicrobial activities were screened for selected compounds against three human cancer cell lines and Mosquito Culex pipiens larvae.
Ovchinnikov, Victor; Cecchini, Marco; Karplus, Martin
2013-01-24
A simple and robust formulation of the path-independent confinement method for the calculation of free energies is presented. The simplified confinement method (SCM) does not require matrix diagonalization or switching off the molecular force field, and has a simple convergence criterion. The method can be readily implemented in molecular dynamics programs with minimal or no code modifications. Because the confinement method is a special case of thermodynamic integration, it is trivially parallel over the integration variable. The accuracy of the method is demonstrated using a model diatomic molecule, for which exact results can be computed analytically. The method is then applied to the alanine dipeptide in vacuum, and to the α-helix ↔ β-sheet transition in a 16-residue peptide modeled in implicit solvent. The SCM requires less effort for the calculation of free energy differences than previous formulations because it does not require computing normal modes. The SCM has a diminished advantage for determining absolute free energy values, because it requires decreasing the MD integration step to obtain accurate results. An approximate confinement procedure is introduced, which can be used to estimate directly the configurational entropy difference between two macrostates, without the need for additional computation of the difference in the free energy or enthalpy. The approximation has convergence properties similar to those of the standard confinement method for the calculation of free energies. The use of the approximation requires about 5 times less wall-clock simulation time than that needed to compute enthalpy differences to similar precision from an MD trajectory. For the biomolecular systems considered in this study, the errors in the entropy approximation are under 10%. Practical applications of the methods to proteins are currently limited to implicit solvent simulations. PMID:23268557
Subtleties in Energy Calculations in the Image Method
ERIC Educational Resources Information Center
Taddei, M. M.; Mendes, T. N. C.; Farina, C.
2009-01-01
In this pedagogical work, we point out a subtle mistake that can be made by undergraduate or graduate students in the computation of the electrostatic energy of a system containing charges and perfect conductors if they naively use the image method. Specifically, we show that naive expressions for the electrostatic energy for these systems…
A calculation of the diffusion energies for adatoms on surfaces of F.C.C. metals
NASA Technical Reports Server (NTRS)
Halicioglu, T.; Pound, G. M.
1979-01-01
The activation energies for diffusion were determined for gold, platinum and iridium adatoms on plane and plane PT surfaces and were found to be in good agreement with the measurements reported by Bassett and Webber. The Lennard-Jones pair potentials were used to model the interatomic forces, and relaxation of the substrate atoms in near proximity to the adatom was considered in detail. The present calculations clarify the mechanism of the observed two-dimensional diffusion of platinum and iridium atoms on a plane PT surface. The results are compared with those obtained using Morse potential functions and different relaxation techniques.
Oxyfluoroborate host glass for upconversion application: phonon energy calculation
NASA Astrophysics Data System (ADS)
Abdel-Baki, Manal; El-Diasty, Fouad
2016-04-01
Reducing the glass phonon energy is an essential procedure to achieve high efficient radiative upconversion process. The degree of covalence of chemical bonds is responsible for the high oscillator strength of intracenter transitions in rare-earth ions. So, conversion covalent to ionic glass character is proposed as a structure-sensitive criterion that controls the phonon energy of the glasses. A series of oxyfluoro aluminum-borate host glasses used for upconversion application is prepared by the conventional melt-quenching technique. Through lithium oxide substitution by lithium fluoride, the ionic-covalent property of Li+ ion successes to regulate the band gap energies of the studied glasses. Furthermore, a new method to determine the glass phonon energy is offered.
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
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.
Horn, Paul R; Head-Gordon, Martin
2016-02-28
In energy decomposition analysis (EDA) of intermolecular interactions calculated via density functional theory, the initial supersystem wavefunction defines the so-called "frozen energy" including contributions such as permanent electrostatics, steric repulsions, and dispersion. This work explores the consequences of the choices that must be made to define the frozen energy. The critical choice is whether the energy should be minimized subject to the constraint of fixed density. Numerical results for Ne2, (H2O)2, BH3-NH3, and ethane dissociation show that there can be a large energy lowering associated with constant density orbital relaxation. By far the most important contribution is constant density inter-fragment relaxation, corresponding to charge transfer (CT). This is unwanted in an EDA that attempts to separate CT effects, but it may be useful in other contexts such as force field development. An algorithm is presented for minimizing single determinant energies at constant density both with and without CT by employing a penalty function that approximately enforces the density constraint. PMID:26931692
NASA Astrophysics Data System (ADS)
Horn, Paul R.; Head-Gordon, Martin
2016-02-01
In energy decomposition analysis (EDA) of intermolecular interactions calculated via density functional theory, the initial supersystem wavefunction defines the so-called "frozen energy" including contributions such as permanent electrostatics, steric repulsions, and dispersion. This work explores the consequences of the choices that must be made to define the frozen energy. The critical choice is whether the energy should be minimized subject to the constraint of fixed density. Numerical results for Ne2, (H2O)2, BH3-NH3, and ethane dissociation show that there can be a large energy lowering associated with constant density orbital relaxation. By far the most important contribution is constant density inter-fragment relaxation, corresponding to charge transfer (CT). This is unwanted in an EDA that attempts to separate CT effects, but it may be useful in other contexts such as force field development. An algorithm is presented for minimizing single determinant energies at constant density both with and without CT by employing a penalty function that approximately enforces the density constraint.
Klopper, W.; Ruscic, B.; Tew, D. P.; Bischoff, F. A.; Wolfsegger, S.; Chemical Sciences and Engineering Division; Univ. Karlsruhe
2008-11-17
The atomization energies of the 105 molecules in the test set of Bakowies [D. Bakowies, J. Chem. Phys. 127 (2007) 084105] have been computed with an estimated standard deviation (from the values compiled in the Active Thermochemical Tables) of {+-}0.1 kJ/mol per valence electron in the molecule. Equilibrium geometries and harmonic vibrational frequencies were calculated at the all-electron CCSD(T)/cc-pCVTZ level, that is, at the level of coupled-cluster theory with singles, doubles and non-iterative triples in a correlation-consistent polarized core-valence triple-zeta basis. Single-point energy calculations were performed at the all-electron CCSD(T) level in a correlation-consistent polarized core-valence quadruple-zeta basis (cc-pCVQZ), and several corrections were added: (i) a correction for the basis-set truncation error, obtained from second-order perturbation theory using Slater-type geminals (MP2-F12 theory), (ii) a correction for the effect of anharmonicity on the zero-point vibrational energy, (iii) a relativistic correction, (iv) a correction for the difference between the full CCSDT model (coupled-cluster theory with singles, doubles and triples) and the CCSD(T) approximation, and (v) a correction for connected quadruple excitations obtained from CCSDT(Q) calculations. The correction for the basis-set truncation error was obtained from MP2-F12 calculations by scaling the MP2 basis-set truncation error by an empirically optimized 'interference factor' of f{sub int} = 0.78. The reference values from the Active Thermochemical Tables for 73 molecules in the test set, the equilibrium geometries, the harmonic vibrational frequencies, and all of the energy corrections represent valuable data for performance assessments of additivity schemes that will be developed in the future, in which the basis-set truncation error will be calculated at the level of coupled-cluster theory using Slater-type geminals (CC-F12 theory). Such a scheme will be free of empirical
Energy levels of isoelectronic impurities by large scale LDA calculations
Li, Jingbo; Wang, Lin-Wang
2002-11-22
Isoelectronic impurity states are localized states induced by stoichiometric single atom substitution in bulk semiconductor. Photoluminescence spectra indicate deep impurity levels of 0.5 to 0.9eV above the top of valence band for systems like: GaN:As, GaN:P, CdS:Te, ZnS:Te. Previous calculations based on small supercells seemingly confirmed these experimental results. However, the current ab initio calculations based on thousand atom supercells indicate that the impurity levels of the above systems are actually much shallower(0.04 to 0.23 eV), and these impurity levels should be compared with photoluminescence excitation spectra, not photoluminescence spectra.
NASA Astrophysics Data System (ADS)
Srivastava, Amit; Granek, Rony
2015-02-01
Motivated by single molecule experiments, we study thermal unfolding pathways of four proteins, chymotrypsin inhibitor, barnase, ubiquitin, and adenylate kinase, using bond network models that combine bond energies and elasticity. The protein elasticity is described by the Gaussian network model (GNM), to which we add prescribed bond binding energies that are assigned to all (nonbackbone) connecting bonds in the GNM of native state and assumed identical for simplicity. Using exact calculation of the Helmholtz free energy for this model, we consider bond rupture single events. The bond designated for rupture is chosen by minimizing the free-energy difference for the process, over all (nonbackbone) bonds in the network. Plotting the free-energy profile along this pathway at different temperatures, we observe a few major partial unfolding, metastable or stable, states, that are separated by free-energy barriers and change role as the temperature is raised. In particular, for adenylate kinase we find three major partial unfolding states, which is consistent with single molecule FRET experiments [Pirchi et al., Nat. Commun. 2, 493 (2011), 10.1038/ncomms1504] for which hidden Markov analysis reveals between three and five such states. Such states can play a major role in enzymatic activity.
A Variational Approach to Enhanced Sampling and Free Energy Calculations
NASA Astrophysics Data System (ADS)
Parrinello, Michele
2015-03-01
The presence of kinetic bottlenecks severely hampers the ability of widely used sampling methods like molecular dynamics or Monte Carlo to explore complex free energy landscapes. One of the most popular methods for addressing this problem is umbrella sampling which is based on the addition of an external bias which helps overcoming the kinetic barriers. The bias potential is usually taken to be a function of a restricted number of collective variables. However constructing the bias is not simple, especially when the number of collective variables increases. Here we introduce a functional of the bias which, when minimized, allows us to recover the free energy. We demonstrate the usefulness and the flexibility of this approach on a number of examples which include the determination of a six dimensional free energy surface. Besides the practical advantages, the existence of such a variational principle allows us to look at the enhanced sampling problem from a rather convenient vantage point.
Variational Approach to Enhanced Sampling and Free Energy Calculations
NASA Astrophysics Data System (ADS)
Valsson, Omar; Parrinello, Michele
2014-08-01
The ability of widely used sampling methods, such as molecular dynamics or Monte Carlo simulations, to explore complex free energy landscapes is severely hampered by the presence of kinetic bottlenecks. A large number of solutions have been proposed to alleviate this problem. Many are based on the introduction of a bias potential which is a function of a small number of collective variables. However constructing such a bias is not simple. Here we introduce a functional of the bias potential and an associated variational principle. The bias that minimizes the functional relates in a simple way to the free energy surface. This variational principle can be turned into a practical, efficient, and flexible sampling method. A number of numerical examples are presented which include the determination of a three-dimensional free energy surface. We argue that, beside being numerically advantageous, our variational approach provides a convenient and novel standpoint for looking at the sampling problem.
Advancing QCD-based calculations of energy loss
NASA Astrophysics Data System (ADS)
Tywoniuk, Konrad
2013-08-01
We give a brief overview of the basics and current developments of QCD-based calculations of radiative processes in medium. We put an emphasis on the underlying physics concepts and discuss the theoretical uncertainties inherently associated with the fundamental parameters to be extracted from data. An important area of development is the study of the single-gluon emission in medium. Moreover, establishing the correct physical picture of multi-gluon emissions is imperative for comparison with data. We will report on progress made in both directions and discuss perspectives for the future.
Turbulent energy exchange: Calculation and relevance for profile prediction
Candy, J.
2013-08-15
The anomalous heat production due to turbulence is neither routinely calculated in nonlinear gyrokinetic simulations nor routinely retained in profile prediction studies. In this work, we develop a symmetrized method to compute the exchange which dramatically reduces the intermittency in the time-dependent moment, thereby improving the accuracy of the time-average. We also examine the practical impact on transport-timescale simulations, and show that the exchange has only a minor impact on profile evolution for a well-studied DIII-D discharge.
Tripathi, Sunil Kumar; Muttineni, Ravikumar; Singh, Sanjeev Kumar
2013-10-01
Molecular docking, free energy calculation and molecular dynamics (MD) simulation studies have been performed, to explore the putative binding modes of 3,5-diaminoindazoles, imidazo(1,2-b)pyridazines and triazolo(1,5-a) pyridazines series of Cyclin-dependent kinase (CDK2) inhibitors. To evaluate the effectiveness of docking protocol in flexible docking, we have selected crystallographic bound compound to validate our docking procedure as evident from root mean square deviations (RMSDs). We found different binding sites namely catalytic, inhibitory phosphorylation, cyclin binding and CKS-binding site of the CDK2 contributing towards the binding of these compounds. Moreover, correlation between free energy of binding and biological activity yielded a statistically significant correlation coefficient. Finally, three representative protein-ligand complexes were subjected to molecular dynamics simulation to determine the stability of the predicted conformations. The low value of the RMSDs between the initial complex structure and the energy minimized final average complex structure suggests that the derived docked complexes are close to equilibrium. We suggest that the phenylacetyl type of substituents and cyclohexyl moiety make the favorable interactions with a number of residues in the active site, and show better inhibitory activity to improve the pharmacokinetic profile of compounds against CDK2. The structure-based drug design strategy described in this study will be highly useful for the development of new inhibitors with high potency and selectivity. PMID:23727278
Energy Activities for the Primary Classroom. Revised.
ERIC Educational Resources Information Center
Tierney, Blue, Comp.
An energy education program at the primary level should help students to understand the nature and importance of energy, consider different energy sources, learn about energy conservation, prepare for energy related careers, and become energy conscious in other career fields. The activities charts, readings, and experiments provided in this…
Improved thrust calculations of active magnetic bearings considering fringing flux
NASA Astrophysics Data System (ADS)
Jang, Seok-Myeong; Kim, Kwan-Ho; Ko, Kyoung-Jin; Choi, Ji-Hwan; Sung, So-Young; Lee, Yong-Bok
2012-04-01
A methodology for deriving fringing permeance in axisymmetric devices such as active thrust magnetic bearings (ATMBs) is presented. The methodology is used to develop an improved equivalent magnetic circuit (EMC) for ATMBs, which considers the fringing effect. This EMC was used to characterize the force between the housing and mover and the dependence of thrust and inductance on the air gap and input current, respectively. These characteristics were validated by comparison with those obtained by the finite element method and in experiments.
The Suppression of Energy Discretization Errors in Multigroup Transport Calculations
Larsen, Edward
2013-06-17
The Objective of this project is to develop, implement, and test new deterministric methods to solve, as efficiently as possible, multigroup neutron transport problems having an extremely large number of groups. Our approach was to (i) use the standard CMFD method to "coarsen" the space-angle grid, yielding a multigroup diffusion equation, and (ii) use a new multigrid-in-space-and-energy technique to efficiently solve the multigroup diffusion problem. The overall strategy of (i) how to coarsen the spatial an energy grids, and (ii) how to navigate through the various grids, has the goal of minimizing the overall computational effort. This approach yields not only the fine-grid solution, but also coarse-group flux-weighted cross sections that can be used for other related problems.
An artificial energy method for calculating flows with shocks
NASA Technical Reports Server (NTRS)
Rose, M. E.
1980-01-01
The artificial-viscosity method, first proposed by von Neumann and Richtmyer, introduces an artificial viscous pressure term in regions of compression such that an increase in entropy occurs in shock transition zones. The paper describes how dissipative flows can be induced by reducing the total energy available for adiabatic processes in shock zones. A class of inviscid fluid flows, called semiflows, is described in which the flows exhibit thermodynamic differences. Induced dissipative flows modify the pressure in regions of compression in a manner analogous to the artificial-viscosity method and for a gas, the effect is equivalent to suitably modifying the gas constant in the equation of state. By employing MacCormack's method and the usual non-adiabatic equations, numerical solutions of a Riemann problem are compared with the modified artificial energy method, showing that the dissipation effect predicted by the analytical formulation is reflected in the numerical method as well.
Ab initio molecular dynamics calculations of ion hydration free energies
Leung, Kevin; Rempe, Susan B.; Lilienfeld, O. Anatole von
2009-05-28
We apply ab initio molecular dynamics (AIMD) methods in conjunction with the thermodynamic integration or '{lambda}-path' technique to compute the intrinsic hydration free energies of Li{sup +}, Cl{sup -}, and Ag{sup +} ions. Using the Perdew-Burke-Ernzerhof functional, adapting methods developed for classical force field applications, and with consistent assumptions about surface potential ({phi}) contributions, we obtain absolute AIMD hydration free energies ({Delta}G{sub hyd}) within a few kcal/mol, or better than 4%, of Tissandier et al.'s [J. Phys. Chem. A 102, 7787 (1998)] experimental values augmented with the SPC/E water model {phi} predictions. The sums of Li{sup +}/Cl{sup -} and Ag{sup +}/Cl{sup -} AIMD {Delta}G{sub hyd}, which are not affected by surface potentials, are within 2.6% and 1.2 % of experimental values, respectively. We also report the free energy changes associated with the transition metal ion redox reaction Ag{sup +}+Ni{sup +}{yields}Ag+Ni{sup 2+} in water. The predictions for this reaction suggest that existing estimates of {Delta}G{sub hyd} for unstable radiolysis intermediates such as Ni{sup +} may need to be extensively revised.
Calculating Capstone Depleted Uranium Aerosol Concentrations from Beta Activity Measurements
Szrom, Fran; Falo, Gerald A.; Parkhurst, MaryAnn; Whicker, Jeffrey J.; Alberth, David P.
2009-03-01
Beta activity measurements were used as surrogate measurements of uranium mass in aerosol samples collected during the field testing phase of the Capstone Depleted Uranium (DU) Aerosol Study. These aerosol samples generated by the perforation of armored combat vehicles were used to characterize the depleted uranium (DU) source term for the subsequent human health risk assessment (HHRA) of Capstone aerosols. Establishing a calibration curve between beta activity measurements and uranium mass measurements is straightforward if the uranium isotopes are in equilibrium with their immediate short-lived, beta-emitting progeny. For DU samples collected during the Capstone study, it was determined that the equilibrium between the uranium isotopes and their immediate short lived, beta-emitting progeny had been disrupted when penetrators had perforated target vehicles. Adjustments were made to account for the disrupted equilibrium and for wall losses in the aerosol samplers. Correction factors for the disrupted equilibrium ranged from 0.16 to 1, and the wall loss correction factors ranged from 1 to 1.92.
Calculations of energy levels and lifetimes of low-lying states of barium and radium
Dzuba, V. A.; Ginges, J. S. M.
2006-03-15
We use the configuration-interaction method and many-body perturbation theory to perform accurate calculations of energy levels, transition amplitudes, and lifetimes of low-lying states of barium and radium. Calculations for radium are needed for the planning of measurements of parity- and time-invariance-violating effects which are strongly enhanced in this atom. Calculations for barium are used to control the accuracy of the calculations.
Free-energy calculation methods for collective phenomena in membranes
NASA Astrophysics Data System (ADS)
Smirnova, Yuliya G.; Fuhrmans, Marc; Barragan Vidal, Israel A.; Müller, Marcus
2015-09-01
Collective phenomena in membranes are those which involve the co-operative reorganization of many molecules. Examples of these are membrane fusion, pore formation, bending, adhesion or fission. The time and length scales, on which these processes occur, pose a challenge for atomistic simulations. Therefore, in order to solve the length scale problem it is popular to introduce a coarse-grained representation. To facilitate sampling of the relevant states additional computational techniques, which encourage the system to explore the free-energy landscape far from equilibrium and visit transition states, are needed. These computational techniques provide insights about the free-energy changes involved in collective transformations of membranes, yielding information about the rate limiting states, the transformation mechanism and the influence of architectural, compositional and interaction parameters. A common approach is to identify an order parameter (or reaction coordinate), which characterizes the pathway of membrane reorganization. However, no general strategy exists to define such an order parameter that can properly describe cooperative reorganizations in membranes. Recently developed methods can overcome this problem of the order-parameter choice and allow us to study collective phenomena in membranes. We will discuss such methods as thermodynamic integration, umbrella sampling, and the string method and results provided by their applications to particle-based simulations, particularly focusing on membrane fusion and pore formation.
Parallel implementation of electronic structure energy, gradient, and Hessian calculations.
Lotrich, V; Flocke, N; Ponton, M; Yau, A D; Perera, A; Deumens, E; Bartlett, R J
2008-05-21
ACES III is a newly written program in which the computationally demanding components of the computational chemistry code ACES II [J. F. Stanton et al., Int. J. Quantum Chem. 526, 879 (1992); [ACES II program system, University of Florida, 1994] have been redesigned and implemented in parallel. The high-level algorithms include Hartree-Fock (HF) self-consistent field (SCF), second-order many-body perturbation theory [MBPT(2)] energy, gradient, and Hessian, and coupled cluster singles, doubles, and perturbative triples [CCSD(T)] energy and gradient. For SCF, MBPT(2), and CCSD(T), both restricted HF and unrestricted HF reference wave functions are available. For MBPT(2) gradients and Hessians, a restricted open-shell HF reference is also supported. The methods are programed in a special language designed for the parallelization project. The language is called super instruction assembly language (SIAL). The design uses an extreme form of object-oriented programing. All compute intensive operations, such as tensor contractions and diagonalizations, all communication operations, and all input-output operations are handled by a parallel program written in C and FORTRAN 77. This parallel program, called the super instruction processor (SIP), interprets and executes the SIAL program. By separating the algorithmic complexity (in SIAL) from the complexities of execution on computer hardware (in SIP), a software system is created that allows for very effective optimization and tuning on different hardware architectures with quite manageable effort. PMID:18500853
Parallel implementation of electronic structure energy, gradient, and Hessian calculations
NASA Astrophysics Data System (ADS)
Lotrich, V.; Flocke, N.; Ponton, M.; Yau, A. D.; Perera, A.; Deumens, E.; Bartlett, R. J.
2008-05-01
ACES III is a newly written program in which the computationally demanding components of the computational chemistry code ACES II [J. F. Stanton et al., Int. J. Quantum Chem. 526, 879 (1992); [ACES II program system, University of Florida, 1994] have been redesigned and implemented in parallel. The high-level algorithms include Hartree-Fock (HF) self-consistent field (SCF), second-order many-body perturbation theory [MBPT(2)] energy, gradient, and Hessian, and coupled cluster singles, doubles, and perturbative triples [CCSD(T)] energy and gradient. For SCF, MBPT(2), and CCSD(T), both restricted HF and unrestricted HF reference wave functions are available. For MBPT(2) gradients and Hessians, a restricted open-shell HF reference is also supported. The methods are programed in a special language designed for the parallelization project. The language is called super instruction assembly language (SIAL). The design uses an extreme form of object-oriented programing. All compute intensive operations, such as tensor contractions and diagonalizations, all communication operations, and all input-output operations are handled by a parallel program written in C and FORTRAN 77. This parallel program, called the super instruction processor (SIP), interprets and executes the SIAL program. By separating the algorithmic complexity (in SIAL) from the complexities of execution on computer hardware (in SIP), a software system is created that allows for very effective optimization and tuning on different hardware architectures with quite manageable effort.
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
Energy Storage. Teachers Guide. Science Activities in Energy.
ERIC Educational Resources Information Center
Jacobs, Mary Lynn, Ed.
Included in this science activities energy package for students in grades 4-10 are 12 activities related to energy storage. Each activity is outlined on the front and back of a single sheet and is introduced by a key question. Most of the activities can be completed in the classroom with materials readily available in any community. Among the…
Sanchez-Parcerisa, D; Cortés-Giraldo, M A; Dolney, D; Kondrla, M; Fager, M; Carabe, A
2016-02-21
In order to integrate radiobiological modelling with clinical treatment planning for proton radiotherapy, we extended our in-house treatment planning system FoCa with a 3D analytical algorithm to calculate linear energy transfer (LET) in voxelized patient geometries. Both active scanning and passive scattering delivery modalities are supported. The analytical calculation is much faster than the Monte-Carlo (MC) method and it can be implemented in the inverse treatment planning optimization suite, allowing us to create LET-based objectives in inverse planning. The LET was calculated by combining a 1D analytical approach including a novel correction for secondary protons with pencil-beam type LET-kernels. Then, these LET kernels were inserted into the proton-convolution-superposition algorithm in FoCa. The analytical LET distributions were benchmarked against MC simulations carried out in Geant4. A cohort of simple phantom and patient plans representing a wide variety of sites (prostate, lung, brain, head and neck) was selected. The calculation algorithm was able to reproduce the MC LET to within 6% (1 standard deviation) for low-LET areas (under 1.7 keV μm(-1)) and within 22% for the high-LET areas above that threshold. The dose and LET distributions can be further extended, using radiobiological models, to include radiobiological effectiveness (RBE) calculations in the treatment planning system. This implementation also allows for radiobiological optimization of treatments by including RBE-weighted dose constraints in the inverse treatment planning process. PMID:26840945
Inconsistencies and ambiguities in calculating enzyme activity: The case of laccase.
Baltierra-Trejo, Eduardo; Márquez-Benavides, Liliana; Sánchez-Yáñez, Juan Manuel
2015-12-01
Laccase is a key enzyme in the degradation of lignin by fungi. Reports indicate that the activity of this enzyme ranges from 3.5 to 484,000 U L(-1). Our aim was to analyze how laccase activity is calculated in the literature, and to determine statistically whether variations in activity are due to biological properties or to inconsistencies in calculation. We found a general lack of consensus on the definition of enzyme activity, and enzymes are sometimes characterized in terms of reaction rate and specific activity. Moreover, enzyme activity is calculated using at least seven different equations. Therefore, it is critical to standardize the calculation of laccase activity in order to compare results directly. PMID:26459230
Heart rate calculation from ensemble brain wave using wavelet and Teager-Kaiser energy operator.
Srinivasan, Jayaraman; Adithya, V
2015-01-01
Electroencephalogram (EEG) signal artifacts are caused by various factors, such as, Electro-oculogram (EOG), Electromyogram (EMG), Electrocardiogram (ECG), movement artifact and line interference. The relatively high electrical energy cardiac activity causes EEG artifacts. In EEG signal processing the general approach is to remove the ECG signal. In this paper, we introduce an automated method to extract the ECG signal from EEG using wavelet and Teager-Kaiser energy operator for R-peak enhancement and detection. From the detected R-peaks the heart rate (HR) is calculated for clinical diagnosis. To check the efficiency of our method, we compare the HR calculated from ECG signal recorded in synchronous with EEG. The proposed method yields a mean error of 1.4% for the heart rate and 1.7% for mean R-R interval. The result illustrates that, proposed method can be used for ECG extraction from single channel EEG and used in clinical diagnosis like estimation for stress analysis, fatigue, and sleep stages classification studies as a multi-model system. In addition, this method eliminates the dependence of additional synchronous ECG in extraction of ECG from EEG signal process. PMID:26737640
Perfetti, Christopher M; Rearden, Bradley T
2014-01-01
This work introduces a new approach for calculating sensitivity coefficients for generalized neutronic responses to nuclear data uncertainties using continuous-energy Monte Carlo methods. The approach presented in this paper, known as the GEAR-MC method, allows for the calculation of generalized sensitivity coefficients for multiple responses in a single Monte Carlo calculation with no nuclear data perturbations or knowledge of nuclear covariance data. The theory behind the GEAR-MC method is presented here, and proof of principle is demonstrated by using the GEAR-MC method to calculate sensitivity coefficients for responses in several 3D, continuous-energy Monte Carlo applications.
THE MAGNETIC ENERGY-HELICITY DIAGRAM OF SOLAR ACTIVE REGIONS
Tziotziou, Kostas; Georgoulis, Manolis K.; Raouafi, Nour-Eddine
2012-11-01
Using a recently proposed nonlinear force-free method designed for single-vector magnetograms of solar active regions, we calculate the instantaneous free magnetic energy and relative magnetic helicity budgets in 162 vector magnetograms corresponding to 42 different active regions. We find a statistically robust, monotonic correlation between the free magnetic energy and the relative magnetic helicity in the studied regions. This correlation implies that magnetic helicity, in addition to free magnetic energy, may be an essential ingredient for major solar eruptions. Eruptive active regions appear well segregated from non-eruptive ones in both free energy and relative helicity with major (at least M-class) flares occurring in active regions with free energy and relative helicity exceeding 4 Multiplication-Sign 10{sup 31} erg and 2 Multiplication-Sign 10{sup 42} Mx{sup 2}, respectively. The helicity threshold agrees well with estimates of the helicity contents of typical coronal mass ejections.
Calculating kinetics parameters and reactivity changes with continuous-energy Monte Carlo
Kiedrowski, Brian C; Brown, Forrest B; Wilson, Paul
2009-01-01
The iterated fission probability interpretation of the adjoint flux forms the basis for a method to perform adjoint weighting of tally scores in continuous-energy Monte Carlo k-eigenvalue calculations. Applying this approach, adjoint-weighted tallies are developed for two applications: calculating point reactor kinetics parameters and estimating changes in reactivity from perturbations. Calculations are performed in the widely-used production code, MCNP, and the results of both applications are compared with discrete ordinates calculations, experimental measurements, and other Monte Carlo calculations.
A Python tool to set up relative free energy calculations in GROMACS.
Klimovich, Pavel V; Mobley, David L
2015-11-01
Free energy calculations based on molecular dynamics (MD) simulations have seen a tremendous growth in the last decade. However, it is still difficult and tedious to set them up in an automated manner, as the majority of the present-day MD simulation packages lack that functionality. Relative free energy calculations are a particular challenge for several reasons, including the problem of finding a common substructure and mapping the transformation to be applied. Here we present a tool, alchemical-setup.py, that automatically generates all the input files needed to perform relative solvation and binding free energy calculations with the MD package GROMACS. When combined with Lead Optimization Mapper (LOMAP; Liu et al. in J Comput Aided Mol Des 27(9):755-770, 2013), recently developed in our group, alchemical-setup.py allows fully automated setup of relative free energy calculations in GROMACS. Taking a graph of the planned calculations and a mapping, both computed by LOMAP, our tool generates the topology and coordinate files needed to perform relative free energy calculations for a given set of molecules, and provides a set of simulation input parameters. The tool was validated by performing relative hydration free energy calculations for a handful of molecules from the SAMPL4 challenge (Mobley et al. in J Comput Aided Mol Des 28(4):135-150, 2014). Good agreement with previously published results and the straightforward way in which free energy calculations can be conducted make alchemical-setup.py a promising tool for automated setup of relative solvation and binding free energy calculations. PMID:26487189
Free energy calculations on snake venom metalloproteinase BaP1.
Lingott, Torsten; Merfort, Irmgard; Steinbrecher, Thomas
2012-06-01
BaP1 is a snake venom metalloproteinase from the venom of Bothrops asper, showing high structural homology with the catalytic domain of human adamalysins and matrix metalloproteinases. It induces the release of cytokines, like interleukin-1 and tumor necrosis factor alpha. Recently, the high-resolution crystal structure of BaP1 with a bound inhibitor became available, representing an interesting model concerning inhibitor design for medicinally important metalloproteinases such as tumor necrosis factor alpha-converting enzyme and MMP13. We here use computational modeling to gain a better understanding about the binding properties of various ligands to BaP1, with a focus on computing ligand binding free energies. The obtained results should be of general significance for future research on medicinally important metalloproteinases. We have investigated the binding of the original inhibitor in detail and calculated its binding strength using MMP/GBSA free energy calculations. Additionally, the binding strengths of alternative ligands have been computed, and two of them are predicted and experimentally verified to strongly inhibit the enzyme. A suggestion for chemical modifications of BaP1 inhibitors could be made to guide future synthesis efforts. Furthermore, a contribution to the proteolytic reaction mechanism of metzincins is given. The pK value of the catalytically active glutamic acid residue 143 has been found to be significantly raised when compared with a free glutamate side chain. Calculations on other matrix metalloproteinases confirmed that this is not confined to BaP1, but seems to be a common feature of metzincins. PMID:22385614
Energy Activities for Junior High Science.
ERIC Educational Resources Information Center
Beaver, David; And Others
This document is a collection of six energy education activities for junior high school science. Its purpose is to help promote knowledge about energy, provide laboratory experiences, provoke inquiry, and relate energy to society through the science curriculum. The six activities are designed to take one to three class periods. Two of the…
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.
LCS/CINDER`90 accelerator tunnel activation calculations for the APT 1700-MeV accelerator tunnel
Court, J.D.; Snow, E.C.; Wilson, W.B.; Pitcher, E.R.
1998-09-01
Calculations have been done to determine the amount of activation in the linac components and tunnel air for the Accelerator Production of Tritium 1700-MeV superconducting linac. Proton transport is accomplished through the use of the LAHET Code System. Particle production and depletion from proton and high-energy neutron reactions, calculated in LAHET, as well as low-energy neutron fluxes calculated by MCNP, are passed to the radionuclide production code CINDER`90 to determine the source terms at various times after irradiation. The upper limit on total air activation based on conservative assumptions, for the entire tunnel air volume, was found to be 4.77 Ci after a nine-month irradiation. This is reduced to 0.09 Ci after a 10-hour cooling off period. The total activation for the full 1-km of beamline components was found to be less than 4 kCi, with the half-lives of the highest contributors ranging from 12 years to 2 minutes. This beamline component activation calculation was done for an irradiation time of 40 years, which is the anticipated lifetime of the superconducting linac.
NASA Astrophysics Data System (ADS)
Shi, De-Heng; Liu, Yu-Fang; Sun, Jin-Feng; Zhu, Zun-Lue; Yang, Xiang-Dong
2006-12-01
The reasonable dissociation limit of the second excited singlet state B1Π of 7LiH molecule is obtained. The accurate dissociation energy and equilibrium geometry of the B1Π state are calculated using a symmetry-adapted-cluster configuration-interaction method in full active space. The whole potential energy curve for the B1Π state is obtained over the internuclear distance ranging from about 0.10 nm to 0.54 nm, and has a least-square fit to the analytic Murrell-Sorbie function form. The vertical excitation energy is calculated from the ground state to the B1Π state and compared with previous theoretical results. The equilibrium internuclear distance obtained by geometry optimization is found to be quite different from that obtained by single-point energy scanning under the same calculation condition. Based on the analytic potential energy function, the harmonic frequency value of the B1Π state is estimated. A comparison of the theoretical calculations of dissociation energies, equilibrium interatomic distances and the analytic potential energy function with those obtained by previous theoretical results clearly shows that the present work is more comprehensive and in better agreement with experiments than previous theories, thus it is an improvement on previous theories.
A Simple Approach for the Calculation of Energy Levels of Light Atoms
ERIC Educational Resources Information Center
Woodyard, Jack R., Sr.
1972-01-01
Describes a method for direct calculation of energy levels by using elementary techniques. Describes the limitations of the approach but also claims that with a minimum amount of labor a student can get greater understanding of atomic physics problems. (PS)
Ghosh, Soumya; Hammes-Schiffer, Sharon
2015-01-01
Electrochemical electron transfer reactions play an important role in energy conversion processes with many technological applications. Electrodes modified by self-assembled monolayers (SAMs) exhibit reduced double layer effects and are used in molecular electronics. An important quantity for calculating the electron transfer rate constant is the reorganization energy, which is associated with changes in the solute geometry and the environment. In this Letter, an approach for calculating the electrochemical reorganization energy for a redox molecule attached to or near a SAM modified electrode is presented. This integral equations formalism polarizable continuum model (IEF-PCM) approach accounts for the detailed electronic structure of the molecule, as well as the contributions from the electrode, SAM, and electronic and inertial solvent responses. The calculated total reorganization energies are in good agreement with experimental data for a series of metal complexes in aqueous solution. This approach will be useful for calculating electron transfer rate constants for molecular electrocatalysts. PMID:26263083
Determining characteristics of melting cheese by activation energy
Technology Transfer Automated Retrieval System (TEKTRAN)
Activation energy of flow (Ea) between 30 and 44 deg C was measured from temperature sweeps of various cheeses to determine its usefulness in predicting rheological behavior upon heating. Seven cheese varieties were heated in a rheometer from 22 to 70 deg C, and Ea was calculated from the resulting ...
Activation energy measurements in rheological analysis of cheese
Technology Transfer Automated Retrieval System (TEKTRAN)
Activation energy of flow (Ea) was calculated from temperature sweeps of cheeses with contrasting characteristics to determine its usefulness in predicting rheological behavior upon heating. Cheddar, Colby, whole milk Mozzarella, low moisture part skim Mozzarella, Parmesan, soft goat, and Queso Fre...
Energy Conservation Activity Packet, Grade 5.
ERIC Educational Resources Information Center
Bakke, Ruth
This activity packet for grade 5 is one of a series developed in response to concern for energy conservation. It contains activities that stress an energy conservation ethic and includes many values clarification activities for grade five. The packet is divided into two parts and provides the teacher with background information, concepts and…
NASA Astrophysics Data System (ADS)
Wu, Zhe-Ying; Qi, Chong; Wyss, Ramon; Liu, Hong-Liang
2015-08-01
Background: The deviation between different model calculations that may occur when one goes toward regions where the masses are unknown is getting increased attention. This is related to the uncertainties of the different models which may have not been fully understood. Purpose: To explore in detail the effect of the isospin dependence of the spin-orbital force in the Woods-Saxon potential on global binding energy and deformation calculations. Method: The microscopic energies and nuclear deformations of about 1850 even-even nuclei are calculated systematically within the macroscopic-microscopic framework using three Woods-Saxon parametrizations, with different isospin dependencies, which were constructed mainly for nuclear spectroscopy calculations. Calculations are performed in the deformation space (β2,γ ,β4) . Both the monopole and doubly stretched quadrupole interactions are considered for the pairing channel. Results: The ground-state deformations obtained by the three calculations are quite similar to each other. Large differences are seen mainly in neutron-rich nuclei and in superheavy nuclei. Systematic calculations on the shape-coexisting second minima are also presented. As for the microscopic energies of the ground states, the results are also very close to each other. Only in a few cases the difference is larger than 2 MeV. The total binding energy is estimated by adding the macroscopic energy provided by the usual liquid drop model with its parameters fitted through the least square root and minimax criteria. Calculations are also compared with the results of other macroscopic-microscopic mass models. Conclusions: All the three calculations give similar values for the deformations, microscopic energies, and binding energies of most nuclei. One may expect to have a better understanding of the isospin dependence of the spin-orbital force with more data on proton- and neutron-rich nuclei.
Wind Energy Finance (WEF): An Online Calculator for Economic Analysis of Wind Projects
Not Available
2004-02-01
This brochure provides an overview of Wind Energy Finance (WEF), a free online cost of energy calculator developed by the National Renewable Energy Laboratory that provides quick, detailed economic evaluation of potential utility-scale wind energy projects. The brochure lists the features of the tool, the inputs and outputs that a user can expect, visuals of the screens and a Cash Flow Results table, and contact information.
NASA Astrophysics Data System (ADS)
Ratsch, Christian; Kaminski, Jakub
In this talk we will present a new approach for the calculation of surface energies of periodic crystal. For non-polar materials slabs (which are terminated by two identical surfaces) the task of calculating the surface energy is trivial. But it is more problematic for polar systems where both terminating surfaces are different, as there is no single established method allowing for equal treatment of a wide range of surface morphologies and orientations. Our proposed new approach addresses this problem. It relies on carefully chosen capping atoms and the assumptions that their bond energy contributions can be used to approximate the total energy of the surface. The choice of the capping atoms is governed by a set of simple guidelines that are applicable for surfaces with different terminations. We present the results for different semiconductor materials and show that our approach leads to surfaces energies with errors as low as 2%. We show that hydrogen is not always the best choice for a capping atom if accurate surface energies are the target of the calculations. Our approach is suitable for high-throughput screening of new material interfaces, as accurate calculations of surface energies can be performed in an unsupervised algorithm. A New Approach for Surface Energy Calculations Applicable to High-throughput Design of New Interfaces.
Variational calculation of highly excited rovibrational energy levels of H2O2.
Polyansky, Oleg L; Kozin, Igor N; Ovsyannikov, Roman I; Małyszek, Paweł; Koput, Jacek; Tennyson, Jonathan; Yurchenko, Sergei N
2013-08-15
Results are presented for highly accurate ab initio variational calculation of the rotation-vibration energy levels of H2O2 in its electronic ground state. These results use a recently computed potential energy surface and the variational nuclear-motion programs WARV4, which uses an exact kinetic energy operator, and TROVE, which uses a numerical expansion for the kinetic energy. The TROVE calculations are performed for levels with high values of rotational excitation, J up to 35. The purely ab initio calculations of the rovibrational energy levels reproduce the observed levels with a standard deviation of about 1 cm(-1), similar to that of the J = 0 calculation, because the discrepancy between theory and experiment for rotational energies within a given vibrational state is substantially determined by the error in the vibrational band origin. Minor adjustments are made to the ab initio equilibrium geometry and to the height of the torsional barrier. Using these and correcting the band origins using the error in J = 0 states lowers the standard deviation of the observed-calculated energies to only 0.002 cm(-1) for levels up to J = 10 and 0.02 cm(-1) for all experimentally known energy levels, which extend up to J = 35. PMID:23611762
T. Downar
2009-03-31
The overall objective of the work here has been to eliminate the approximations used in current resonance treatments by developing continuous energy multi-dimensional transport calculations for problem dependent self-shielding calculations. The work here builds on the existing resonance treatment capabilities in the ORNL SCALE code system.
Energy and power limits for microbial activity
NASA Astrophysics Data System (ADS)
LaRowe, D.; Amend, J.
2014-12-01
The goal of this presentation is to describe a quantitative framework for determining how energy limits microbial activity, biomass and, ultimately, biogeochemical processes. Although this model can be applied to any environment, its utility is demonstrated in marine sediments, which are an attractive test habitat because they encompass a broad spectrum of energy levels, varying amounts of biomass and are ubiquitous. The potential number of active microbial cells in Arkonas Basin (Baltic Sea) sediments are estimated as a function of depth by quantifying the amount of energy that is available to them and the rate at which it is supplied: power. The amount of power supplied per cubic centimeter of sediment is determined by calculating the Gibbs energy of fermentation and sulfate reduction in combination with the rate of particulate organic carbon, POC, degradation. The Reactive Continuum Model (Boudreau and Ruddick, 1991), RCM, is used to determine the rate at which POC is made available for microbial consumption. The RCM represents POC as containing a range of different types of organic compounds whose ability to be consumed by microorganisms varies as a function of the age of the sediment and on the distribution of compound types that were initially deposited. The sediment age model and RCM parameters determined by (Mogollon et al., 2012) are used. The power available for fermentation and sulfate reduction coupled to H2 and acetate oxidation varies from 10-8 W cm-3 at the sediment water interface to between 10-11 - 10-12 W cm-3 at 3.5 meters below the seafloor, mbsf. Using values of maintenance powers for each of these catabolic activities taken from the literature, the total number of active cells in these sediments similarly decreases from just less than 108 cell cm-3 at the SWI to 4.6 x 104 cells cm-3 at 3.5 mbsf. The number of moles of POC decreases from 2.6 x 10-5 to 9.5 x 10-6, also becoming more recalcitrant with depth. Boudreau, B. P. and Ruddick, B. R
METHODOLOGICAL NOTES: Energy density calculations for ball-lightning-like luminous silicon balls
NASA Astrophysics Data System (ADS)
Paiva, Gerson S.; Ferreira, Joacy V.; Bastos, Cristiano C.; dos Santos, Marcus V.; Pavão, Antonio C.
2010-05-01
The energy density of a luminous silicon ball [Phys. Rev. Lett. 98 048501 (2007)] is calculated for a model with a metal core surrounded by an atmosphere of silicon oxides. Experimental data combined with the molecular orbital calculations of the oxidation enthalpy lead to a mean energy density of 3.9 MJ m-3, which is within the range of estimates from other ball lightning models. This result provides good evidence to support the silicon-based model.
George, K.; Schweizer, T.
2008-01-01
This report details the methodology used by DOE to calculate levelized cost of wind energy and demonstrates the variation in COE estimates due to different financing assumptions independent of wind generation technology.
The Effect of Calculator-Based Ranger Activities on Students' Graphing Ability.
ERIC Educational Resources Information Center
Kwon, Oh Nam
2002-01-01
Addresses three issues of Calculator-based Ranger (CBR) activities on graphing abilities: (a) the effect of CBR activities on graphing abilities; (b) the extent to which prior knowledge about graphing skills affects graphing ability; and (c) the influence of instructional styles on students' graphing abilities. Indicates that CBR activities are…
Nuclear data processing for energy release and deposition calculations in the MC21 Monte Carlo code
Trumbull, T. H.
2013-07-01
With the recent emphasis in performing multiphysics calculations using Monte Carlo transport codes such as MC21, the need for accurate estimates of the energy deposition-and the subsequent heating - has increased. However, the availability and quality of data necessary to enable accurate neutron and photon energy deposition calculations can be an issue. A comprehensive method for handling the nuclear data required for energy deposition calculations in MC21 has been developed using the NDEX nuclear data processing system and leveraging the capabilities of NJOY. The method provides a collection of data to the MC21 Monte Carlo code supporting the computation of a wide variety of energy release and deposition tallies while also allowing calculations with different levels of fidelity to be performed. Detailed discussions on the usage of the various components of the energy release data are provided to demonstrate novel methods in borrowing photon production data, correcting for negative energy release quantities, and adjusting Q values when necessary to preserve energy balance. Since energy deposition within a reactor is a result of both neutron and photon interactions with materials, a discussion on the photon energy deposition data processing is also provided. (authors)
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.
Douillard, Jean-Marc; Salles, Fabrice; Henry, Marc; Malandrini, Harold; Clauss, Frédéric
2007-01-15
The surface energies of talc and chlorite is computed using a simple model, which uses the calculation of the electrostatic energy of the crystal. It is necessary to calculate the atomic charges. We have chosen to follow Henry's model of determination of partial charges using scales of electronegativity and hardness. The results are in correct agreement with a determination of the surface energy obtained from an analysis of the heat of immersion data. Both results indicate that the surface energy of talc is lower than the surface energy of chlorite, in agreement with observed behavior of wettability. The influence of Al and Fe on this phenomenon is discussed. Surface energy of this type of solids seems to depend more strongly on the geometry of the crystal than on the type of atoms pointing out of the surface; i.e., the surface energy depends more on the physics of the system than on its chemistry. PMID:17081554
PERTURB: A program for calculating vibrational energies by generalized algebraic quantization
NASA Astrophysics Data System (ADS)
Fried, Laurence E.; Ezra, Gregory S.
1988-09-01
We describe PERTURB, a special purpose algebraic manipulation program which calculates vibrational eigenvalues in coupled oscillator systems. PERTURB implements the method of generalized algebraic quantization (AQ), in which Van Vleck perturbation theory is formulated in a mock phase space. The phase space formulation enables quantum and classical perturbation theory to be treated on the same footing, and allows the systematic calculation of corrections to classical perturbation results in powers of h̷. Generalized AQ is a powerful and efficient technique for calculating semiclassical vibrational energy levels. In many cases, including just the first correction to classical perturbation theory yields highly accurate energies.
Continuous-energy eigenvalue sensitivity coefficient calculations in TSUNAMI-3D
Perfetti, C. M.; Rearden, B. T.
2013-07-01
Two methods for calculating eigenvalue sensitivity coefficients in continuous-energy Monte Carlo applications were implemented in the KENO code within the SCALE code package. The methods were used to calculate sensitivity coefficients for several test problems and produced sensitivity coefficients that agreed well with both reference sensitivities and multigroup TSUNAMI-3D sensitivity coefficients. The newly developed CLUTCH method was observed to produce sensitivity coefficients with high figures of merit and a low memory footprint, and both continuous-energy sensitivity methods met or exceeded the accuracy of the multigroup TSUNAMI-3D calculations. (authors)
Development of a SCALE Tool for Continuous-Energy Eigenvalue Sensitivity Coefficient Calculations
Perfetti, Christopher M; Rearden, Bradley T
2013-01-01
Two methods for calculating eigenvalue sensitivity coefficients in continuous-energy Monte Carlo applications were implemented in the KENO code within the SCALE code package. The methods were used to calculate sensitivity coefficients for several criticality safety problems and produced sensitivity coefficients that agreed well with both reference sensitivities and multigroup TSUNAMI-3D sensitivity coefficients. The newly developed CLUTCH method was observed to produce sensitivity coefficients with high figures of merit and low memory requirements, and both continuous-energy sensitivity methods met or exceeded the accuracy of the multigroup TSUNAMI-3D calculations.
Development of a SCALE Tool for Continuous-Energy Eigenvalue Sensitivity Coefficient Calculations
NASA Astrophysics Data System (ADS)
Perfetti, Christopher M.; Rearden, Bradley T.
2014-06-01
Two methods for calculating eigenvalue sensitivity coefficients in continuous-energy Monte Carlo applications were implemented in the KENO code within the SCALE code package. The methods were used to calculate sensitivity coefficients for several criticality safety problems and produced sensitivity coefficients that agreed well with both reference sensitivities and multigroup TSUNAMI-3D sensitivity coefficients. The newly developed CLUTCH method was observed to produce sensitivity coefficients with high figures of merit and low memory requirements, and both continuous-energy sensitivity methods met or exceeded the accuracy of the multigroup TSUNAMI-3D calculations.
Rucker, Robert; Oelschlaeger, Peter; Warshel, Arieh
2010-01-01
DNA polymerase β (pol β) is a small eukaryotic enzyme with the ability to repair short single-stranded DNA gaps that has found use as a model system for larger replicative DNA polymerases. For all DNA polymerases, the factors determining their catalytic power and fidelity are the interactions between the bases of the base pair, amino acids near the active site, and the two magnesium ions. In this report, we study effects of all three aspects on human pol β transition state (TS) binding free energies by reproducing a consistent set of experimentally determined data for different structures. Our calculations comprise the combination of four different base pairs (incoming pyrimidine nucleotides incorporated opposite both matched and mismatched purines) with four different pol β structures (wild type and three separate mutations of ionized residues to alanine). We decompose the incoming deoxynucleoside 5′-triphosphate-TS, and run separate calculations for the neutral base part and the highly charged triphosphate part, using different dielectric constants in order to account for the specific electrostatic environments. This new approach improves our ability to predict the effect of matched and mismatched base pairing and of mutations in DNA polymerases on fidelity and may be a useful tool in studying the potential of DNA polymerase mutations in the development of cancer. It also supports our point of view with regards to the origin of the structural control of fidelity, allowing for a quantified description of the fidelity of DNA polymerases. PMID:19842163
Intrinsic free energy in active nematics
NASA Astrophysics Data System (ADS)
Thampi, Sumesh P.; Doostmohammadi, Amin; Golestanian, Ramin; Yeomans, Julia M.
2015-10-01
Basing our arguments on the theory of active liquid crystals, we demonstrate, both analytically and numerically, that the activity can induce an effective free energy which enhances ordering in extensile systems of active rods and in contractile suspensions of active discs. We argue that this occurs because any ordering fluctuation is enhanced by the flow field it produces. A phase diagram in the temperature-activity plane compares ordering due to a thermodynamic free energy to that resulting from the activity. We also demonstrate that activity can drive variations in concentration, but for a different physical reason that relies on the separation of hydrodynamic and diffusive time scales.
Simple energy-calculation method for solar industrial-process-heat steam systems
NASA Astrophysics Data System (ADS)
Gee, R. C.
1983-01-01
Designing a solar industrial process heat (IPH) system, sizing its components and predicting its annual energy delivery requires a method for calculating solar system performance. A calculation method that is accurate, easy to use, accounts for the impact of all important system parameters, and does not require use of a computer is described. Only simple graphs and a hand calculator are required to predict annual collector field performance and annual system losses. The energy calculation method is applicable to a variety of solar system configurations. The calculation method applied only to parabolic trough steam generation systems that do not employ thermal storage is described. Both flash tank and unfired boiler steam systems are covered.
Simple energy-calculation method for solar industrial-process-heat steam systems
Gee, R.
1983-01-01
Designing a solar industrial-process heat (IPH) system, sizing its components and predicting its annual energy delivery requires a method for calculating solar system performance. A calculation method that is accurate, easy to use, accounts for the impact of all important system parameters, and does not require use of a computer is described. Only simple graphs and a hand calculator are required to predict annual collector field performance and annual system losses. The energy-calculation method is applicable to a variety of solar-system configurations. The calculation method applied only to parabolic-trough steam-generation systems that do not employ thermal storage is described. Both flash tank and unfired-boiler steam systems are covered.
NASA Astrophysics Data System (ADS)
Montgomery, Jason
2007-12-01
Scattering resonances play a key role in many chemical processes, including unimolecular and bimolecular reactions and photodissociation. A significant theoretical emphasis over the past several decades has been placed on accurate resonance calculations for polyatomic systems. In spite of such efforts, a quantum treatment of molecular systems which exhibit a high density of states and strong coordinate coupling near dissociation remains a formidable task. The research described herein employs improved quantum mechanical methods to calculate a representation of nuclear motion, both bound and unbound, which is used subsequently to calculate accurate resonance energies and lifetimes for two triatomic systems: the neon trimer and ozone. Specifically, theory and results are given regarding the construction of an optimal, L2 eigenbasis using techniques such as the discrete variable representation, the energy selected basis (ESB) method, and iterative diagonalization methods. A new energy selection method is also developed and implemented for the neon trimer. Subsequent resonance calculations are described which make use of the artificial boundary inhomogeneity (ABI) method, adapted to work with the above mentioned ESB and hyperspherical coordinates. The ABI method is used to calculate a set of linearly independent wavefunctions (LIWs) at a given energy for the representation of the scattering wavefunction. Resonance parameters are obtained by imposing scattering boundary conditions on a linear combination of LIWs and solving for the S-matrix, S, its energy derivative, dS/dE, and the Smith lifetime matrix, Q. When available, comparisons are made with previously reported calculations.
Meson self-energies calculated by the relativistic particle-hole-antiparticle representation
Nakano, M.; Noda, N.; Mitsumori, T.; Koide, K.; Kouno, H.; Hasegawa, A.; Liu, L.
1997-12-01
A new formulation of meson self-energies is introduced for {sigma},{omega},{pi},{rho},{delta}, and {eta} mesons on the basis of the particle-hole-antiparticle representation. We have studied the difference between the meson self-energy (MSE) of this representation and the MSE of the traditional density-Feynman (DF) representation. It is shown that the new formulation describes exactly the physical processes such as particle-hole excitations or particle-antiparticle excitations, and that, on the other hand, the meson self-energy based on the DF representation includes unphysical components. By numerical calculations, the meson self-energies describing the particle-hole excitations are shown to be close to each other for most of the meson self-energy in low momentum (R{lt}500 MeV) and low energy (R{sub 0}{lt}200 MeV). This fact implies that former calculations using the low momentum and low-energy part do not change greatly. The density part of the density-Feynman representation has been shown to have a resonant structure around the energy of particle-antiparticle excitation, which causes a large difference between the two representations in the meson spectrum calculations. Our investigation concludes that the former calculations based on the density-Feynman representation are not invalidated in many cases, but the particle-hole-antiparticle representation is more appropriate to treat exactly the physical processes. {copyright} {ital 1997} {ital The American Physical Society}
ERIC Educational Resources Information Center
Barbiric, Dora; Tribe, Lorena; Soriano, Rosario
2015-01-01
In this laboratory, students calculated the nutritional value of common foods to assess the energy content needed to answer an everyday life application; for example, how many kilometers can an average person run with the energy provided by 100 g (3.5 oz) of beef? The optimized geometries and the formation enthalpies of the nutritional components…
Electron affinities (EAs) and free energies for electron attachment have been calculated for 42 polynuclear aromatic hydrocarbons and related molecules by a variety of theoretical models, including Koopmans' theorem methods and the L1E method from differences in energy between th...
Yang, W.; Wu, H.; Cao, L.
2012-07-01
More and more MOX fuels are used in all over the world in the past several decades. Compared with UO{sub 2} fuel, it contains some new features. For example, the neutron spectrum is harder and more resonance interference effects within the resonance energy range are introduced because of more resonant nuclides contained in the MOX fuel. In this paper, the wavelets scaling function expansion method is applied to study the resonance behavior of plutonium isotopes within MOX fuel. Wavelets scaling function expansion continuous-energy self-shielding method is developed recently. It has been validated and verified by comparison to Monte Carlo calculations. In this method, the continuous-energy cross-sections are utilized within resonance energy, which means that it's capable to solve problems with serious resonance interference effects without iteration calculations. Therefore, this method adapts to treat the MOX fuel resonance calculation problem natively. Furthermore, plutonium isotopes have fierce oscillations of total cross-section within thermal energy range, especially for {sup 240}Pu and {sup 242}Pu. To take thermal resonance effect of plutonium isotopes into consideration the wavelet scaling function expansion continuous-energy resonance calculation code WAVERESON is enhanced by applying the free gas scattering kernel to obtain the continuous-energy scattering source within thermal energy range (2.1 eV to 4.0 eV) contrasting against the resonance energy range in which the elastic scattering kernel is utilized. Finally, all of the calculation results of WAVERESON are compared with MCNP calculation. (authors)
CCSD(T)/CBS fragment-based calculations of lattice energy of molecular crystals.
Červinka, Ctirad; Fulem, Michal; Růžička, Květoslav
2016-02-14
A comparative study of the lattice energy calculations for a data set of 25 molecular crystals is performed using an additive scheme based on the individual energies of up to four-body interactions calculated using the coupled clusters with iterative treatment of single and double excitations and perturbative triples correction (CCSD(T)) with an estimated complete basis set (CBS) description. The CCSD(T)/CBS values on lattice energies are used to estimate sublimation enthalpies which are compared with critically assessed and thermodynamically consistent experimental values. The average absolute percentage deviation of calculated sublimation enthalpies from experimental values amounts to 13% (corresponding to 4.8 kJ mol(-1) on absolute scale) with unbiased distribution of positive to negative deviations. As pair interaction energies present a dominant contribution to the lattice energy and CCSD(T)/CBS calculations still remain computationally costly, benchmark calculations of pair interaction energies defined by crystal parameters involving 17 levels of theory, including recently developed methods with local and explicit treatment of electronic correlation, such as LCC and LCC-F12, are also presented. Locally and explicitly correlated methods are found to be computationally effective and reliable methods enabling the application of fragment-based methods for larger systems. PMID:26874495
CCSD(T)/CBS fragment-based calculations of lattice energy of molecular crystals
NASA Astrophysics Data System (ADS)
Červinka, Ctirad; Fulem, Michal; Růžička, Květoslav
2016-02-01
A comparative study of the lattice energy calculations for a data set of 25 molecular crystals is performed using an additive scheme based on the individual energies of up to four-body interactions calculated using the coupled clusters with iterative treatment of single and double excitations and perturbative triples correction (CCSD(T)) with an estimated complete basis set (CBS) description. The CCSD(T)/CBS values on lattice energies are used to estimate sublimation enthalpies which are compared with critically assessed and thermodynamically consistent experimental values. The average absolute percentage deviation of calculated sublimation enthalpies from experimental values amounts to 13% (corresponding to 4.8 kJ mol-1 on absolute scale) with unbiased distribution of positive to negative deviations. As pair interaction energies present a dominant contribution to the lattice energy and CCSD(T)/CBS calculations still remain computationally costly, benchmark calculations of pair interaction energies defined by crystal parameters involving 17 levels of theory, including recently developed methods with local and explicit treatment of electronic correlation, such as LCC and LCC-F12, are also presented. Locally and explicitly correlated methods are found to be computationally effective and reliable methods enabling the application of fragment-based methods for larger systems.
NASA Astrophysics Data System (ADS)
Packard, T. T.; Osma, N.; Fernández-Urruzola, I.; Codispoti, L. A.; Christensen, J. P.; Gómez, M.
2015-05-01
Oceanic depth profiles of plankton respiration are described by a power function, RCO2 = (RCO2)0 (z/z
NASA Astrophysics Data System (ADS)
Packard, T. T.; Osma, N.; Fernández-Urruzola, I.; Codispoti, L. A.; Christensen, J. P.; Gómez, M.
2014-11-01
Oceanic depth profiles of plankton respiration are described by a power function, RCO2 = (RCO2)0(z/z0)b similar to the vertical carbon flux profile. Furthermore, because both ocean processes are closely related, conceptually and mathematically, each can be calculated from the other. The exponent (b), always negative, defines the maximum curvature of the respiration depth-profile and controls the carbon flux. When b is large, the C flux (FC) from the epipelagic ocean is low and the nutrient retention efficiency (NRE) is high allowing these waters to maintain high productivity. The opposite occurs when b is small. This means that the attenuation of respiration in ocean water columns is critical in understanding and predicting both vertical FC as well as the capacity of epipelagic ecosystems to retain their nutrients. The NRE is a new metric defined as the ratio of nutrient regeneration in a seawater layer to the nutrients introduced into that layer via FC. A depth-profile of FC is the integral of water column respiration. This relationship facilitates calculating ocean sections of FC from water column respiration. In a FC section across the Peru upwelling system we found a FC maximum extending down to 400 m, 50 km off the Peru coast. Finally, coupling respiratory electron transport system activity to heterotrophic oxidative phosphorylation promoted the calculation of an ocean section of heterotrophic energy production (HEP). It ranged from 250 to 500 J d-1 m-3 in the euphotic zone, to less than 5 J d-1 m-3 below 200 m on this ocean section.
Complete calculation of electroweak corrections for polarized Møller scattering at high energies
NASA Astrophysics Data System (ADS)
Zykunov, V. A.
2009-09-01
A complete calculation of electroweak radiative corrections to observables of polarized Møller scattering at high energies was performed. This calculation took explicitly into account contributions caused by hard bremsstrahlung. A FORTRAN code that permitted including radiative corrections to high-energy Møller scattering under arbitrary electron-detection conditions was written. It was shown that the electroweak corrections caused by hard bremsstrahlung were rather strongly dependent on the choice of experimental cuts and changed substantially the polarization asymmetry in the region of high energies and over a broad interval of scattering angles.
ERIC Educational Resources Information Center
Wai, C. M.; Hutchinson, S. G.
1989-01-01
Discusses the calculation of free energy in reactions between silicon dioxide and carbon. Describes several computer programs for calculating the free energy minimization and their uses in chemistry classrooms. Lists 16 references. (YP)
McKechnie, Scott; Booth, George H.; Cohen, Aron J.; Cole, Jacqueline M.
2015-05-21
The best practice in computational methods for determining vertical ionization energies (VIEs) is assessed, via reference to experimentally determined VIEs that are corroborated by highly accurate coupled-cluster calculations. These reference values are used to benchmark the performance of density functional theory (DFT) and wave function methods: Hartree-Fock theory, second-order Møller-Plesset perturbation theory, and Electron Propagator Theory (EPT). The core test set consists of 147 small molecules. An extended set of six larger molecules, from benzene to hexacene, is also considered to investigate the dependence of the results on molecule size. The closest agreement with experiment is found for ionization energies obtained from total energy difference calculations. In particular, DFT calculations using exchange-correlation functionals with either a large amount of exact exchange or long-range correction perform best. The results from these functionals are also the least sensitive to an increase in molecule size. In general, ionization energies calculated directly from the orbital energies of the neutral species are less accurate and more sensitive to an increase in molecule size. For the single-calculation approach, the EPT calculations are in closest agreement for both sets of molecules. For the orbital energies from DFT functionals, only those with long-range correction give quantitative agreement with dramatic failing for all other functionals considered. The results offer a practical hierarchy of approximations for the calculation of vertical ionization energies. In addition, the experimental and computational reference values can be used as a standardized set of benchmarks, against which other approximate methods can be compared.
McKechnie, Scott; Booth, George H; Cohen, Aron J; Cole, Jacqueline M
2015-05-21
The best practice in computational methods for determining vertical ionization energies (VIEs) is assessed, via reference to experimentally determined VIEs that are corroborated by highly accurate coupled-cluster calculations. These reference values are used to benchmark the performance of density functional theory (DFT) and wave function methods: Hartree-Fock theory, second-order Møller-Plesset perturbation theory, and Electron Propagator Theory (EPT). The core test set consists of 147 small molecules. An extended set of six larger molecules, from benzene to hexacene, is also considered to investigate the dependence of the results on molecule size. The closest agreement with experiment is found for ionization energies obtained from total energy difference calculations. In particular, DFT calculations using exchange-correlation functionals with either a large amount of exact exchange or long-range correction perform best. The results from these functionals are also the least sensitive to an increase in molecule size. In general, ionization energies calculated directly from the orbital energies of the neutral species are less accurate and more sensitive to an increase in molecule size. For the single-calculation approach, the EPT calculations are in closest agreement for both sets of molecules. For the orbital energies from DFT functionals, only those with long-range correction give quantitative agreement with dramatic failing for all other functionals considered. The results offer a practical hierarchy of approximations for the calculation of vertical ionization energies. In addition, the experimental and computational reference values can be used as a standardized set of benchmarks, against which other approximate methods can be compared. PMID:26001454
Surface Segregation Energies of BCC Binaries from Ab Initio and Quantum Approximate Calculations
NASA Technical Reports Server (NTRS)
Good, Brian S.
2003-01-01
We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy method. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameterization. Quantum approximate segregation energies are computed with and without atomistic relaxation. The ab initio calculations are performed without relaxation for the most part, but predicted relaxations from quantum approximate calculations are used in selected cases to compute approximate relaxed ab initio segregation energies. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with other quantum approximate and ab initio theoretical work, and available experimental results.
Li, Z.; Pan, Y.K.; Tao, F.M.
1996-01-15
Bond function basis sets combined with the counterpoise procedure are used to calculate the molecular dissociation energies D{sub e} of 24 diatomic molecules and ions. The calculated values of D{sub e} are compared to those without bond functions and/or counterpoise corrections. The equilibrium bond lengths r{sub e}, and harmonic frequencies w{sub e} are also calculated for a few selected molecules. The calculations at the fourth-order-Moller-Plesset approximation (MP4) have consistently recovered about 95-99% of the experimental values for D{sub e}, compared to as low as 75% without use of bond functions. The calculated values of r{sub 3} are typically 0.01 {Angstrom} larger than the experimental values, and the calculated values of w{sub e} are over 95% of the experimental values. 37 refs., 2 tabs.
Solar Energy Education. Renewable energy activities for biology
Not Available
1982-01-01
An instructional aid for teachers is presented that will allow biology students the opportunity to learn about renewable energy sources. Some of the school activities include using leaves as collectors of solar energy, solar energy stored in wood, and a fuel value test for green and dry woods. A study of organic wastes as a source of fuel is included. (BCS)
Solar energy education. Renewable energy activities for general science
Not Available
1985-01-01
Renewable energy topics are integrated with the study of general science. The literature is provided in the form of a teaching manual and includes such topics as passive solar homes, siting a home for solar energy, and wind power for the home. Other energy topics are explored through library research activities. (BCS)
NASA Astrophysics Data System (ADS)
Li, Ming; Kapusta, Joseph I.
2016-08-01
We generalize calculations of the energy-momentum tensor for classical gluon fields in the boost-invariant McLerran-Venugopalan model using the small-τ power series expansion method. Results to all orders for the energy density and pressures are given in the leading Q2 approximation and with the inclusion of estimated running coupling effects. The energy density and transverse pressure decrease monotonically with time while the longitudinal pressure starts from a negative value and increases towards zero.
Validation of a new formula for calculating the energy requirements of burn patients.
Allard, J P; Pichard, C; Hoshino, E; Stechison, S; Fareholm, L; Peters, W J; Jeejeebhoy, K N
1990-01-01
In a previous study, we analyzed the roles of the estimated basal energy expenditure (EBEE) calculated by the Harris-Benedict equation, caloric intake (CI), total burn surface area (TBSA), body temperature (Temp) and number of post-burn days (PBD) in order to estimate the resting energy expenditure (EE) of burn patients. By multiple regression analysis we found that the measured EE (MEE) is best approximated by the following formula: -4343 + (10.5 x %TBSA) + (0.23 x CI) + (0.84 x EBEE) + (114 x Temp (degree C)) - (4.5 x PBD), r = 0.82, p less than 0.001, (Toronto formula (TF)). To validate this, 10 patients with a mean TBSA of 49.1 +/- 5.5% had their resting MEE done by indirect calorimetry when fed either by the TF or by 2 x EBEE. The caloric intake when on 2 x EBEE was 3260 +/- 45 kcal/day which significantly exceeded the MEE (2765 +/- 101 kcal/day, p less than 0.001). The caloric intake when on the TF was 2542 +/- 52 kcal/day and matches the MEE which was 2537 +/- 86 kcal/day (NS). These results show that the TF matches the MEE very closely. With the addition of a factor of activity for the 24-hr EE, it can be used to accurately feed individual burn patients. PMID:2112617
Energy-efficiency testing activities of the Mobile Energy Laboratory
Parker, G.B.
1991-01-01
This report summarizes energy-efficiency testing activities during the first and second quarters of fiscal year 1990 applying the Mobile Energy Laboratory (MEL) testing capabilities. Four MELs, developed by the US Department of Energy (DOE) Federal Energy Management Program (FEMP), are administered by Pacific Northwest Laboratory (PNL) for energy testing and program support functions at federal facilities. The using agencies principally fund MEL applications, while DOE/FEMP funds program administration and capability enhancement activities. This report fulfills the requirements established in the MEL Use Plan (PNL-6861) for semiannual reporting on energy-efficiency testing activities using the MEL capabilities. The MEL Use Committee, formally established in 1989, developed the MEL Use Plan and meets semiannually to establish priorities for energy-efficient testing applications using the MEL capabilities.
Monte Carlo calculation of the energy response characteristics of a RadFET radiation detector
NASA Astrophysics Data System (ADS)
Belicev, P.; Spasic Jokic, V.; Mayer, S.; Milosevic, M.; Ilic, R.; Pesic, M.
2010-07-01
The Metal -Oxide Semiconductor Field-Effect-Transistor (MOSFET, RadFET) is frequently used as a sensor of ionizing radiation in nuclear-medicine, diagnostic-radiology, radiotherapy quality-assurance and in the nuclear and space industries. We focused our investigations on calculating the energy response of a p-type RadFET to low-energy photons in range from 12 keV to 2 MeV and on understanding the influence of uncertainties in the composition and geometry of the device in calculating the energy response function. All results were normalized to unit air kerma incident on the RadFET for incident photon energy of 1.1 MeV. The calculations of the energy response characteristics of a RadFET radiation detector were performed via Monte Carlo simulations using the MCNPX code and for a limited number of incident photon energies the FOTELP code was also used for the sake of comparison. The geometry of the RadFET was modeled as a simple stack of appropriate materials. Our goal was to obtain results with statistical uncertainties better than 1% (fulfilled in MCNPX calculations for all incident energies which resulted in simulations with 1 - 2×109 histories.
Ghosh, Soumya; Hammes-Schiffer, Sharon
2015-01-02
Electrochemical electron transfer reactions play an important role in energy conversion processes with many technological applications. Electrodes modified by self-assembled monolayers (SAMs) are useful because the double layer effects are reduced. An important quantity for calculating the electron transfer rate constant is the reorganization energy, which is associated with changes in solute geometry and solvent configuration. In this Letter, an approach for calculating the electrochemical solvent reorganization energy for a redox molecule attached to or near a SAM modified electrode is presented. This integral equations formalism polarizable continuum model (IEF-PCM) approach accounts for the detailed electronic structure of the molecule, as well as the contributions from the electrode, SAM, and electronic and inertial solvent responses. The calculated total reorganization energies are in good agreement with experimental data for a series of metal complex in aqueous solution. This approach will be useful for calculating electron transfer rate constants for molecular electrocatalysts. This work was supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences.
Calculation of the characteristics of clinical high-energy photon beams with EGS5-MPI
NASA Astrophysics Data System (ADS)
Shimizu, M.; Morishita, Y.; Kato, M.; Kurosawa, T.; Tanaka, T.; Takata, N.; Saito, N.
2014-03-01
A graphite calorimeter has been developed as a Japanese primary standard of absorbed dose to water in the high-energy photon beams from a clinical linac. To obtain conversion factors for the graphite calorimeter, the beam characteristics of the high-energy photon beams from the clinical linac at National Metrology Institute of Japan were calculated with the EGS5 Monte Carlo simulation code. To run the EGS5 code on High Performance Computing machines that have more than 1000 CPU cores, we developed the EGS5 parallelisation package "EGS5-MPI" by implementing a message-passing interface. We calculated the photon energy spectra, which are in good agreement with those previously calculated by D. Sheikh-Bagheri and D. W. O. Rogers (Med. Phys. 29 3). We also estimated the percentage-depth-dose distributions of photon beams from the linac using the calculated photon energy spectra. These calculated percentage-depth-dose distributions were compared with our measured distributions and were found they are in good agreement as well. We will calculate conversion factors for the graphite calorimeter using our results.
Comparison between calculation and measurement of energy deposited by 800 MeV protons
Loewe, W.E.
1980-04-03
The High Energy Transport Code, HETC, was obtained from the Radiation Shielding Information Center (RSIC) at Oak Ridge National Laboratory and altered as necessary to run on a CDC 7600 using the LTSS software in use at LLNL. HETC was then used to obtain calculated estimates of energy deposited, for comparison with a series of benchmark experiments done by LLNL. These experiments used proton beams of various energies incident on well-defined composite targets in good geometry. In this report, two aspects of the comparison between calculated and experimental energy depositions from an 800 MeV proton beam are discussed. Both aspects involve the fact that workers at SAI had previously used their version of HETC to calculate this experiment and reported their comparison with the measured data. The first aspect addressed is that their calculated data and LLNL calculations do not agree, suggesting an error in the conversion process from the RSIC code. The second aspect is not independent of the first, but is of sufficient importance to merit separate emphasis. It is that the SAI calculations agree well with experiments at the detector plate located some distance from the shower plate, whereas the LLNL calculations show a clearcut discrepancy there in comparison with the experiment. A contract was let in January 1980 by LLNL with SAI in order to obtain full details on the two cited aspects of the comparison between calculated and experimental energy depositions from an 800 MeV proton beam. The ensuing discussion is based on the final report of that contracted work.
Solar Energy Project, Activities: General Solar Topics.
ERIC Educational Resources Information Center
Tullock, Bruce, Ed.; And Others
This guide contains lesson plans and outlines of activities which introduce students to concepts and issues relating to solar energy. Lessons frequently presented in the context of solar energy as it relates to contemporary energy problems. Each unit presents an introduction; objectives; necessary skills and knowledge; materials; method;…
AlaScan: A Graphical User Interface for Alanine Scanning Free-Energy Calculations.
Ramadoss, Vijayaraj; Dehez, François; Chipot, Christophe
2016-06-27
Computation of the free-energy changes that underlie molecular recognition and association has gained significant importance due to its considerable potential in drug discovery. The massive increase of computational power in recent years substantiates the application of more accurate theoretical methods for the calculation of binding free energies. The impact of such advances is the application of parent approaches, like computational alanine scanning, to investigate in silico the effect of amino-acid replacement in protein-ligand and protein-protein complexes, or probe the thermostability of individual proteins. Because human effort represents a significant cost that precludes the routine use of this form of free-energy calculations, minimizing manual intervention constitutes a stringent prerequisite for any such systematic computation. With this objective in mind, we propose a new plug-in, referred to as AlaScan, developed within the popular visualization program VMD to automate the major steps in alanine-scanning calculations, employing free-energy perturbation as implemented in the widely used molecular dynamics code NAMD. The AlaScan plug-in can be utilized upstream, to prepare input files for selected alanine mutations. It can also be utilized downstream to perform the analysis of different alanine-scanning calculations and to report the free-energy estimates in a user-friendly graphical user interface, allowing favorable mutations to be identified at a glance. The plug-in also assists the end-user in assessing the reliability of the calculation through rapid visual inspection. PMID:27214306
Hep, J.; Konecna, A.; Krysl, V.; Smutny, V.
2011-07-01
This paper describes the application of effective source in forward calculations and the adjoint method to the solution of fast neutron fluence and activation detector activities in the reactor pressure vessel (RPV) and RPV cavity of a VVER-440 reactor. Its objective is the demonstration of both methods on a practical task. The effective source method applies the Boltzmann transport operator to time integrated source data in order to obtain neutron fluence and detector activities. By weighting the source data by time dependent decay of the detector activity, the result of the calculation is the detector activity. Alternatively, if the weighting is uniform with respect to time, the result is the fluence. The approach works because of the inherent linearity of radiation transport in non-multiplying time-invariant media. Integrated in this way, the source data are referred to as the effective source. The effective source in the forward calculations method thereby enables the analyst to replace numerous intensive transport calculations with a single transport calculation in which the time dependence and magnitude of the source are correctly represented. In this work, the effective source method has been expanded slightly in the following way: neutron source data were performed with few group method calculation using the active core calculation code MOBY-DICK. The follow-up neutron transport calculation was performed using the neutron transport code TORT to perform multigroup calculations. For comparison, an alternative method of calculation has been used based upon adjoint functions of the Boltzmann transport equation. Calculation of the three-dimensional (3-D) adjoint function for each required computational outcome has been obtained using the deterministic code TORT and the cross section library BGL440. Adjoint functions appropriate to the required fast neutron flux density and neutron reaction rates have been calculated for several significant points within the RPV
NASA Astrophysics Data System (ADS)
Dieterich, Johannes M.; Werner, Hans-Joachim; Mata, Ricardo A.; Metz, Sebastian; Thiel, Walter
2010-01-01
Energy and free energy barriers for acetaldehyde conversion in aldehyde oxidoreductase are determined for three reaction pathways using quantum mechanical/molecular mechanical (QM/MM) calculations on the solvated enzyme. Ab initio single-point QM/MM energies are obtained at the stationary points optimized at the DFT(B3LYP)/MM level. These ab initio calculations employ local correlation treatments [LMP2 and LCCSD(T0)] in combination with augmented triple- and quadruple-zeta basis sets, and the final coupled cluster results include MP2-based corrections for basis set incompleteness and for the domain approximation. Free energy perturbation (FEP) theory is used to generate free energy profiles at the DFT(B3LYP)/MM level for the most important reaction steps by sampling along the corresponding reaction paths using molecular dynamics. The ab initio and FEP QM/MM results are combined to derive improved estimates of the free energy barriers, which differ from the corresponding DFT(B3LYP)/MM energy barriers by about 3 kcal mol-1. The present results confirm the qualitative mechanistic conclusions from a previous DFT(B3LYP)/MM study. Most favorable is a three-step Lewis base catalyzed mechanism with an initial proton transfer from the cofactor to the Glu869 residue, a subsequent nucleophilic attack that yields a tetrahedral intermediate (IM2), and a final rate-limiting hydride transfer. The competing metal center activated pathway has the same final step but needs to overcome a higher barrier in the initial step on the route to IM2. The concerted mechanism has the highest free energy barrier and can be ruled out. While confirming the qualitative mechanistic scenario proposed previously on the basis of DFT(B3LYP)/MM energy profiles, the present ab initio and FEP QM/MM calculations provide corrections to the barriers that are important when aiming at high accuracy.
Scaling equations to model variations in operating conditions for activation calculations
Ho, S.K.
1994-11-01
A simple approximate scheme of calculating the variations of fusion activation-induced radioactive inventories with reactor operating conditions (power and time) by scaling equations is formulated. Application of the scheme to evaluate radiological dose-contributing activation products for HT-9 ferritic steel first wall shows decent agreements with numerical results.
NASA Astrophysics Data System (ADS)
Li, Yongxiu; Zhang, Saiqun; Zhang, John Z. H.; He, Xiao
2016-05-01
Accurate description of the conformational energies of the amino acids is essential for molecular dynamics simulation of protein structures. In this study, we compute the relative energies at 51 conformations for a trialanine tetrapeptide at different levels of theory. The computed energies at various theoretical levels, including the semiempirical DFTB method, HF, DFT, MP2 and CCSD(T), are compared with each other. The calculated energies from density-fitting local CCSD(T)/CBS (complete basis set) calculations are taken as the benchmark. The accuracy of the theoretical methods is highly dependent on the electronic correlation and dispersion corrections as well as the size of the basis sets. The involvement of the empirical dispersion energies in HF and DFT methods consistently improves their performance. Considering both the accuracy and computational efficiency, the Minnesota density functional M06-L-D and M06-2X-D are efficient and accurate for modeling of trialanine structures.
Gaylord, R F
2005-07-26
In order to determine if the measured beta activity for a solution containing potassium was exactly as predicted, particularly since the CES gas counter is not calibrated specifically with K-40, an experiment was conducted to compare measured activities from two radioanalytical methods (gamma spectroscopy and gas proportional counting) to calculated activities across a range of potassium concentrations. Potassium, being ubiquitous and naturally radioactive, is a well-known and common interference in gross beta counting methods. By measuring the observed beta activity due to K-40 in potassium-containing solutions across a wide range of concentrations, it was found that the observed beta activity agrees well with the beta activity calculated from the potassium concentration measured by standard inorganic analytical techniques, such as ICP-OES, and that using the measured potassium concentration to calculate the expected beta activity, and comparing this to the observed beta activity to determine if potassium can account for all the observed activity in a sample, is a valid technique. It was also observed that gamma spectroscopy is not an effective means of measuring K-40 activity below approximately 700 pCi/L, which corresponds to a solution with approximately 833 mg/L total potassium. Gas proportional counting for gross beta activity has a much lower detection limit, typically 20-50 picoCi/L for a liquid low in total dissolved solids, which corresponds to a potassium concentration of approximately 30-70 ppm K.
NASA Astrophysics Data System (ADS)
Shi, De-Heng; Liu, Yu-Fang; Sun, Jin-Feng; Yang, Xiang-Dong; Zhu, Zun-Lue
2006-05-01
The reasonable dissociation limit of the A1Σ+ state 7LiH molecule is obtained. The accurate dissociation energy and the equilibrium geometry of this state are calculated using a symmetry-adapted-cluster configuration-interaction method in complete active space for the first time. The whole potential energy curve and the dipole moment function for the A1Σ+ state are calculated over a wide internuclear separation range from about 0.1 to 1.4 nm. The calculated equilibrium geometry and dissociation energy of this potential energy curve are of Re=0.2487 nm and De=1.064 eV, respectively. The unusual negative values of the anharmonicity constant and the vibration-rotational coupling constant are of ωeχe=-4.7158cm-1 and αe=-0.08649cm-1, respectively. The vertical excitation energy from the ground to the A1Σ+ state is calculated and the value is of 3.613 eV at 0.15875 nm (the equilibrium position of the ground state). The highly anomalous shape of this potential energy curve, which is exceptionally flat over a wide radial range around the equilibrium position, is discussed in detail. The harmonic frequency value of 502.47cm-1 about this state is approximately estimated. Careful comparison of the theoretical determinations with those obtained by previous theories about the A1Σ+ state dissociation energy clearly shows that the present calculations are much closer to the experiments than previous theories, thus represents an improvement.
Hybrid energy harvesting using active thermal backplane
NASA Astrophysics Data System (ADS)
Kim, Hyun-Wook; Lee, Dong-Gun
2016-04-01
In this study, we demonstrate the concept of a new hybrid energy harvesting system by combing solar cells with magneto-thermoelectric generator (MTG, i.e., thermal energy harvesting). The silicon solar cell can easily reach high temperature under normal operating conditions. Thus the heated solar cell becomes rapidly less efficient as the temperature of solar cell rises. To increase the efficiency of the solar cell, air or water-based cooling system is used. To surpass conventional cooling devices requiring additional power as well as large working space for air/water collectors, we develop a new technology of pairing an active thermal backplane (ATB) to solar cell. The ATB design is based on MTG technology utilizing the physics of the 2nd order phase transition of active ferromagnetic materials. The MTG is cost-effective conversion of thermal energy to electrical energy and is fundamentally different from Seebeck TEG devices. The ATB (MTG) is in addition to being an energy conversion system, a very good conveyor of heat through both conduction and convection. Therefore, the ATB can provide dual-mode for the proposed hybrid energy harvesting. One is active convective and conductive cooling for heated solar cell. Another is active thermal energy harvesting from heat of solar cell. These novel hybrid energy harvesting device have potentially simultaneous energy conversion capability of solar and thermal energy into electricity. The results presented can be used for better understanding of hybrid energy harvesting system that can be integrated into commercial applications.
Energy and Energy Conservation Activities for High School Students.
ERIC Educational Resources Information Center
Bottinelli, Charles A., Ed.; Dow, John O., Ed.
This manual contains fifteen energy activities suitable for high school physical and environmental science and mathematics classrooms. The activities are independent, each having its own objectives, introduction, and background information. A special section of each activity is written for the instructor and contains limits, sample data, and…
NASA Astrophysics Data System (ADS)
Tan, Zhenyu; Dong, Lei; Tang, Fule
2012-08-01
The calculations of the characteristic quantities of low-energy electron (⩽20 keV) irradiation to the five typical spacecraft dielectrics, i.e. epoxy, kapton, mylar, polyethylene, and teflon, have been performed by means of Monte Carlo method. These characteristic quantities include the electron backscattering coefficient, the depth distributions of both energy deposition and deposited electrons, and the maximum penetration depth of deposited electrons in the dielectrics. A Monte Carlo model has been specifically constructed for simulating the transport of low-energy electrons in spacecraft dielectrics (organic materials). In this model, the description of the inelastic scattering of energetic electron is based on the dielectric approach developed previously and the Born-Ochkur's exchange correction is included. Especially, the optical energy loss functions of organic materials can be obtained using an empirical evaluation. In addition, the mean cross section based on the Mott model is proposed for calculating the elastic scattering of energetic electrons in organic materials for high simulation efficiency. The constructed Monte Carlo model has been examined by a series of calculations and comparisons with the reported experiments and other theoretical results. For the dielectrics under consideration and in the energy range of E0 ⩽ 20 keV, the calculated electron backscattering coefficients and the extrapolated range of deposited electrons are listed at selected energies in numerical form for convenient use, and an empirical expression of estimating the extrapolated range in the energy range of 1-20 keV is given. The distribution characteristics of both energy deposition and deposited electrons are presented, and it is found that kapton and mylar present the close characteristic quantities for each other, which is of significance for the choice of the dielectrics in design of spacecraft. The characteristic quantity calculations presented in this work are a
Calculation of intensity of high energy muon groups observed deep underground
NASA Technical Reports Server (NTRS)
Vavilov, Y. N.; Dedenko, L. G.
1985-01-01
The intensity of narrow muon groups observed in Kolar Gold Field (KGF) at the depth of 3375 m.w.e. was calculated in terms of quark-gluon strings model for high energy hadron - air nuclei interactions by the method of direct modeling of nuclear cascade in the air and muon propagation in the ground for normal primary cosmic ray composition. The calculated intensity has been found to be approx. 10 to the 4 times less than one observed experimentally.