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

Simple finite field nuclear relaxation method for calculating vibrational contribution to degenerate four-wave mixing

A simple extended finite field nuclear relaxation procedure for calculating vibrational contributions to degenerate four-wave mixing (also known as the intensity-dependent refractive index) is presented. As a by-product one also obtains the static vibrationally averaged linear polarizability, as well as the first and second hyperpolarizability. The methodology is validated by illustrative calculations on the water molecule. Further possible extensions are suggested

Automatic generation of active coordinates for quantum dynamics calculations: application to the dynamics of benzene photochemistry

A new practical method to generate a subspace of active coordinates for quantum dynamics calculations is presented. These reduced coordinates are obtained as the normal modes of an analytical quadratic representation of the energy difference between excited and ground states within the complete active space self-consistent field method. At the Franck-Condon point, the largest negative eigenvalues of this Hessian correspond to the photoactive modes: those that reduce the energy difference and lead to the conical intersection; eigenvalues close to 0 correspond to bath modes, while modes with large positive eigenvalues are photoinactive vibrations, which increase the energy difference. The efficacy of quantum dynamics run in the subspace of the photoactive modes is illustrated with the photochemistry of benzene, where theoretical simulations are designed to assist optimal control experiments

Relaxation of selective constraints causes independent selenoprotein extinction in insect genomes

BACKGROUND: Selenoproteins are a diverse family of proteins notable for the presence of the 21st amino acid, selenocysteine. Until very recently, all metazoan genomes investigated encoded selenoproteins, and these proteins had therefore been believed to be essential for animal life. Challenging this assumption, recent comparative analyses of insect genomes have revealed that some insect genomes appear to have lost selenoprotein genes. METHODOLOGY/PRINCIPAL FINDINGS: In this paper we investigate in detail the fate of selenoproteins, and that of selenoprotein factors, in all available arthropod genomes. We use a variety of in silico comparative genomics approaches to look for known selenoprotein genes and factors involved in selenoprotein biosynthesis. We have found that five insect species have completely lost the ability to encode selenoproteins and that selenoprotein loss in these species, although so far confined to the Endopterygota infraclass, cannot be attributed to a single evolutionary event, but rather to multiple, independent events. Loss of selenoproteins and selenoprotein factors is usually coupled to the deletion of the entire no-longer functional genomic region, rather than to sequence degradation and consequent pseudogenisation. Such dynamics of gene extinction are consistent with the high rate of genome rearrangements observed in Drosophila. We have also found that, while many selenoprotein factors are concomitantly lost with the selenoproteins, others are present and conserved in all investigated genomes, irrespective of whether they code for selenoproteins or not, suggesting that they are involved in additional, non-selenoprotein related functions. CONCLUSIONS/SIGNIFICANCE: Selenoproteins have been independently lost in several insect species, possibly as a consequence of the relaxation in insects of the selective constraints acting across metazoans to maintain selenoproteins. The dispensability of selenoproteins in insects may be related to the fundamental differences in antioxidant defense between these animals and other metazoans.

On the Properties and Cosmology of Q-balls

The properties and cosmological importance of a class of non-topological solitons, Q-balls, are studied. Aspects of Q-ball solutions and Q-ball cosmology discussed in the literature are reviewed. Q-balls are particularly considered in the Minimal Supersymmetric Standard Model with supersymmetry broken by a hidden sector mechanism mediated by either gravity or gauge interactions. Q-ball profiles, charge-energy relations and evaporation rates for realistic Q-ball profiles are calculated for general polynomial potentials and for the gravity mediated scenario. In all of the cases, the evaporation rates are found to increase with decreasing charge. Q-ball collisions are studied by numerical means in the two supersymmetry breaking scenarios. It is noted that the collision processes can be divided into three types: fusion, charge transfer and elastic scattering. Cross-sections are calculated for the different types of processes in the different scenarios. The formation of Q-balls from the fragmentation of the Aflieck-Dine -condensate is studied by numerical and analytical means. The charge distribution is found to depend strongly on the initial energy-charge ratio of the condensate. The final state is typically noted to consist of Q- and anti-Q-balls in a state of maximum entropy. By studying the relaxation of excited Q-balls the rate at which excess energy can be emitted is calculated in the gravity mediated scenario. The Q-ball is also found to withstand excess energy well without significant charge loss. The possible cosmological consequences of these Q-ball properties are discussed.

Cumulative cultural dynamics and the coevolution of cultural innovation and transmission: an ESS model for panmictic and structured populations.

When individuals in a population can acquire traits through learning, each individual may express a certain number of distinct cultural traits. These traits may have been either invented by the individual himself or acquired from others in the population. Here, we develop a game theoretic model for the accumulation of cultural traits through individual and social learning. We explore how the rates of innovation, decay, and transmission of cultural traits affect the evolutionary stable (ES) levels of individual and social learning and the number of cultural traits expressed by an individual when cultural dynamics are at a steady-state. We explore the evolution of these phenotypes in both panmictic and structured population settings. Our results suggest that in panmictic populations, the ES level of learning and number of traits tend to be independent of the social transmission rate of cultural traits and is mainly affected by the innovation and decay rates. By contrast, in structured populations, where interactions occur between relatives, the ES level of learning and the number of traits per individual can be increased (relative to the panmictic case) and may then markedly depend on the transmission rate of cultural traits. This suggests that kin selection may be one additional solution to Rogers's paradox of nonadaptive culture.

Equivalence of piecewise-linear approximation and Lagrangian relaxation for network revenue management

The network revenue management (RM) problem arises in airline, hotel, media,and other industries where the sale products use multiple resources. It can be formulatedas a stochastic dynamic program but the dynamic program is computationallyintractable because of an exponentially large state space, and a number of heuristicshave been proposed to approximate it. Notable amongst these -both for their revenueperformance, as well as their theoretically sound basis- are approximate dynamic programmingmethods that approximate the value function by basis functions (both affinefunctions as well as piecewise-linear functions have been proposed for network RM)and decomposition methods that relax the constraints of the dynamic program to solvesimpler dynamic programs (such as the Lagrangian relaxation methods). In this paperwe show that these two seemingly distinct approaches coincide for the network RMdynamic program, i.e., the piecewise-linear approximation method and the Lagrangianrelaxation method are one and the same.

Crossing the Rio Grande: Migrations, business cycles and the welfare state

In this paper we study the macroeconomic effects of an inflow oflow-skilled workers into an economy where there is capital accumulation and two types of agents. We find that there are substantial dynamic effects following unexpected migrations with adjustments that resemble those triggered by a sudden disruption of the capital stock. We look at the interrelations between these dynamic effects and three different fiscal systems for the redistribution of income and find that these schemes can change the dynamics and lead to prolonged periods of adjustments. Theaggregate welfare implications are sensitive to the welfare system: while there are welfare gains without redistribution, these gains may be turned into costs when the state engages in redistribution.

Minimum price variations, time priority and quotes dynamics

We analyze the impact of a minimum price variation (tick) and timepriority on the dynamics of quotes and the trading costs when competitionfor the order flow is dynamic. We find that convergence to competitiveoutcomes can take time and that the speed of convergence is influencedby the tick size, the priority rule and the characteristics of the orderarrival process. We show also that a zero minimum price variation is neveroptimal when competition for the order flow is dynamic. We compare thetrading outcomes with and without time priority. Time priority is shownto guarantee that uncompetitive spreads cannot be sustained over time.However it can sometimes result in higher trading costs. Empiricalimplications are proposed. In particular, we relate the size of thetrading costs to the frequency of new offers and the dynamics of theinside spread to the state of the book.

Modelling state dependence and feedback effects between poverty, employment and parental home emancipation among European youth

Youth is one of the phases in the life-cycle when some of the most decisivelife transitions take place. Entering the labour market or leaving parentalhome are events with important consequences for the economic well-beingof young adults. In this paper, the interrelationship between employment,residential emancipation and poverty dynamics is studied for eight Europeancountries by means of an econometric model with feedback effects. Resultsshow that youth poverty genuine state dependence is positive and highly significant.Evidence proves there is a strong causal effect between poverty andleaving home in Scandinavian countries, however, time in economic hardshipdoes not last long. In Southern Europe, instead, youth tend to leave theirparental home much later in order to avoid falling into a poverty state that ismore persistent. Past poverty has negative consequences on the likelihood ofemployment.

Homeostatic and circadian contribution to EEG and molecular state variables of sleep regulation.

STUDY OBJECTIVES: Besides their well-established role in circadian rhythms, our findings that the forebrain expression of the clock-genes Per2 and Dbp increases and decreases, respectively, in relation to time spent awake suggest they also play a role in the homeostatic aspect of sleep regulation. Here, we determined whether time of day modulates the effects of elevated sleep pressure on clock-gene expression. Time of day effects were assessed also for recognized electrophysiological (EEG delta power) and molecular (Homer1a) markers of sleep homeostasis. DESIGN: EEG and qPCR data were obtained for baseline and recovery from 6-h sleep deprivation starting at ZT0, -6, -12, or -18. SETTING: Mouse sleep laboratory. PARTICIPANTS: Male mice. INTERVENTIONS: Sleep deprivation. RESULTS: The sleep-deprivation induced changes in Per2 and Dbp expression importantly varied with time of day, such that Per2 could even decrease during sleep deprivations occurring at the decreasing phase in baseline. Dbp showed similar, albeit opposite dynamics. These unexpected results could be reliably predicted assuming that these transcripts behave according to a driven damped harmonic oscillator. As expected, the sleep-wake distribution accounted for a large degree of the changes in EEG delta power and Homer1a. Nevertheless, the sleep deprivation-induced increase in delta power varied also with time of day with higher than expected levels when recovery sleep started at dark onset. CONCLUSIONS: Per2 and delta power are widely used as exclusive state variables of the circadian and homeostatic process, respectively. Our findings demonstrate a considerable cross-talk between these two processes. As Per2 in the brain responds to both sleep loss and time of day, this molecule is well positioned to keep track of and to anticipate homeostatic sleep need. CITATION: Curie T; Mongrain V; Dorsaz S; Mang GM; Emmenegger Y; Franken P. Homeostatic and circadian contribution to EEG and molecular state variables of sleep regulation. SLEEP 2013;36(3):311-323.

Dynamics of 14C-labeled glucose and ammonium in saline arable soils

Organic matter dynamics and nutrient availability in saline agricultural soils of the State of Guanajuato might provide information for remediation strategies. 14C labeled glucose with or without 200 mg kg-1 of NH4+-N soil was added to two clayey agricultural soils with different electrolytic conductivity (EC), i.e. 0.94 dS m-1 (low EC; LEC) and 6.72 dS m-1 (high EC; HEC), to investigate the effect of N availability and salt content on organic material decomposition. Inorganic N dynamics and production of CO2 and 14CO2 were monitored. Approximately 60 % of the glucose-14C added to LEC soil evolved as 14CO2, but only 20 % in HEC soil after the incubation period of 21 days. After one day, < 200 mg 14C was extractable from LEC soil, but > 500 mg 14C from HEC soil. No N mineralization occurred in the LEC and HEC soils and glucose addition reduced the concentrations of inorganic N in unamended soil and soil amended with NH4+-N. The NO2- and NO3- concentrations were on average higher in LEC than in HEC soil, with exception of NO2- in HEC amended with NH4+-N. It was concluded that increases in soil EC reduced mineralization of the easily decomposable C substrate and resulted in N-depleted soil.

Soil water dynamics related to the degree of compaction of two brazilian oxisols under no-tillage

Soil water properties are related to crop growth and environmental aspects and are influenced by the degree of soil compaction. The objective of this study was to determine the water infiltration and hydraulic conductivity of saturated soil under field conditions in terms of the compaction degree of two Oxisols under a no-tillage (NT). Two commercial fields were studied in the state of Rio Grande do Sul, Brazil: one a Haplortox after 14 years under NT; the other a Hapludox after seven years under NT. Maps (50 x 30 m) of the levels of mechanical penetration resistance (PR) were drawn based on the kriging method, differentiating three compaction degrees (CD): high, intermediate and low. In each CD area, the infiltration rate (initial and steady-state) and cumulative water infiltration were measured using concentric rings, with six replications, and the saturated hydraulic conductivity (K(θs)) was determined using the Guelph permeameter. Statistical evaluation was performed based on a randomized design, using the least significant difference (LSD) test and regression analysis. The steady-state infiltration rate was not influenced by the compaction degree, with mean values of 3 and 0.39 cm h-1 in the Haplortox and the Hapludox, respectively. In the Haplortox, saturated soil hydraulic conductivity was 26.76 cm h-1 at a low CD and 9.18 cm h-1 at a high CD, whereas in the Hapludox, this value was 5.16 cm h-1 and 1.19 cm h-1 for the low and high CD, respectively. The compaction degree did not affect the initial and steady-state water infiltration rate, nor the cumulative water infiltration for either soil type, although the values were higher for the Haplortox than the Hapludox.

Resting-state temporal synchronization networks emerge from connectivity topology and heterogeneity.

Spatial patterns of coherent activity across different brain areas have been identified during the resting-state fluctuations of the brain. However, recent studies indicate that resting-state activity is not stationary, but shows complex temporal dynamics. We were interested in the spatiotemporal dynamics of the phase interactions among resting-state fMRI BOLD signals from human subjects. We found that the global phase synchrony of the BOLD signals evolves on a characteristic ultra-slow (<0.01Hz) time scale, and that its temporal variations reflect the transient formation and dissolution of multiple communities of synchronized brain regions. Synchronized communities reoccurred intermittently in time and across scanning sessions. We found that the synchronization communities relate to previously defined functional networks known to be engaged in sensory-motor or cognitive function, called resting-state networks (RSNs), including the default mode network, the somato-motor network, the visual network, the auditory network, the cognitive control networks, the self-referential network, and combinations of these and other RSNs. We studied the mechanism originating the observed spatiotemporal synchronization dynamics by using a network model of phase oscillators connected through the brain's anatomical connectivity estimated using diffusion imaging human data. The model consistently approximates the temporal and spatial synchronization patterns of the empirical data, and reveals that multiple clusters that transiently synchronize and desynchronize emerge from the complex topology of anatomical connections, provided that oscillators are heterogeneous.

Force volume and stiffness tomography investigation on the dynamics of stiff material under bacterial membranes.

The determination of the characteristics of micro-organisms in clinical specimens is essential for the rapid diagnosis and treatment of infections. A thorough investigation of the nanoscale properties of bacteria can prove to be a fundamental tool. Indeed, in the latest years, the importance of high resolution analysis of the properties of microbial cell surfaces has been increasingly recognized. Among the techniques available to observe at high resolution specific properties of microscopic samples, the Atomic Force Microscope (AFM) is the most widely used instrument capable to perform morphological and mechanical characterizations of living biological systems. Indeed, AFM can routinely study single cells in physiological conditions and can determine their mechanical properties with a nanometric resolution. Such analyses, coupled with high resolution investigation of their morphological properties, are increasingly used to characterize the state of single cells. In this work, we exploit the capabilities and peculiarities of AFM to analyze the mechanical properties of Escherichia coli in order to evidence with a high spatial resolution the mechanical properties of its structure. In particular, we will show that the bacterial membrane is not mechanically uniform, but contains stiffer areas. The force volume investigations presented in this work evidence for the first time the presence and dynamics of such structures. Such information is also coupled with a novel stiffness tomography technique, suggesting the presence of stiffer structures present underneath the membrane layer that could be associated with bacterial nucleoids.