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

A new efficient approach to the direct restricted active space self-consistent field method

An implicitly parallel method for integral-block driven restricted active space self-consistent field (RASSCF) algorithms is presented. The approach is based on a model space representation of the RAS active orbitals with an efficient expansion of the model subspaces. The applicability of the method is demonstrated with a RASSCF investigation of the first two excited states of indole

Variational calculation of static and dynamic vibrational nonlinear optical properties

The vibrational configuration interaction method used to obtain static vibrational (hyper)polarizabilities is extended to dynamic nonlinear optical properties in the infinite optical frequency approximation. Illustrative calculations are carried out on H2 O and N H3. The former molecule is weakly anharmonic while the latter contains a strongly anharmonic umbrella mode. The effect on vibrational (hyper)polarizabilities due to various truncations of the potential energy and property surfaces involved in the calculation are examined

Calculation of Franck-Condon factors including anharmonicity: simulation of the C2 H4+ X̃2B 3u←C2 H4X̃1 Ag band in the photoelectron spectrum of ethylene

Our new simple method for calculating accurate Franck-Condon factors including nondiagonal (i.e., mode-mode) anharmonic coupling is used to simulate the C2H4+X2B 3u←C2H4X̃1 Ag band in the photoelectron spectrum. An improved vibrational basis set truncation algorithm, which permits very efficient computations, is employed. Because the torsional mode is highly anharmonic it is separated from the other modes and treated exactly. All other modes are treated through the second-order perturbation theory. The perturbation-theory corrections are significant and lead to a good agreement with experiment, although the separability assumption for torsion causes the C2 D4 results to be not as good as those for C2 H4. A variational formulation to overcome this circumstance, and deal with large anharmonicities in general, is suggested

A comparative analysis of two methods for the calculation of electric-field-induced perturbations to molecular vibration

Two common methods of accounting for electric-field-induced perturbations to molecular vibration are analyzed and compared. The first method is based on a perturbation-theoretic treatment and the second on a finite-field treatment. The relationship between the two, which is not immediately apparent, is made by developing an algebraic formalism for the latter. Some of the higher-order terms in this development are documented here for the first time. As well as considering vibrational dipole polarizabilities and hyperpolarizabilities, we also make mention of the vibrational Stark effec

The curvature of the conical intersection seam: an approximate second-order analysis

We present a method for analyzing the curvature (second derivatives) of the conical intersection hyperline at an optimized critical point. Our method uses the projected Hessians of the degenerate states after elimination of the two branching space coordinates, and is equivalent to a frequency calculation on a single Born-Oppenheimer potential-energy surface. Based on the projected Hessians, we develop an equation for the energy as a function of a set of curvilinear coordinates where the degeneracy is preserved to second order (i.e., the conical intersection hyperline). The curvature of the potential-energy surface in these coordinates is the curvature of the conical intersection hyperline itself, and thus determines whether one has a minimum or saddle point on the hyperline. The equation used to classify optimized conical intersection points depends in a simple way on the first- and second-order degeneracy splittings calculated at these points. As an example, for fulvene, we show that the two optimized conical intersection points of C2v symmetry are saddle points on the intersection hyperline. Accordingly, there are further intersection points of lower energy, and one of C2 symmetry - presented here for the first time - is found to be the global minimum in the intersection space

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

One- and two-center physical space partitioning of the energy in the density functional theory

A conceptually new approach is introduced for the decomposition of the molecular energy calculated at the density functional theory level of theory into sum of one- and two-atomic energy components, and is realized in the "fuzzy atoms" framework. (Fuzzy atoms mean that the three-dimensional physical space is divided into atomic regions having no sharp boundaries but exhibiting a continuous transition from one to another.) The new scheme uses the new concept of "bond order density" to calculate the diatomic exchange energy components and gives them unexpectedly close to the values calculated by the exact (Hartree-Fock) exchange for the same Kohn-Sham orbitals

Proximal tibia volumetric bone mineral density is correlated to the magnitude of local acceleration in male long distance runners.

Introduction: The beneficial effect of physical exercise on bone mineral density (BMD) is at least partly explained by the forces exerted directly on the bones. Male runners present generally higher BMD than sedentary individuals. We postulated that the proximal tibia BMD is related to the running distance as well as to the magnitude of the shocks (while running) in male runners. Methods: A prospective study (three yearly measurements) included 81 healthy male subjects: 16 sedentary lean subjects and three groups of runners (5-30 km/week, n=19; 30-50 km/week, n=29; 50-100 km/week, n=17). Several measurements were performed at the proximal tibia level: volumetric BMD (vBMD), cortical index (CI) i.e. an index of cortical bone thickness and peak accelerations (an index of shocks during heel strike) while running (measured by a 3-D accelerometer). A general linear model assessed the prediction of vBMD or CI by a) simple effects (running distance, peak accelerations, time) and b) interactions (for instance if vBMD prediction by peak acceleration depends on running distance). Results: CI and vBMD a) increase with running distance to reach a plateau over 30 km/wk, b) are positively associated with peak accelerations over 30 km/week. Discussion: Running may be associated with high peak accelerations in order to have beneficial effects on BMD. More important strains are needed to be associated with the same increase in BMD during running sessions of short duration than those of long duration. Conclusion: CI and vBMD are associated with the magnitude of the shocks during heel strike in runners. Key words: Bone mineral density, strains, physical exercise, running distance.

Low dispersing species : population genetic processes in space and time : the genus "Trochulus" (Gastropoda: Hygromiidae) as a case study

RÉSUMÉ Une espèce est rarement composée d'une population unique. Parce que les individus ont des capacités de dispersion limitées et que les paysages sont des mosaïques d'habitats, la plupart des espèces sont plutôt composées de sous-populations connectées par la migration. Cette variation spatiale influence directement la distribution de la variabilité génétique dans et entre les populations. Durant ce travail, nous avons abordé certains des processus populationnels qui ont joué un rôle supposé dans l'apparition de nouvelles espèces au sein du genre Trochulus. Plus précisément, nous avons tenté d'évaluer les impacts respectifs de l'isolement passé (facteurs historiques) et présent (facteurs locaux). Nous avons d'abord pu montrer que les faibles capacités de dispersion des escargots terrestres ont directement influencé leur histoire évolutive à toutes les échelles spatiales et temporelles. En réduisant l'effet homogénéisant de la migration, une faible dispersion maintient dans les populations les traces génétiques d'évènements passés. A l'échelle de la distribution globale de Trochulus villosus, ces traces ont permis de reconstruire une histoire faite d'isolements et d'expansions de populations. En combinant des données génétiques avec une modélisation de la niche climatique passée, il a été possible de proposer un scénario significativement meilleur que toutes les hypothèses alternatives que nous avons testées. A l'échelle locale par contre, l'héritage historique est difficile à distinguer de la dynamique actuelle. Ce fut le cas des lignées mitochondriales du complexe sericeus-hispidus : les deux principales lignées étaient phylogénétiquement éloignées, avaient eu des démographies passées différentes et corrélaient avec des différences morphologiques. D'un autre côté, le flux de gène nucléaire était fort, contredisant l'idée de deux espèces cryptiques isolées reproductivement. Pour pouvoir conclure à la présence ou non de deux espèces, il nous a manqué des informations locales sur la dynamique des populations et les conditions écologiques que l'on trouve dans la région d'étude. Enfin, nous avons pu souligner que la connectivité entre populations d'escargots est soumise à la qualité des habitats et à leur organisation spatiale. Les escargots sont dépendants d'un habitat et s'y adaptent, comme l'indiquent la présence de «poils » uniquement sur la coquille d'espèces vivant dans des habitats humides ou la corrélation entre morphologie et habitat au sein du complexe sericeus-hispidus. Logiquement donc, les escargots migrent préférentiellement au travers d'habitats favorables comme l'a montré la réduction de flux de gènes au travers des prairies chez T. villosus (une espèce forestière). De ces données, nous pouvons supposer que les populations d'escargots en particulier, et des espèces à faible dispersion en général, ont de fortes chances d'être affectées par les changements climatiques, avec de probables implications pour leurs histoires évolutives. SUMMARY : Species rarely consists in a single population. Because individuals have limited dispersal abilities, because landscapes are habitat patchworks, most species are made of several subpopulations connected by migration. This spatial variation has consequences on the distribution of genetic diversity within and between populations, creating a structure among the populations. During the present work, we investigated some of the population processes assumed to have played an important role on the speciation within the genus Trochulus. More specifically, we questioned the respective impacts of past (historical factors) or present (local factors) population isolations. We first could show that the poor dispersal abilities of land snails have had profound impacts on their evolutionary histories at all spatial and temporal scales. Low dispersal maintains a strong signature of past events in the populations by minimising the homogenising effects of geneflow. At the scale of Trochulus villosus global distribution, they allowed to retrieve the detailed history of this species population isolations and expansions. Combining a large genetic dataset with paleo-climatic niche modelling ended up with a historical scenario significantly better than all traditional alternatives we tested. At local scale on the contrary, past events become difficult to tease apart from ongoing processes. This was the case for the divergent mitochondria) lineages within the sericeus-hispidus complex: the two principal lineages appeared to be phylogenetically distant, to have experienced different demographic histories and to correlate with morphological differences. On the other hand, nuclear (present day) geneflow was high, contradicting the idea of two reproductively isolated cryptic species. Information on the local population dynamics and environmental conditions are lacking to be able to decide whether past isolation has indeed resulted here in new species. Finally, we emphasised the importance of the habitat types present in a landscape as well as their spatial organisation for the population connectivity of land snails. These species are tightly dependent on a habitat and adapt to it as shown by thé occurrence of hair-like structures only in species living in humid environments or by the correlation between shell morphology and habitat in the sericeus-hispidus complex. As a result, land snails preferentially migrate through favourable habitats: Trochulus villosus, a forest species, had its geneflow significantly reduced across meadows. From these data, we can hypothesise that the populations of land snails in particular and of low dispersing species in general are likely to be strongly affected by the ongoing climate changes, with potential major consequences on their evolutionary histories.

Long-distance wound signalling in arabidopsis

The leaves of all plants use elaborate and inducible defence systems to protect themselves. A wide variety of such defences are known and they include defence chemicals such as alkaloids, phenolics and terpenes, physical structures ranging from fibre cells to silica deposits, and a wide variety of defence proteins many of which target digestive processes in herbivores. It has long been known that the defence responses of plants under attack by insects are not restricted to the site of attack. Instead, if a leaf is damaged, defence can be triggered in other parts of the plant body, for example in distal leaves or even in roots and flowers. This raises the question of what are the organ-to-organ signals that coordinate this process. Several hypotheses have been proposed. These include the long-distance transfer of chemical signals through the plant vasculature, hydraulic signals that may transit through the xylem, and electrical signals that would move through living tissues such as the phloem. Much evidence for each of these scenarios has been published. In this thesis we took advantage of the fact that many plant defence responses are regulated by a signal transduction pathway based on a molecule called jasmonic acid. We used this molecule, one of its derivatives (jasmonoyl-isoleucine), and some of the genes it regulates as markers. Using these we investigated the possible role of the electrical signals in the leaf- to-leaf activation of the jasmonate pathway. We found that feeding insects stimulate easily detected electrical activity in the leaves of Arabidopsis thaliana and we used non-invasive surface electrodes to record this activity. This approach showed that jasmonate pathway activity and the electrical activity provoked by mechanical wounding occurred within identical spatial boundaries. Measurements of the apparent speed of surface potentials agreed well with previous velocity estimates for the speed of leaf-to-leaf signals that activate the jasmonate pathway. Using this knowledge we were able to investigate the effects of current injection into Arabidopsis leaves. This resulted in the strong expression of many jasmonate-regulated genes. All these results showed that electrical activity and the activation of jasmonate signalling were highly correlated. In order to test for possible causal links between the two processes, we conducted a small-scale reverse genetic screen on a series of T-DNA insertion mutants in ion channel genes and in other genes encoding proteins such as proton pumps. This screen, which was based on surface potential measurements, revealed that mutations in genes related to ionotropic glutamate receptors in animals had impaired electrical activity after wounding. Combining mutation of two of these glutamate-receptor-like genes in a double mutant reduced the response of leaves to current injection. When a leaf of this double mutant was wounded it failed to transmit a long-distance signal to a distal leaf. This result distinguished the double mutant from the wild-type plant and provides the first genetic evidence that electrical signalling is necessary to coordinate defence responses between organs in plants. - Les feuilles des plantes disposent de systèmes de défense inductibles très élaborés. Un grand nombre de ces systèmes de défenses sont connus et sont basés sur des composés chimiques comme les alcaloïdes, les composés phénoliques ou les terpènes, des systèmes physiques allant de la production de cellules fibreuses aux cristaux de silice ainsi qu'un grand nombre de protéines de défense ciblant le processus digestif des herbivores. Il est connu dépuis longtemps que la réponse défensive de la plante face à l'attaque pas un insecte n'est pas seulement localisée au niveau de la zone d'attaque. A la place, si une feuille est attaquée, les systèmes de défense peuvent être activés ailleurs dans la plante, comme par exemple dans d'autres feuilles, les racines ou même les fleurs. Ces observations soulèvent la question de la nature des signaux d'organes à organes qui régulent ces systèmes. Plusieurs hypothèses ont été formulées; une ou plusieures molécules pourraient être véhiculées dans la plante grâce au système vasculaire, un signal hydraulique transmis au travers du xylème ou encore des signaux électriques transmis par les cellules comme dans le phloème par exemple. De nombreuses études ont été publiées sur ces différentes hypothèses. Dans ce travail de thèse, nous avons choisi d'utiliser à notre avantage le fait que de nombreuses réponses de défense de la plante sont régulées par une même voie de signalisation utilisant l'acide jasmonique. Nous avons utilisé comme marqueurs cette molécule, un de ses dérivés (le jasmonoyl-isoleucine) ainsi que certains des gènes que l'acide jasmonique régule. Nous avons alors testé l'implication de la transmission de signaux électriques dans l'activation de la voie du jasmonate de feuille à feuille. Nous avons découvert que les insectes qui se nourrissent de feuilles d'Arabidopsis thaliana activent un signal électrique que nous avons pu mesurer grâce à une technique non invasive d'électrodes de surface. Les enregistrements ont montré que la génération de signaux électriques et l'activation de la voie du jasmonate avaient lieu aux mêmes endroits. La mesure de la vitesse de déplacement des impulsions électriques correspond aux estimations faites concernant l'activation de la voie du jasmonate. Grâce à cela, nous avons pu tester l'effet d'injection de courant électrique dans les feuilles d'Arabidopsis. La conséquence a été une forte expression de nombreux gènes de la voie du jasmonate, suggérant une forte corrélation entre l'activité électrique et l'activation de la voie du jasmonate. Afin de tester le lien de cause entre ces deux phénomènes, nous avons entrepris un criblage génétique sur une série de mutants d'insertion à l'ADN-T dans des gènes de canaux ioniques et d'autres gènes d'intérêt comme les gènes des pompes à protons. Ce criblage, basé sur la mesure de potentiels de surface, a permis de montrer que plusieurs mutations de gènes liés aux récepteurs au glutamate ionotropique présentent une baisse drastique de leurs activités électriques après une blessure mécanique des feuilles par rapport au type sauvage. Par la combinaison de deux mutations de ces récepteurs au glutamate en un double mutant, on obtient une réponse à la stimulation électrique encore plus faible. Quand une feuille du double mutant est blessée, elle est incapable de transmettre un signal à longue distance vers une feuille éloignée. Ce résultat permet de distinguer le double mutant de la plante sauvage et amène la première preuve génétique que l'activité électrique est nécessaire pour coordonner les réponses de défense entre les organes chez les plantes.

Improved earthquake response via simulation and integrated space- and ground-based technologies: the TREMOR proposal

Earthquakes occurring around the world each year cause thousands ofdeaths, millions of dollars in damage to infrastructure, and incalculablehuman suffering. In recent years, satellite technology has been asignificant boon to response efforts following an earthquake and itsafter-effects by providing mobile communications between response teamsand remote sensing of damaged areas to disaster management organizations.In 2007, an international team of students and professionals assembledduring theInternational Space University’s Summer Session Program in Beijing, Chinato examine how satellite and ground-based technology could be betterintegrated to provide an optimised response in the event of an earthquake.The resulting Technology Resources for Earthquake MOnitoring and Response(TREMOR) proposal describes an integrative prototype response system thatwill implement mobile satellite communication hubs providing telephone anddata links between response teams, onsite telemedicine consultation foremergency first-responders, and satellite navigation systems that willlocate and track emergency vehicles and guide search-and-rescue crews. Aprototype earthquake simulation system is also proposed, integratinghistorical data, earthquake precursor data, and local geomatics andinfrastructure information to predict the damage that could occur in theevent of an earthquake. The backbone of these proposals is a comprehensiveeducation and training program to help individuals, communities andgovernments prepare in advance. The TREMOR team recommends thecoordination of these efforts through a centralised, non-governmentalorganization.

A Review of the State's Leased Space and the Leasing Process, General Services Enterprise of the Department of Administrative Services, July 1, 1998 through June 30, 2003

Other Audit Reports - State Leasing