A novel method for automated classification of epileptiform activity in the human electroencephalogram-based on independent component analysis.

Diagnosis of several neurological disorders is based on the detection of typical pathological patterns in the electroencephalogram (EEG). This is a time-consuming task requiring significant training and experience. Automatic detection of these EEG patterns would greatly assist in quantitative analysis and interpretation. We present a method, which allows automatic detection of epileptiform events and discrimination of them from eye blinks, and is based on features derived using a novel application of independent component analysis. The algorithm was trained and cross validated using seven EEGs with epileptiform activity. For epileptiform events with compensation for eyeblinks, the sensitivity was 65 +/- 22% at a specificity of 86 +/- 7% (mean +/- SD). With feature extraction by PCA or classification of raw data, specificity reduced to 76 and 74%, respectively, for the same sensitivity. On exactly the same data, the commercially available software Reveal had a maximum sensitivity of 30% and concurrent specificity of 77%. Our algorithm performed well at detecting epileptiform events in this preliminary test and offers a flexible tool that is intended to be generalized to the simultaneous classification of many waveforms in the EEG.

Ab initio absorption spectrum of NiO combining molecular dynamics with the embedded cluster approach in a discrete reaction field

We developed a procedure that combines three complementary computational methodologies to improve the theoretical description of the electronic structure of nickel oxide. The starting point is a Car-Parrinello molecular dynamics simulation to incorporate vibrorotational degrees of freedom into the material model. By means ofcomplete active space self-consistent field second-order perturbation theory (CASPT2) calculations on embedded clusters extracted from the resulting trajectory, we describe localized spectroscopic phenomena on NiO with an efficient treatment of electron correlation. The inclusion of thermal motion into the theoretical description allowsus to study electronic transitions that, otherwise, would be dipole forbidden in the ideal structure and results in a natural reproduction of the band broadening. Moreover, we improved the embedded cluster model by incorporating self-consistently at the complete active space self-consistent field (CASSCF) level a discrete (or direct) reaction field (DRF) in the cluster surroundings. The DRF approach offers an efficient treatment ofelectric response effects of the crystalline embedding to the electronic transitions localized in the cluster. We offer accurate theoretical estimates of the absorption spectrum and the density of states around the Fermi level of NiO, and a comprehensive explanation of the source of the broadening and the relaxation of the charge transferstates due to the adaptation of the environment

Social meets molecular: Combining phylogenetic and latent class analyses to understand HIV-1 transmission in Switzerland.

Switzerland has a complex human immunodeficiency virus (HIV) epidemic involving several populations. We examined transmission of HIV type 1 (HIV-1) in a national cohort study. Latent class analysis was used to identify socioeconomic and behavioral groups among 6,027 patients enrolled in the Swiss HIV Cohort Study between 2000 and 2011. Phylogenetic analysis of sequence data, available for 4,013 patients, was used to identify transmission clusters. Concordance between sociobehavioral groups and transmission clusters was assessed in correlation and multiple correspondence analyses. A total of 2,696 patients were infected with subtype B, 203 with subtype C, 196 with subtype A, and 733 with recombinant subtypes (mainly CRF02_AG and CRF01_AE). Latent class analysis identified 8 patient groups. Most transmission clusters of subtype B were shared between groups of gay men (groups 1-3) or between the heterosexual groups "heterosexual people of lower socioeconomic position" (group 4) and "injection drug users" (group 8). Clusters linking homosexual and heterosexual groups were associated with "older heterosexual and gay people on welfare" (group 5). "Migrant women in heterosexual partnerships" (group 6) and "heterosexual migrants on welfare" (group 7) shared non-B clusters with groups 4 and 5. Combining approaches from social and molecular epidemiology can provide insights into HIV-1 transmission and inform the design of prevention strategies.

Development of a cost-effective method for platelet-rich plasma (PRP) preparation for topical wound healing.

Platelet-rich plasma (PRP) is a volume of plasma fraction of autologous blood having platelet concentrations above baseline whole-blood values due to processing and concentration. PRP is used in various surgical fields to enhance soft-tissue and bone healing by delivering supra-physiological concentrations of autologous platelets at the site of tissue damage. These preparations may provide a good cellular source of various growth factors and cytokines, and modulate tissue response to injury. Common clinically available materials for blood preparations combined with a two-step centrifugation protocol at 280g each, to ensure cellular component integrity, provided platelet preparations which were concentrated 2-3 fold over total blood values. Costs were shown to be lower than those of other methods which require specific equipment and high-cost disposables, while safety and traceability can be increased. PRP can be used for the treatment of wounds of all types including burns and also of split-thickness skin graft donor sites, which are frequently used in burn management. The procedure can be standardized and is easy to adapt in clinical settings with minimal infrastructure, thus enabling large numbers of patients to benefit from a form of cellular therapy.

Spatial relationship between the atmospheric circulation and the precipitation measured in the western Swiss Alps by means of the analogue method

An adaptation technique based on the synoptic atmospheric circulation to forecast local precipitation, namely the analogue method, has been implemented for the western Swiss Alps. During the calibration procedure, relevance maps were established for the geopotential height data. These maps highlight the locations were the synoptic circulation was found of interest for the precipitation forecasting at two rain gauge stations (Binn and Les Marécottes) that are located both in the alpine Rhône catchment, at a distance of about 100 km from each other. These two stations are sensitive to different atmospheric circulations. We have observed that the most relevant data for the analogue method can be found where specific atmospheric circulation patterns appear concomitantly with heavy precipitation events. Those skilled regions are coherent with the atmospheric flows illustrated, for example, by means of the back trajectories of air masses. Indeed, the circulation recurrently diverges from the climatology during days with strong precipitation on the southern part of the alpine Rhône catchment. We have found that for over 152 days with precipitation amount above 50 mm at the Binn station, only 3 did not show a trajectory of a southerly flow, meaning that such a circulation was present for 98% of the events. Time evolution of the relevance maps confirms that the atmospheric circulation variables have significantly better forecasting skills close to the precipitation period, and that it seems pointless for the analogue method to consider circulation information days before a precipitation event as a primary predictor. Even though the occurrence of some critical circulation patterns leading to heavy precipitation events can be detected by precursors at remote locations and 1 week ahead (Grazzini, 2007; Martius et al., 2008), time extrapolation by the analogue method seems to be rather poor. This would suggest, in accordance with previous studies (Obled et al., 2002; Bontron and Obled, 2005), that time extrapolation should be done by the Global Circulation Model, which can process atmospheric variables that can be used by the adaptation method.

ASIC1a oligomerization structure-function of its extracellular vestibule and functional characterization of αS243Pβɣ ENaC

Quatre cristaux du canal ASIC1a ont été publiés et soutiennent une stoechiométrie trimérique. Cependant, ces données contredisant de précédentes analyses fonctionnelles effectuées sur des canaux de la même famille, notre intérêt fut porté sur l'oligomérisation d'ASIC1a. Dans ce sens, un nouvel essai couplant la méthode d'analyse par substitution de cystéines (SCAM) avec l'utilisation de réactifs sulfhydryls bifonctionnels (crosslinkers) a été mis en place. Le but étant de stabiliser, puis sélectionner les canaux fonctionnels, pour ensuite les séparer selon leur taille par SDS-PAGE. Grâce à cette technique, nous avons démontré que le complexe stabilisé a une taille coïncidant avec une organisation tétramérique. En plus de son oligomérisation, le chemin emprunté par les ions pour traverser le canal n'est pas clairement défini dans ces structures. De ce fait, utilisant une approche électrophysiologique, nous avons étudié le lien entre la structure et la fonction du vestibule extracellulaire d'ASIC1a. Dans ce but, nous nous sommes intéressés l'accessibilité de cystéines spécifiques localisées dans ce vestibule pour des réactifs méthanethiosulfonates (MTS). Ainsi, nous avons pu corréler les cinétiques de modification de ces cystéines par les MTS avec les effets sur le courant sodique, et donc avoir des informations supplémentaires sur la voie empruntée par les ions. De plus, la simulation informatique de liaison de ces réactifs illustre le remplissage total de ce vestibule. Fonctionnellement, cette interaction ne perturbe pas le passage de ions, c'est pourquoi il nous apparaît probable que le vestibule présente une taille plus large que celle illustrée par les cristaux. Dans un deuxième temps, notre intérêt fut porté sur ENaC. Ce canal est composé des trois sous-unités (a, ß et y) et est exprimé dans divers épithéliums, dont les tubules des reins. Il participe à l'homéostasie sodique et est essentiellement régulé par voie hormonale via l'aldostérone et la Vasopressine, mais également par des sérines protéases ou le Na+. Nous avons étudié la répercussion fonctionnelle de la mutation aS243P, découverte chez un nouveau-né prématuré atteint de pseudohypoaldostéronisme de type 1. Cette maladie autosomale récessive se caractérise, généralement, par une hyponatrémie liée à d'importantes pertes de sel dans les urines, une hyperkaliémie, ainsi qu'un niveau élevé d'aldostérone. Tout d'abord aucune des expériences biochimiques et électrophysiologiques n'a pu démontrer un défaut d'expression ou une forte diminution de l'activité soutenant les données cliniques. Cependant, en challengeant aS243PßyENaC avec une forte concentration de Na+ externe, une hypersensibilité de canal fut observée. En effet, ni les phénomènes régulateurs de « feedback inhibition » ou de « Na+ self-inhibition » n'étaient semblables au canal sauvage. De ce fait, ils apparaissaient exacerbés en présence de la mutation, amenant ainsi à une diminution de la réabsorption de Na+. Ceci corrobore entièrement l'hyponatrémie diagnostiquée. Le rein d'un prématuré étant immature, la quantité de Na+ atteignant la partie distale du néphron est plus élevée, du fait que les autres mécanismes de réabsorption en amont ne sont probablement pas encore en place. Cette hypothèse est renforcée par l'existence d'un frère présentant la même mutation, mais qui, né à terme, ne présentait aucun signe d'hyponatrémie. - The main topic of my thesis is the structure-function relationship of the ENaC/Deg family of ion channels, namely the Acid-Sensing Ion Channel ASIC1a and the Epithelial Na Channel ENaC. The primary part of this research is dedicated to the structure of ASIC1a. Four channel crystals have been published, which support a trimeric stoichiometry, although these data contradict previous functional experiments on other ENaC/Deg members. We are therefore interested in ASIC1a oligomerization and have set up a new assay combining the Substituted- Cysteine Accessibility Method (SCAM) with Afunctional sulfhydryl reagents (crosslinkers) allowing its study. The aim was to first stabilize the channels, then select those that are functional and then resolve them according to their size on SDS-PAGE. We demonstrated that the stabilized complex has a molecular weight corresponding to a tetrameric stoichiometry. In addition to our interest in the oligomerization of the ENaC/Deg family of ion channels, we also wanted to investigate the thus far undefined way of permeation for these channels. Therefore, taking the advantage of a more electrophysiological approach, we studied the accessibility of specific cysteines for methanethiosulfonate reagents (MTS) and were able to correlate the MTS association kinetics on cysteine residues with Na+ currents. These results have given us an insight into ion permeation and our functional evidence indicates that the extracellular is larger than that depicted by the crystal structures. As a side project, we focused on ENaC, which is made up of three subunits (a, ß and y) and is expressed in various epithelia, especially in the distal nephron of the kidneys. It plays a role in Na+ homeostasis and is essentially regulated by hormones via aldosterone and vasopressin, but also by serine proteases or Na+. We have studied the functional impact of the aS243P mutation, discovered in a premature baby suffering from pseudohypoaldosteronism of type 1. This autosomal recessive disease is characterized by hyponatremia, hyperkalemia and high aldosterone levels. Firstly, neither biochemical nor electrophysiological experiments indicated an expression defect or a strong decrease in activity. However, challenging aS243PßyENaC with increased external Na+ concentration showed channel hypersensitivity. Indeed, both the "feedback inhibition" and the "Na+ self-inhibition" regulatory mechanisms are impaired, leading to a decrease in Na+ reabsorption, entirely supports the diagnosis. The kidneys in preterm infants are immature and Na+ levels reaching the distal nephron are higher than normally observed. We hypothesize that the upstream reabsorption machinery is unlikely to be sufficiently matured and this assumption is supported by an asymptomatic sibling carrying the same mutation, but born at term. - La cellule, unité fonctionnelle du corps humain, est délimitée par une membrane plasmique servant de barrière biologique entre les milieux intra et extracellulaires. Une communication entre cellules est indispensable pour un fonctionnement adéquat. Sa survie dépend, entre autres, du maintien de la teneur en ions dans chacun des milieux qui doivent pouvoir être réabsorbés, ou sécrétés, selon les besoins. Les protéines insérées dans la membrane forment un canal et sont un moyen de communication permettant spécifiquement à des ions tel que le sodium (Na+) de traverser. Le Na+ se trouve dans la plupart des aliments et le sel, et est spécifiquement réabsorbé au niveau des reins grâce au canal sodique épithélial ENaC. Cette réabsorption se fait de l'urine primaire vers l'intérieur de la cellule, puis est transporté vers le sang. Pour maintenir un équilibre, une régulation de ce canal est nécessaire. En effet, des dysfonctionnements impliquant la régulation ou l'activité d'ENaC lui-même sont à l'origine de maladies telles que la mucoviscidose, l'hypertension ou encore, le pseudohypoaldostéronisme (PHA). Cette maladie est caractérisée, notamment, par d'importantes pertes de sel dans les urines. Des pédiatres ont diagnostiqué un PHA chez un nouveau-né, ce dernier présentant une modification du canal ENaC, nous avons recréé cette protéine afin d'étudier l'impact de ce changement sur son activité. Nous avons démontré que la régulation d'ENaC était effectivement perturbée, conduisant ainsi à une forte réduction de la réabsorption sodique. Afin de développer des molécules capables de moduler l'activité de protéines. Il est nécessaire d'en connaître la structure. Celle du canal sodique sensible à l'acidification ASIC1, un canal cousin d'ENaC, est connue. Ces données structurales contredisant cependant les analyses fonctionnelles, nous nous sommes penchés une nouvelle fois sur ASIC1. Une protéine est une macromolécule biologique composée d'une chaîne d'acides aminés (aa). De l'enchaînement d'aa à la protéine fonctionnelle, quatre niveaux de structuration existent. Chaque aa donne une indication quant au repliement et plus particulièrement la cystéine. Arborant un groupe sulfhydryle (SH) capable de former une liaison spécifique et stable avec un autre SH, celle-ci est souvent impliquée dans la structure tridimensionnelle de la protéine. Ce type de liaison intervient également dans la stabilisation de la structure quaternaire, qui est l'association de plusieurs protéines identiques (homomère), ou pas (hétéromère). Dans cette partie, nous avons remplacé des aa par des cystéines à des endroits spécifiques. Le but était de stabiliser plusieurs homomères d'ASICl ensemble avec des réactifs créant des ponts entre deux SH. Ainsi, nous avons pu déterminer le nombre de protéines ASIC1 participant à la formation d'un canal fonctionnel. Nos résultats corroborent les données fonctionnelles soutenant un canal tétramérique. Nous avons également étudié l'accessibilité de ces nouvelles cystéines afin d'obtenir des informations supplémentaires sur la structure du chemin emprunté par le Na+ à travers ASIC1 et plus particulièrement du vestibule extracellulaire.

Response-to-comments about: "Is it really the method for carbon dioxide determination in human postmortem cardiac gas samples using Headspace-Gas Chromatography-Mass Spectrometry valid?" from T. Saffaj and B. Ihssane

Saffaj et al. recently criticized our method of monitoring carbon dioxide in human postmortem cardiac gas samples using Headspace-Gas Chromatography-Mass Spectrometry. According to the authors, their demonstration, based on the latest SFSTP guidelines (established after 2007 [1,2]) fitted for the validation of drug monitoring bioanalytical methods, has put in evidence potential errors. However, our validation approach was built using SFSTP guidelines established before 2007 [3-6]. We justify the use of these guidelines because of the post-mortem context of the study (and not clinical) and the gaseous state of the sample (and not solid or liquid). Using these guidelines, our validation remains correct.

Combining MAR elements and transposable vectors for more reliable transgene integration and expression

Integration without cytotoxic effects and long-term expression of a transgene constitutes a major challenge in gene therapy and biotechnology applications. In this context, transposons represent an attractive system for gene transfer because of their ability to promote efficient integration of a transgene in a variety of cell lines. However, the transgene integration can lead to insertional mutagenesis and/or unstable transgene expression by epigenetic modifications. These unwanted events may be limited by the use of chromatin control elements called MARs (matrix attachment regions). Indeed, the insertion of these DNA elements next to the transgene usually results in higher and more stable expression by maintaining transgene chromatin in an active configuration and preventing gene silencing. In this study, we tested if the inclusion of the MAR 1-68 in the piggyBac transposon system may lead to efficient and safer transgene integration and ensure reliable stable and long-term expression of a transgene. The MAR-containing transposon construct was tested in CHO cells, for biotechnology applications, and in mesoangioblast cells that can differentiate into muscle cells and are important candidates for potential stem cell therapies of myopathies. We showed that the addition of the MAR 1 -68 in the piggyBac transposon did not interfere with transposition, thereby maintaining high frequency of transgene integrations in these cells. Moreover, the MAR allowed higher transgene expression from fewer transposon integration events. We also found that enriched transgene-expressing cell populations could be obtained without the need of selection pressure. Since antibiotic-enforced selection protocols often result in a higher integrated copy number and mosaic expression patterns, this strategy could benefit many applications in which a low copy number of integrated transgenes and antibiotic-free conditions are desired. In addition, the intramuscular transplantation of mouse tibialis anterior muscles with mesoangioblasts containing the transposon led to widespread and sustained myofiber transgene expression after differentiation of these cells in vivo. These findings indicated that piggyBac vectors may provide a viable approach to achieve stable gene transfer in the context of Duchenne muscular dystrophy therapy. - L'intégration sans effets cytotoxiques et l'expression à long terme d'un transgène constituent un défi majeur en thérapie génique et en biotechnologie. Dans ce contexte, les transposons représentent un système attrayant pour le transfert de gènes en raison de leur capacité à promouvoir l'intégration efficace d'un transgène dans une variété de lignées cellulaires. Toutefois, l'intégration d'un transgène peut conduire à une mutagénèse insertionnelle et/ou à une expression instable due au silençage du transgène suite à des modifications épigénétiques. Ces événements indésirables de silençage génique peuvent être diminués par l'utilisation d'éléments de contrôle de la chromatine appelés MAR (matrix attachment region). En effet, l'insertion de ces éléments d'ADN à proximité du transgène se traduit généralement par une expression plus élevée et plus stable de celui-ci, en permettant le maintien d'une chromatine dans une configuration active autour du transgène et en empêchant l'inactivation du gène. Dans cette étude, nous avons testé si l'inclusion du MAR 1-68 dans le système transposon piggyBac peut améliorer l'efficacité d'intégration de façon sécuritaire et l'expression à long terme d'un transgène. Le transposon contenant l'élément MAR a été testé dans les cellules CHO, couramment utilisées en biotechnologie, et dans des cellules progénitrices appelées mésoangioblastes, qui peuvent se différencier en cellules musculaires, et qui constituent ainsi des candidats prometteurs pour la thérapie à partir de cellules souches de patients souffrant de myopathie. Nous avons montré que l'addition du MAR 1-68 dans le transposon piggyBac n'interfère pas avec la transposition et permet de maintenir une fréquence élevée d'intégration du transgène dans ces deux types cellulaires. De plus, il semble que cette association mène à une meilleure expression du transgène à partir de peu d'événements d'intégration du transposon. En outre, ces populations enrichies en cellules exprimant de façon stable le transgène ont pu être obtenues sans avoir recours à une pression de sélection. Etant donné que les protocoles de sélection basée sur l'utilisation d'antibiotiques conduisent souvent à un nombre plus élevé de copies intégrées et à la variégation de l'expression du transgène et qu'ils impliquent une longue culture in vitro, cette stratégie pourrait profiter à des applications pour lesquelles on souhaite un faible nombre de copies intégrées et/ou l'utilisation d'antibiotiques n'est pas souhaitable. De plus, la transplantation intramusculaire de mésoangioblastes contenant le transposon dans le muscle tibial antérieur de souris a conduit, après la différentiation de ces cellules in vivo, à une expression constante et étendue du transgène dans les myofibres. Ces résultats indiquent que les vecteurs piggyBac pourraient fournir une approche viable pour assurer un transfert de gènes stables dans le contexte d'un traitement de la dystrophic musculaire de Duchenne.

A systematic and feasible method for computing nuclear contributions to electrical properties of polyatomic molecules

An analytic method to evaluate nuclear contributions to electrical properties of polyatomic molecules is presented. Such contributions control changes induced by an electric field on equilibrium geometry (nuclear relaxation contribution) and vibrational motion (vibrational contribution) of a molecular system. Expressions to compute the nuclear contributions have been derived from a power series expansion of the potential energy. These contributions to the electrical properties are given in terms of energy derivatives with respect to normal coordinates, electric field intensity or both. Only one calculation of such derivatives at the field-free equilibrium geometry is required. To show the useful efficiency of the analytical evaluation of electrical properties (the so-called AEEP method), results for calculations on water and pyridine at the SCF/TZ2P and the MP2/TZ2P levels of theory are reported. The results obtained are compared with previous theoretical calculations and with experimental values

BrainCheck - a very brief tool to detect incipient cognitive decline: optimized case-finding combining patient- and informant-based data.

INTRODUCTION: Optimal identification of subtle cognitive impairment in the primary care setting requires a very brief tool combining (a) patients' subjective impairments, (b) cognitive testing, and (c) information from informants. The present study developed a new, very quick and easily administered case-finding tool combining these assessments ('BrainCheck') and tested the feasibility and validity of this instrument in two independent studies. METHODS: We developed a case-finding tool comprised of patient-directed (a) questions about memory and depression and (b) clock drawing, and (c) the informant-directed 7-item version of the Informant Questionnaire on Cognitive Decline in the Elderly (IQCODE). Feasibility study: 52 general practitioners rated the feasibility and acceptance of the patient-directed tool. Validation study: An independent group of 288 Memory Clinic patients (mean ± SD age = 76.6 ± 7.9, education = 12.0 ± 2.6; 53.8% female) with diagnoses of mild cognitive impairment (n = 80), probable Alzheimer's disease (n = 185), or major depression (n = 23) and 126 demographically matched, cognitively healthy volunteer participants (age = 75.2 ± 8.8, education = 12.5 ± 2.7; 40% female) partook. All patient and healthy control participants were administered the patient-directed tool, and informants of 113 patient and 70 healthy control participants completed the very short IQCODE. RESULTS: Feasibility study: General practitioners rated the patient-directed tool as highly feasible and acceptable. Validation study: A Classification and Regression Tree analysis generated an algorithm to categorize patient-directed data which resulted in a correct classification rate (CCR) of 81.2% (sensitivity = 83.0%, specificity = 79.4%). Critically, the CCR of the combined patient- and informant-directed instruments (BrainCheck) reached nearly 90% (that is 89.4%; sensitivity = 97.4%, specificity = 81.6%). CONCLUSION: A new and very brief instrument for general practitioners, 'BrainCheck', combined three sources of information deemed critical for effective case-finding (that is, patients' subject impairments, cognitive testing, informant information) and resulted in a nearly 90% CCR. Thus, it provides a very efficient and valid tool to aid general practitioners in deciding whether patients with suspected cognitive impairments should be further evaluated or not ('watchful waiting').

Adaptação para a língua portuguesa e validação do Lunney Scoring Method for Rating Accuracy of Nursing Diagnoses

O Lunney Scoring Method for Rating Accuracy of Nursing Diagnoses (LSM) é uma escala de diferencial semântico que foi desenvolvida por Lunney para estimar a acurácia dos diagnósticos de enfermagem. O objetivo deste estudo foi adaptar o LSM para a língua portuguesa e avaliar as sua propriedades psicométricas. A escala original foi traduzida para o português, revertida para o inglês e as duas versões em inglês foram comparadas para ajustar a versão em português que passou a ser denominada Escala de Acurácia de Diagnóstico de Enfermagem de Lunney - EADE. Quatro enfermeiras foram orientadas sobre a EADE e a aplicaram em 159 diagnósticos formulados para 26 pacientes de três estudos primários com base nos registros de entrevista e exame físico de cada paciente. Os índices Kappa de Cohen mostraram ausência de concordância entre as avaliadoras, o que indica que o instrumento adaptado não tem confiabilidade satisfatória. Em virtude desse resultado, não foi realizada estimativa de validade.

Select-divide-and-conquer method for large-scale configuration interaction

A select-divide-and-conquer variational method to approximate configuration interaction (CI) is presented. Given an orthonormal set made up of occupied orbitals (Hartree-Fock or similar) and suitable correlation orbitals (natural or localized orbitals), a large N-electron target space S is split into subspaces S0,S1,S2,...,SR. S0, of dimension d0, contains all configurations K with attributes (energy contributions, etc.) above thresholds T0={T0egy, T0etc.}; the CI coefficients in S0 remain always free to vary. S1 accommodates KS with attributes above T1≤T0. An eigenproblem of dimension d0+d1 for S0+S 1 is solved first, after which the last d1 rows and columns are contracted into a single row and column, thus freezing the last d1 CI coefficients hereinafter. The process is repeated with successive Sj(j≥2) chosen so that corresponding CI matrices fit random access memory (RAM). Davidson's eigensolver is used R times. The final energy eigenvalue (lowest or excited one) is always above the corresponding exact eigenvalue in S. Threshold values {Tj;j=0, 1, 2,...,R} regulate accuracy; for large-dimensional S, high accuracy requires S 0+S1 to be solved outside RAM. From there on, however, usually a few Davidson iterations in RAM are needed for each step, so that Hamiltonian matrix-element evaluation becomes rate determining. One μhartree accuracy is achieved for an eigenproblem of order 24 × 106, involving 1.2 × 1012 nonzero matrix elements, and 8.4×109 Slater determinants

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