Relativistic Dirac-Fock-Slater program to calculate potential-energy curves for diatomic molecules

A LCAO-MO (linear combination of atomic orbitals - molecular orbitals) relativistic Dirac-Fock-Slater program is presented, which allows one to calculate accurate total energies for diatomic molecules. Numerical atomic Dirac-Fock-Slater wave functions are used as basis functions. All integrations as well as the solution of the Poisson equation are done fully numerical, with a relative accuracy of 10{^-5} - 10{^-6}. The details of the method as well as first results are presented here.

Asymptotically Zero Energy Computing Using Split-Level Charge Recovery Logic

The dynamic power requirement of CMOS circuits is rapidly becoming a major concern in the design of personal information systems and large computers. In this work we present a number of new CMOS logic families, Charge Recovery Logic (CRL) as well as the much improved Split-Level Charge Recovery Logic (SCRL), within which the transfer of charge between the nodes occurs quasistatically. Operating quasistatically, these logic families have an energy dissipation that drops linearly with operating frequency, i.e., their power consumption drops quadratically with operating frequency as opposed to the linear drop of conventional CMOS. The circuit techniques in these new families rely on constructing an explicitly reversible pipelined logic gate, where the information necessary to recover the energy used to compute a value is provided by computing its logical inverse. Information necessary to uncompute the inverse is available from the subsequent inverse logic stage. We demonstrate the low energy operation of SCRL by presenting the results from the testing of the first fully quasistatic 8 x 8 multiplier chip (SCRL-1) employing SCRL circuit techniques.

Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients BK of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, ΔEK=(E-HKK)BK2/(1-BK2), with BK determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Δdis, is approximated by the sum of ΔEKS of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound ES is computed by CI in a selected space S, and EM=E S+ΔEdis+δE, where δE is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm-1) is achieved in a model space M of 1.4× 109 CSFs (1.1 × 1012 determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 × 1012 CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since ΔEdis can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of ES is taken up in a companion paper

Variational calculation of vibrational linear and nonlinear optical properties

A variational approach for reliably calculating vibrational linear and nonlinear optical properties of molecules with large electrical and/or mechanical anharmonicity is introduced. This approach utilizes a self-consistent solution of the vibrational Schrödinger equation for the complete field-dependent potential-energy surface and, then, adds higher-level vibrational correlation corrections as desired. An initial application is made to static properties for three molecules of widely varying anharmonicity using the lowest-level vibrational correlation treatment (i.e., vibrational Møller-Plesset perturbation theory). Our results indicate when the conventional Bishop-Kirtman perturbation method can be expected to break down and when high-level vibrational correlation methods are likely to be required. Future improvements and extensions are discussed

Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential VA on a one-atom density ρ A or ρ B. To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale

Potential energy surface intersections in triatomic molecules

The different types of surface intersection which may occur in linear configurations of triatomic molecules are reviewed, particularly with regard to the way in which the degeneracy is split as the molecule bends. The Renner-Teller effect in states of symmetry Π, Δ, Φ, etc., and intersections between Σ and Π, Σ and Δ, and Π and Δ states are discussed. A general method of modelling such intersecting potential surfaces is proposed, as a development of the model previously used by Murrell and Carter and co-workers for single-valued surfaces. Some of the lower energy surfaces of H2O, NH2, O3, C3, and HNO are discussed as examples.

Effect of energy and protein supplementation on phosphorus utilization in lactating dairy cows

Two experiments were undertaken in which grass silage was used in conjunction with a series of different concentrate types designed to examine the effect of carbohydrate source, protein level and degradability on total dietary phosphorus (P) utilization with emphasis on P pollution. Twelve Holstein-Friesian dairy cows in early to mid-lactation were used in an incomplete changeover design with four periods consisting of 4 weeks each. Phosphorus intake ranged from 54 to 80 g/day and faecal P represented the principal route by which ingested P was disposed of by cows, with insignificant amounts being voided in urine. A positive linear relationship between faecal P and P intake was established. In Experiment 1, P utilization was affected by dietary carbohydrate type, with an associated output of 3.3 g faecal P/g milk P produced for all treatments except those utilizing low degradable starch and low protein supplements, where a mean value of 2.8 g faecal P/g milk P was observed. In Experiment 2, where two protein levels and three protein degradabilities were examined, the efficiency of P utilization for milk P production was not affected by either level or degradability of crude protein (CP) but a significant reduction in faecal P excretion due to lower protein and P intake was observed. In general, P utilization in Experiment 2 was substantially improved compared to the Experiment 1, with an associated output of 1.8 g faecal P/g milk P produced. The improved utilization of P in Experiment 2 could be due to lower P content of the diets offered and higher dry matter (DM) intake. For dairy cows weighing 600 kg, consuming 17-18 kg DM/day and producing about 25 kg milk, P excretion in faeces and hence P pollution to the environment might be minimized without compromising lactational performance by formulating diets to supply about 68 g P/day, which is close to recent published recommended requirements for P.

Structure and biogenesis of the capsular F1 antigen from Yersinia pestis: Preserved folding energy drives fiber formation

Most gram-negative pathogens express fibrous adhesive virulence organelles that mediate targeting to the sites of infection. The F1 capsular antigen from the plague pathogen Yersinia pestis consists of linear fibers of a single subunit (Caf1) and serves as a prototype for nonpilus organelles assembled via the chaperone/usher pathway. Genetic data together with high-resolution X-ray structures corresponding to snapshots of the assembly process reveal the structural basis of fiber formation. Comparison of chaperone bound Caf1 subunit with the subunit in the fiber reveals a novel type of conformational change involving the entire hydrophobic core of the protein. The observed conformational change suggests that the chaperone traps a high-energy folding intermediate of Caf1. A model is proposed in which release of the subunit allows folding to be completed, driving fiber formation.

A case of four-coordinate silver(I) adopting a high energy planar geometry. Experiment and theory

The ligands PhL and MeL are obtained by condensing 2-formylpyridine with benzil dihydrazone and diacetyl dihydrazone, respectively, in 2: 1 molar proportion. With silver( I), PhL yields a double-stranded dinuclear cationic helicate 1 in which the metal is tetrahedral but MeL gives a cationic one-dimensional polymeric complex 2 where silver( I) is distorted square planar and the ligand backbone is nearly planar. In both complexes, metal: ligand ratio is 1: 1. Ab initio calculations on the ligands at the HF/6-31+G* level reveal that while PhL strongly prefers a helical conformation, MeL has a natural inclination to remain in a planar conformation. Density functional theory calculations on model silver( I) complexes show that formation of the linear polymer in the case of MeL is also an important factor in imposing the planar geometry of Ag(I) in 2.

Potential energy surface and moleculardynamics of MbNO: existence of an unsuspected FeON minimum

Ligands such as CO, O2, or NO are involved in the biological function of myoglobin. Here we investigate the energetics and dynamics of NO interacting with the Fe(II) heme group in native myoglobin using ab initio and molecular dynamics simulations. At the global minimum of the ab initio potential energy surface (PES), the binding energy of 23.4 kcal/mol and the Fe-NO structure compare well with the experimental results. Interestingly, the PES is found to exhibit two minima: There exists a metastable, linear Fe-O-N minimum in addition to the known, bent Fe-N-O global minimum conformation. Moreover, the T-shaped configuration is found to be a saddle point, in contrast to the corresponding minimum for NO interacting with Fe(III). To use the ab initio results for finite temperature molecular dynamics simulations, an analytical function was fitted to represent the Fe-NO interaction. The simulations show that the secondary minimum is dynamically stable up to 250 K and has a lifetime of several hundred picoseconds at 300 K. The difference in the topology of the heme-NO PES from that assumed previously (one deep, single Fe-NO minimum) suggests that it is important to use the full PES for a quantitative understanding of this system. Why the metastable state has not been observed in the many spectroscopic studies of myoglobin interacting with NO is discussed, and possible approaches to finding it are outlined.

Energy and protein interactions and their effect on nitrogen excretion in dairy cows

The principal driver of nitrogen (N) losses from the body including excretion and secretion in milk is N intake. However, other covariates may also play a role in modifying the partitioning of N. This study tests the hypothesis that N partitioning in dairy cows is affected by energy and protein interactions. A database containing 470 dairy cow observations was collated from calorimetry experiments. The data include N and energy parameters of the diet and N utilization by the animal. Univariate and multivariate meta-analyses that considered both within and between study effects were conducted to generate prediction equations based on N intake alone or with an energy component. The univariate models showed that there was a strong positive linear relationships between N intake and N excretion in faeces, urine and milk. The slopes were 0.28 faeces N, 0.38 urine N and 0.20 milk N. Multivariate model analysis did not improve the fit. Metabolizable energy intake had a significant positive effect on the amount of milk N in proportion to faeces and urine N, which is also supported by other studies. Another measure of energy considered as a covariate to N intake was diet quality or metabolizability (the concentration of metabolizable energy relative to gross energy of the diet). Diet quality also had a positive linear relationship with the proportion of milk N relative to N excreted in faeces and urine. Metabolizability had the largest effect on faeces N due to lower protein digestibility of low quality diets. Urine N was also affected by diet quality and the magnitude of the effect was higher than for milk N. This research shows that including a measure of diet quality as a covariate with N intake in a model of N execration can enhance our understanding of the effects of diet composition on N losses from dairy cows. The new prediction equations developed in this study could be used to monitor N losses from dairy systems.

Further investigations of the quartz optically stimulated luminescence components using linear modulation

The optically stimulated luminescence (OSL) from quartz is known to be the sum of several components with different rates of charge loss, originating from different trap types. The OSL components are clearly distinguished using the linear modulation (LM OSL) technique. A variety of pre-treatment and measurement conditions have been used on sedimentary samples in conjunction with linearly modulated optical stimulation to study in detail the behaviour of the OSL components of quartz. Single aliquots of different quartz samples have been found to contain typically five or six common LM OSL components when stimulated at View the MathML source. The components have been parameterised in terms of thermal stability (i.e. E and s), photoionisation cross-section energy dependence and dose response. The results of studies concerning applications of component-resolved LM OSL measurements on quartz are also presented. These include the detection of partial bleaching in young samples, use of ‘stepped wavelength’ stimulation to observe OSL from single components and attempts to extend the age range of quartz OSL dating.

Short communication: plasma concentration of glucose-dependent insulinotropic polypeptide may regulate milk energy production in lactating dairy cows

In dairy cows, an increase in plasma concentration of glucose-dependent insulinotropic polypeptide (GIP) is associated with an increase in metabolizable energy intake, but the role of GIP in energy partitioning of dairy cattle is not certain. The objective of this study was to examine the relationship between plasma GIP concentrations and energy partitioning toward milk production. Four mid-lactation, primiparous, rumenfistulated Holstein-Friesian cows were fed a control diet of 55% forage and 45% concentrate [dry matter (DM) basis] in a 4 × 4 Latin square design with 4-wk periods. The 4 treatments were (1) control diet fed at 1000 and 1600 h, and (2) once-daily (1000 h) feeding, (3) twice daily (1000 and 1600 h) feeding, and (4) 4 times/d (1000, 1600, 2200 and 0400 h) feeding of the control diet plus 1 dose (1.75 kg on a DM basis at 0955 h) into the rumen of supplemental vegetable proteins (Amino Green; SCA NuTec Ltd., Thirsk, UK). Measurements of respiratory exchange and energy balance were obtained over 4 d during the last week of each period while cows were housed in open-circuit respiration chambers. Blood was collected from the jugular vein every 30 min for 12 h, using indwelling catheters, starting at 0800 h on d 20 of each period. Plasma GIP concentration was measured in samples pooled over each 5 consecutive blood samplings. The relationships between plasma GIP, DM intake, heat production, respiratory quotient, milk yield, and milk energy output were analyzed using linear correlation procedures, with metabolizable intake as a partial variant. Plasma GIP concentration was not correlated with heat production, or milk yield, but was positively correlated with milk energy yield (correlation coefficient = 0.67) and negatively correlated with RQ (correlation coefficient = −0.72). The correlations between GIP and RQ and milk energy output do not imply causality, but suggest that a role for GIP may exist in the regulation of energy metabolism in dairy cows.

Energy-balance mechanisms underlying consistent large-scale temperature responses in warm and cold climates

Climate simulations show consistent large-scale temperature responses including amplified land–ocean contrast, high-latitude/low-latitude contrast, and changes in seasonality in response to year-round forcing, in both warm and cold climates, and these responses are proportional and nearly linear across multiple climate states. We examine the possibility that a small set of common mechanisms controls these large-scale responses using a simple energy-balance model to decompose the temperature changes shown in multiple lgm and abrupt4 × CO 2 simulations from the CMIP5 archive. Changes in the individual components of the energy balance are broadly consistent across the models. Although several components are involved in the overall temperature responses, surface downward clear-sky longwave radiation is the most important component driving land–ocean contrast and high-latitude amplification in both warm and cold climates. Surface albedo also plays a significant role in promoting high-latitude amplification in both climates and in intensifying the land–ocean contrast in the warm climate case. The change in seasonality is a consequence of the changes in land–ocean and high-latitude/low-latitude contrasts rather than an independent temperature response. This is borne out by the fact that no single component stands out as being the major cause of the change in seasonality, and the relative importance of individual components is different in cold and warm climates.

A comparison of duality and energy a posteriori estimates for $\mathrm {L}_{\infty }(0,T;\mathrm {L}_2(\varOmega ))$ in parabolic problems

We use the elliptic reconstruction technique in combination with a duality approach to prove a posteriori error estimates for fully discrete backward Euler scheme for linear parabolic equations. As an application, we combine our result with the residual based estimators from the a posteriori estimation for elliptic problems to derive space-error indicators and thus a fully practical version of the estimators bounding the error in the $ \mathrm {L}_{\infty }(0,T;\mathrm {L}_2(\varOmega ))$ norm. These estimators, which are of optimal order, extend those introduced by Eriksson and Johnson in 1991 by taking into account the error induced by the mesh changes and allowing for a more flexible use of the elliptic estimators. For comparison with previous results we derive also an energy-based a posteriori estimate for the $ \mathrm {L}_{\infty }(0,T;\mathrm {L}_2(\varOmega ))$-error which simplifies a previous one given by Lakkis and Makridakis in 2006. We then compare both estimators (duality vs. energy) in practical situations and draw conclusions.