An explorative study of inflammation-related proteins associated with kidney injury in male heat-stressed workers.

Chronic kidney disease of non-traditional origin (CKDnt) and acute kidney injury (AKI) often affect heat-stressed Mesoamerican manual workers. Inflammation occurs concurrently with AKI in this population, but its role remains unknown. To explore links between inflammation and kidney injury in heat stress, we compared levels of inflammation-related proteins in cutters with and without increasing serum creatinine levels during sugarcane harvest. These sugarcane cutters have previously been identified to be repeatedly exposed to severe heat stress during the five month harvest season. A nested case-control study was conducted among male Nicaraguan sugarcane cutters in a CKDnt hotspot. Cases (n = 30) were defined as having an increase in creatinine of ≥0.3 mg/dL across the five-month harvest. Controls (n = 57) had stable creatinine levels. Ninety-two inflammation-related proteins in serum were measured before and after harvest using Proximity Extension Assays. Mixed linear regression was used to identify differences in protein concentrations between cases and controls before harvest, differential trends during harvest, and association between protein concentrations and the urine kidney injury markers Kidney Injury Molecule (KIM)-1, Monocyte Chemoattractant Protein (MCP)-1 and albumin. One protein, chemokine (C-C motif) ligand 23 (CCL23), was elevated among cases at pre-harvest. Changes in seven inflammation-related proteins (CCL19, CCL23, colony-stimulating factor 1 [CSF1], hepatocyte and fibroblast growth factors [HGF and FGF23], and tumor necrosis factor beta [TNFB] and TNF-related activation-induced cytokine [TRANCE]) were associated with case status and at least two out of three urine kidney injury markers (KIM-1, MCP-1 and albumin). Several of these have been implicated in myofibroblast activation, which likely is an important step in kidney interstitial fibrotic disease such as CKDnt. This study provides an initial exploration of immune system determinants of, and activation during, kidney injury experienced during prolonged heat stress.

Large difference but high correlation between creatinine and cystatin C estimated glomerular filtration rate in Mesoamerican sugarcane cutters.

OBJECTIVES: To explore the relationship between creatinine and cystatin C based estimated glomerular filtration rate (eGFR) in actively working sugarcane cutters. METHODS: This cohort study included 458 sugarcane cutters from Nicaragua and El Salvador. Serum samples were taken before and at end of harvest seasons and analysed for creatinine and cystatin C. Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) formulas were used to calculate eGFRs based on creatinine (eGFRcr), cystatin C (eGFRcys) and both creatinine and cystatin C (eGFRcrcys) at each time point. Bland-Altman plots and paired t-tests were used to compare the difference between eGFRcr and eGFRcys, and the difference in eGFRs between before and at end of the harvest seasons. RESULTS: The mean eGFRcr was higher than eGFRcys in both cohorts; absolute difference 22 mL/min/1.73 m2 (95% CI 21 to 23) in Nicaragua and 13 mL/min/1.73 m2 (95% CI 11 to 15) in El Salvador. Correlations between eGFRcr and eGFRcys were high, with r=0.69, 0.77 and 0.67 in Nicaragua at pre-harvest, end-harvest and cross-harvest, and r=0.89, 0.89 and 0.49 in El Salvador. CONCLUSIONS: Creatinine increases among heat-stressed workers reflect reduced glomerular filtration as estimated using eGFRcys, a marker independent of muscle mass and metabolism. The discrepancy between eGFRcr and eGFRcys may indicate reduced glomerular filtration of larger molecules and/or systemic bias in CKD-EPI performance in this population.

Effectiveness and Safety of the European Society of Cardiology 0-/1-h Troponin Rule-Out Protocol: The Design of the ESC-TROP Multicenter Implementation Study.

Chest pain is one of the most common complaints at the emergency department (ED), and it is commonly the perceived likelihood of acute coronary syndrome (ACS) that drives management. Guidelines from the European Society of Cardiology (ESC) recommend the use of a 0-/1-h high-sensitivity cardiac troponin T (hs-cTnT) protocol to rule out or in ACS, but this is mostly based on observational studies. The aim of the ESC-TROP trial is to determine the safety and effectiveness of the ESC 0-/1-h hs-cTnT protocol when implemented in routine care. Adult chest pain patients at 5 EDs in the Skåne Region, Sweden, are included in the trial. The 0-/1-h hs-cTnT ESC protocol supplemented with clinical assessment and electrocardiography (ECG) is implemented at 3 EDs, and the other 2 EDs act as concurrent controls. Outcomes will be evaluated during the 10 months after the implementation and the corresponding 10 months of the previous year. The 2 co-primary outcomes are (a) acute myocardial infarction (AMI) and all-cause death within 30 days in patients discharged from the ED, and (b) ED length of stay of the same patients. Secondary outcomes include the proportion of chest pain patients discharged from the ED and the number of ruled-out patients undergoing objective testing within 30 days. The ESC-TROP trial will determine the performance and applicability of the 0-/1-h hs-cTnT ESC protocol supplemented with clinical assessment and ECG when implemented in routine ED care. It will provide evidence whether 0-/1-h hs-cTnT testing is safe, effective, and feasible, and whether widespread implementation as recommended by ESC guidelines should be supported.

Preventing kidney injury among sugarcane workers: promising evidence from enhanced workplace interventions.

OBJECTIVES: To assess if improvement of working conditions related to heat stress was associated with improved kidney health outcomes among sugarcane harvest workers in Chichigalpa, Nicaragua, a region heavily affected by the epidemic of chronic kidney disease of non-traditional origin. METHODS: Based on our findings during the 2017-2018 harvest (harvest 1), recommendations that enhanced the rest schedule and improved access to hydration and shade were given before the 2018-2019 harvest (harvest 2). Actual work conditions during harvest 2 were then observed. Serum creatinine (SCr) was measured before and at end-harvest, and cross-harvest changes in estimated glomerular filtration rate (eGFR) and incident kidney injury (IKI, ie, SCr increase by ≥0.30 mg/dL or ≥1.5 times the baseline value) were compared between harvest 1 and harvest 2 for three jobs with different physical workloads using regression modelling. Workers who left during harvest were contacted at home, to address the healthy worker selection effect. RESULTS: In burned cane cutters, mean cross-harvest eGFR decreased 6 mL/min/1.73 m2 (95% CI 2 to 9 mL/min/1.73 m2) less and IKI was 70% (95% CI 90% to 50%) lower in harvest 2 as compared with harvest 1 data. No such improvements were seen among seed cutters groups with less successful intervention implementation. CONCLUSION: Kidney injury risk was again elevated in workers with strenuous jobs. The results support further efforts to prevent kidney injury among sugarcane workers, and other heat-stressed workers, by improving access to water, rest and shade. The distinction between design and implementation of such interventions should be recognised.

Pregnancy-induced changes in serum concentrations of perfluoroalkyl substances and the influence of kidney function.

BACKGROUND: Epidemiological associations between maternal concentrations of perfluoroalkyl substances (PFAS) and birth weight are inconsistent. There is concern that studies based on samples collected in late pregnancy may be confounded by kidney function but studies of the relation between pregnancy-induced changes in PFAS and kidney function are lacking. Our aims were to investigate changes in serum concentrations of perfluorononanoic acid (PFNA), perfluorooctanoic acid (PFOA), perfluorooctane sulfonate (PFOS) and perfluorohexane sulfonate (PFHxS) from early to late pregnancy and to explore relations to changes in glomerular filtration rate (GFR) and glomerular pore size. METHODS: We conducted the study in a cohort of 73 pregnancies of normal-weight Swedish women without gestational diabetes and preeclampsia, enrolled 2009-2014. Blood was collected in median weeks 11 and 36, respectively, and analysed PFAS using liquid chromatography-tandem-mass-spectrometry. We estimated GFR based on creatinine and cystatin C and used the ratio eGFRcystatin C/eGFRcreatinine to indicate glomerular pore size. We used Wilcoxon signed-rank test to compare early and late measures and partial Spearman rank correlations to explore relations between changes in PFAS and kidney function. RESULTS: Median concentrations of PFNA, PFOA and PFOS decreased by 15-21% but changes were uncorrelated to changes in kidney function (partial R = - 0.06-0.11). The observed increase in median PFHxS concentration of 69% was likely an artefact of systematic measurement error caused by coeluting endogenous inferences. CONCLUSIONS: Serum concentrations of PFNA, PFOA and PFOS decrease during pregnancy but the magnitudes of change are unrelated to parallel changes in eGFR and glomerular pore size, suggesting that changes in these indicators of kidney function are not important confounders in studies of PFAS and birth weight in pregnancies without gestational diabetes and preeclampsia.

ReSpect: Relativistic spectroscopy DFT program package.

With the increasing interest in compounds containing heavier elements, the experimental and theoretical community requires computationally efficient approaches capable of simultaneous non-perturbative treatment of relativistic, spin-polarization, and electron correlation effects. The ReSpect program has been designed with this goal in mind and developed to perform relativistic density functional theory (DFT) calculations on molecules and solids at the quasirelativistic two-component (X2C Hamiltonian) and fully relativistic four-component (Dirac-Coulomb Hamiltonian) level of theory, including the effects of spin polarization in open-shell systems at the Kramers-unrestricted self-consistent field level. Through efficient algorithms exploiting time-reversal symmetry, biquaternion algebra, and the locality of atom-centered Gaussian-type orbitals, a significant reduction of the methodological complexity and computational cost has been achieved. This article summarizes the essential theoretical and technical advances made in the program, supplemented by example calculations. ReSpect allows molecules with >100 atoms to be efficiently handled at the four-component level of theory on standard central processing unit-based commodity clusters, at computational costs that rarely exceed a factor of 10 when compared to the non-relativistic realm. In addition to the prediction of band structures in solids, ReSpect offers a growing list of molecular spectroscopic parameters that range from electron paramagnetic resonance parameters (g-tensor, A-tensor, and zero-field splitting), via (p)NMR chemical shifts and nuclear spin-spin couplings, to various linear response properties using either conventional or damped-response time-dependent DFT (TDDFT): excitation energies, frequency-dependent polarizabilities, and natural chiroptical properties (electronic circular dichroism and optical rotatory dispersion). In addition, relativistic real-time TDDFT electron dynamics is another unique feature of the program. Documentation, including user manuals and tutorials, is available at the program's website http://www.respectprogram.org.

Insufficient mixing of thawed serum samples leading to erroneous results - experience from a field study and use of a correction procedure.

Incorrect analysis results that are close to expected might not be recognized in scientific studies or routine patient care. In two field studies we obtained unexpected results in a large number of samples. The present study aimed to identify the source of error in the samples from these studies and to validate a method to obtain correct results. Pre-analytical procedures were scrutinized, giving no indications of inappropriate pre-analytical sample handling in the field or during transport in a tropical climate. Using a new set of samples from volunteers in simulation experiments, we observed the known concentration gradient of analytes sampled in gel as well as plain tubes after freezer storage and thawing. Experiments demonstrated that mixing of samples by vortexing alone was not sufficient to disrupt the gradient formed by freezing and thawing, which appeared to cause the problem encountered when we in field studies analyzed and biobanked large sample sets by robot pipetting. A correction procedure was introduced, in which the obtained value of an analyte was multiplied by a correction factor calculated for each sample using the expected sodium level (140 mmol/L) divided by the measured sodium value. When it was validated on results from the simulation experiments, we repeatedly found that the correction lead to results very close to true values for analytes of different size and charge. Usefulness of the procedure was demonstrated when applied to a large set of field study results.

Workload and cross-harvest kidney injury in a Nicaraguan sugarcane worker cohort.

OBJECTIVES: To examine the association between workload and kidney injury in a fieldworker cohort with different levels of physically demanding work over a sugarcane harvest, and to assess whether the existing heat prevention efforts at a leading occupational safety and health programme are sufficient to mitigate kidney injury. METHODS: Biological and questionnaire data were collected before (n=545) and at the end (n=427) of harvest among field support staff (low workload), drip irrigation workers (moderate), seed cutters (high) and burned sugarcane cutters (very high). Dropouts were contacted (87%) and reported the reason for leaving work. Cross-harvest incident kidney injury (IKI) was defined as serum creatinine increase ≥0.30 mg/dL or ≥1.5 times the baseline value, or among dropouts reporting kidney injury leading to leaving work. RESULTS: Mean cross-harvest estimated glomerular filtration rate change was significantly associated with workload, increasing from 0 mL/min/1.73 m2 in the low-moderate category to -5 mL/min/1.73 m2 in the high and -9 mL/min/1.73 m2 in the very high workload group. A similar pattern occurred with IKI, where low-moderate workload had 2% compared with 27% in the very high workload category. A healthy worker selection effect was detected, with 32% of dropouts reporting kidney injury. Fever and C reactive protein elevation were associated with kidney injury. CONCLUSIONS: Workers considered to have the highest workload had more cross-harvest kidney damage than workers with less workload. Work practices preventing heat stress should be strengthened and their role in preventing kidney damage examined further. Future occupational studies on chronic kidney disease of unknown aetiology should account for a healthy worker effect by pursuing those lost to follow-up.

Low-level internalization of cystatin E/M affects legumain activity and migration of melanoma cells.

The ratio between proteases and their inhibitors is unbalanced in cancer. The cysteine protease inhibitor cystatin C is internalized by some cancer cells, which affects cellular properties. Here we aimed to investigate if uptake of cystatin C and the related inhibitor cystatin E/M occur in melanoma cell lines and to evaluate to what extent the uptake affects the legumain activity that is typically increased in melanoma. First we studied the basic expression, secretion, and intracellular content of all type 2 cystatins as well as expression and activity of their possible target enzymes legumain and cathepsin B in MDA-MB-435S, A375, and C8161 melanoma cells. Legumain activity was measureable in all cell lines, and of the potential legumain inhibitors, cystatin C, E/M, and F, cystatin C was the one mainly produced. All cells internalized cystatin C added to culture media, leading to increased intracellular cystatin C levels by 120-200%. Cystatin E/M was internalized as well but at a modest rate. The effects on intracellular legumain activity were nevertheless pronounced, probably because the cells lacked this inhibitor, and its affinity for legumain is 100-fold higher than that of cystatin C. Likewise, the low-degree uptake resulted in reduced migration and invasion of A375 cells in Matrigel to an extent comparable with the W106F variant of cystatin C with optimal uptake properties and resulting in much higher intracellular levels. Thus, cystatin E/M appears to be a good candidate to efficiently down-regulate the increased legumain activity, possibly important for the malignant phenotype of melanoma cells.

Computing UV/vis spectra using a combined molecular dynamics and quantum chemistry approach: bis-triazin-pyridine (BTP) ligands studied in solution.

We report a combined computational and experimental study to investigate the UV/vis spectra of 2,6-bis(5,6-dialkyl-1,2,4-triazin-3-yl)pyridine (BTP) ligands in solution. In order to study molecules in solution using theoretical methods, force-field parameters for the ligand-water interaction are adjusted to ab initio quantum chemical calculations. Based on these parameters, molecular dynamics (MD) simulations are carried out from which snapshots are extracted as input to quantum chemical excitation-energy calculations to obtain UV/vis spectra of BTP ligands in solution using time-dependent density functional theory (TDDFT) employing the Tamm-Dancoff approximation (TDA). The range-separated CAM-B3LYP functional is used to avoid large errors for charge-transfer states occurring in the electronic spectra. In order to study environment effects with theoretical methods, the frozen-density embedding scheme is applied. This computational procedure allows to obtain electronic spectra calculated at the (range-separated) DFT level of theory in solution, revealing solvatochromic shifts upon solvation of up to about 0.6 eV. Comparison to experimental data shows a significantly improved agreement compared to vacuum calculations and enables the analysis of relevant excitations for the line shape in solution.

The importance of current contributions to shielding constants in density-functional theory.

The sources of error in the calculation of nuclear-magnetic-resonance shielding constants determined by density-functional theory are examined. Highly accurate Kohn-Sham wave functions are obtained from coupled-cluster electron density functions and used to define accurate-but current independent-density-functional shielding constants. These new reference values, in tandem with high-accuracy coupled-cluster shielding constants, provide a benchmark for the assessment of errors in common density-functional approximations. In particular the role of errors arising in the diamagnetic and paramagnetic terms is investigated, with particular emphasis on the role of current-dependence in the latter. For carbon and nitrogen the current correction is found to be, in some cases, larger than 10 ppm. This indicates that the absence of this correction in general purpose exchange-correlation functionals is one of the main sources of error in shielding calculations using density functional theory. It is shown that the current correction improves the shielding performance of many popular approximate DFT functionals.

Cystatin C properties crucial for uptake and inhibition of intracellular target enzymes.

To elucidate the molecular requirements for cancer cell internalization of the extracellular cysteine protease inhibitor cystatin C, 12 variants of the protein were produced and used for uptake experiments in MCF-7 cells. Variants with alterations in the cysteine cathepsin binding region ((Δ1-10)-, K5A-, R8G-, (R8G,L9G,V10G)-, (R8G,L9G,V10G,W106G)-, and W106G-cystatin C) were internalized to a very low extent compared with the wild-type inhibitor. Substitutions of N39 in the legumain binding region (N39K- and N39A-cystatin C) decreased the internalization and (R24A,R25A)-cystatin C, with substitutions of charged residues not involved in enzyme inhibition, was not taken up at all. Two variants, W106F- and K75A-cystatin C, showed that the internalization can be positively affected by engineering of the cystatin molecule. Microscopy revealed vesicular co-localization of internalized cystatin C with the lysosomal marker proteins cathepsin D and legumain. Activities of both cysteine cathepsins and legumain, possible target enzymes associated with cancer cell invasion and metastasis, were down-regulated in cell homogenates following cystatin C uptake. A positive effect on regulation of intracellular enzyme activity by a cystatin variant selected from uptake properties was illustrated by incubating cells with W106F-cystatin C. This resulted in more efficient down-regulation of intracellular legumain activity than when cells were incubated with wild-type cystatin C. Uptake experiments in prostate cancer cells corroborated that the cystatin C internalization is generally relevant and confirmed an increased uptake of W106F-cystatin C, in PC3 cells. Thus, intracellular cysteine proteases involved in cancer-promoting processes might be controled by cystatin uptake.

Sampling microcanonical ensembles of trajectories using harmonic approximation in internal coordinates.

In this paper, we modify quasiclassical harmonic sampling of microcanonical ensembles of trajectories by using the curvilinear internal coordinates. The harmonic approximation in the curvilinear normal coordinates provides a more realistic description of the PES than in the conventional rectilinear ones at finite displacements. Therefore, the sampling of vibrations in the internal coordinates significantly improves the quality of the sampling in a block-box fashion, providing more realistic displacements and reducing the errors in the potential energy. In particular, the sampling of large-amplitude torsion vibrations, which is non-realistic in the Cartesian modes, becomes accurate or acceptable in the curvilinear modes.

Differentiable but exact formulation of density-functional theory.

The universal density functional F of density-functional theory is a complicated and ill-behaved function of the density-in particular, F is not differentiable, making many formal manipulations more complicated. While F has been well characterized in terms of convex analysis as forming a conjugate pair (E, F) with the ground-state energy E via the Hohenberg-Kohn and Lieb variation principles, F is nondifferentiable and subdifferentiable only on a small (but dense) subset of its domain. In this article, we apply a tool from convex analysis, Moreau-Yosida regularization, to construct, for any ε > 0, pairs of conjugate functionals ((ε)E, (ε)F) that converge to (E, F) pointwise everywhere as ε → 0(+), and such that (ε)F is (Fréchet) differentiable. For technical reasons, we limit our attention to molecular electronic systems in a finite but large box. It is noteworthy that no information is lost in the Moreau-Yosida regularization: the physical ground-state energy E(v) is exactly recoverable from the regularized ground-state energy (ε)E(v) in a simple way. All concepts and results pertaining to the original (E, F) pair have direct counterparts in results for ((ε)E, (ε)F). The Moreau-Yosida regularization therefore allows for an exact, differentiable formulation of density-functional theory. In particular, taking advantage of the differentiability of (ε)F, a rigorous formulation of Kohn-Sham theory is presented that does not suffer from the noninteracting representability problem in standard Kohn-Sham theory.

Analytic calculations of hyper-Raman spectra from density functional theory hyperpolarizability gradients.

We present the first analytic calculations of the geometrical gradients of the first hyperpolarizability tensors at the density-functional theory (DFT) level. We use the analytically calculated hyperpolarizability gradients to explore the importance of electron correlation effects, as described by DFT, on hyper-Raman spectra. In particular, we calculate the hyper-Raman spectra of the all-trans and 11-cis isomers of retinal at the Hartree-Fock (HF) and density-functional levels of theory, also allowing us to explore the sensitivity of the hyper-Raman spectra on the geometrical characteristics of these structurally related molecules. We show that the HF results, using B3LYP-calculated vibrational frequencies and force fields, reproduce the experimental data for all-trans-retinal well, and that electron correlation effects are of minor importance for the hyper-Raman intensities.

Analytic cubic and quartic force fields using density-functional theory.

We present the first analytic implementation of cubic and quartic force constants at the level of Kohn-Sham density-functional theory. The implementation is based on an open-ended formalism for the evaluation of energy derivatives in an atomic-orbital basis. The implementation relies on the availability of open-ended codes for evaluation of one- and two-electron integrals differentiated with respect to nuclear displacements as well as automatic differentiation of the exchange-correlation kernels. We use generalized second-order vibrational perturbation theory to calculate the fundamental frequencies of methane, ethane, benzene, and aniline, comparing B3LYP, BLYP, and Hartree-Fock results. The Hartree-Fock anharmonic corrections agree well with the B3LYP corrections when calculated at the B3LYP geometry and from B3LYP normal coordinates, suggesting that the inclusion of electron correlation is not essential for the reliable calculation of cubic and quartic force constants.

Internal-to-Cartesian back transformation of molecular geometry steps using high-order geometric derivatives.

In geometry optimizations and molecular dynamics calculations, it is often necessary to transform a geometry step that has been determined in internal coordinates to Cartesian coordinates. A new method for performing such transformations, the high-order path-expansion (HOPE) method, is here presented. The new method treats the nonlinear relation between internal and Cartesian coordinates by means of automatic differentiation. The method is reliable, applicable to any system of internal coordinates, and computationally more efficient than the traditional method of iterative back transformations. As a bonus, the HOPE method determines not just the Cartesian step vector but also a continuous step path expressed in the form of a polynomial, which is useful for determining reaction coordinates, for integrating trajectories, and for visualization.

Point-of-care biomarkers for prediction of kidney function trajectory among sugarcane cutters: a comparative test accuracy study.

OBJECTIVES: Heat-stressed Mesoamerican workers, such as sugarcane cutters, suffer from high rates of chronic kidney disease of non-traditional origin (CKDnt). We aimed to identify easily available early markers of rapid kidney function decline in a population at high risk of CKDnt. DESIGN: The accuracy of different biomarkers measured during harvest for prediction of cross-harvest kidney function decline were assessed in an exploratory study group, and the performance of the most promising biomarker was then assessed in an independent confirmation group. SETTING: Male sugarcane cutters in El Salvador and Nicaragua. PARTICIPANTS: 39 male Salvadoran sugarcane cutters sampled fortnightly at ≤9 occasions before and after work shift during harvest. 371 male Nicaraguan sugarcane cutters were sampled as part of routine monitoring during two harvests. Cutters worked at high physical intensity at wet-bulb globe temperatures mostly above 29°C for 6-8 hours per day 6 days a week during the 5-6 months harvest season. PRIMARY OUTCOMES: Change in estimated glomerular filtration rate (CKD Epidemiology Collaboration) across the harvest season (ΔeGFRcross-harvest). RESULTS: Dipstick leukocyturia after work shift in the El Salvadoran group was the most promising marker, explaining >25% of ΔeGFRcross-harvest variance at 8/9 occasions during harvest. Leukocyturia was associated with experiencing fever, little or dark urine, cramps, headache, dizziness and abdominal pain in the preceding 2-week period. Decreasing blood haemoglobin (Hb) and eGFR during harvest were also predictive of ΔeGFRcross-harvest. In the Nicaraguan confirmation dataset, those having ≥++ leukocyturia at any sampling during harvest had a 13 mL/min/1.73 m2 (95% CI 10 to 16 mL/min/1.73 m2) worse ΔeGFRcross-harvest than those without recorded leukocyturia. CONCLUSION: Leukocyturia and Hb, both measurable with point-of-care methods, may be early indicators for kidney injury and risk for eGFR decline among heat-stressed male workers, thereby facilitating individual-level prevention and research aiming to understand the causes of CKDnt.

Electronic spectroscopy of UO2(2+), NUO(+) and NUN: an evaluation of time-dependent density functional theory for actinides.

The performance of the time-dependent density functional theory (TDDFT) approach has been evaluated for the electronic spectrum of the UO(2)(2+), NUO(+) and NUN molecules. Different exchange-correlation functionals (LDA, PBE, BLYP, B3LYP, PBE0, M06, M06-L, M06-2X, CAM-B3LYP) and the SAOP model potential have been investigated, as has the relative importance of the adiabatic local density approximation (ALDA) to the exchange-correlation kernel. The vertical excitation energies have been compared with reference data obtained using accurate wave-function theory (WFT) methods.

The ab initio calculation of molecular electric, magnetic and geometric properties.

We give an account of some recent advances in the development of ab initio methods for the calculation of molecular response properties, involving electric, magnetic, and geometric perturbations. Particular attention is given to properties in which the basis functions depend explicitly both on time and on the applied perturbations such as perturbations involving nuclear displacements or external magnetic fields when London atomic orbitals are used. We summarize a general framework based on the quasienergy for the calculation of arbitrary-order molecular properties using the elements of the density matrix in the atomic-orbital basis as the basic variables. We demonstrate that the necessary perturbed density matrices of arbitrary order can be determined from a set of linear equations that have the same formal structure as the set of linear equations encountered when determining the linear response equations (or time-dependent self-consistent-field equations). Additional components needed to calculate properties involving perturbation-dependent basis sets are flexible one- and two-electron integral techniques for geometric or magnetic-field differentiated integrals; in Kohn-Sham density-functional theory (KS-DFT), we also need to calculate derivatives of the exchange-correlation functional. We describe a recent proposal for evaluating these contributions based on automatic differentiation. Within this framework, it is now possible to calculate any molecular property for an arbitrary self-consistent-field reference state, including two- and four-component relativistic self-consistent-field wave functions. Examples of calculations that can be performed with this formulation are presented.