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.

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.

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.

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.

Pathophysiological Mechanisms by which Heat Stress Potentially Induces Kidney Inflammation and Chronic Kidney Disease in Sugarcane Workers.

BACKGROUND: Chronic kidney disease of non-traditional origin (CKDnt) is common among Mesoamerican sugarcane workers. Recurrent heat stress and dehydration is a leading hypothesis. Evidence indicate a key role of inflammation. METHODS: Starting in sports and heat pathophysiology literature, we develop a theoretical framework of how strenuous work in heat could induce kidney inflammation. We describe the release of pro-inflammatory substances from a leaky gut and/or injured muscle, alone or in combination with tubular fructose and uric acid, aggravation by reduced renal blood flow and increased tubular metabolic demands. Then, we analyze longitudinal data from >800 sugarcane cutters followed across harvest and review the CKDnt literature to assess empirical support of the theoretical framework. RESULTS: Inflammation (CRP elevation and fever) and hyperuricemia was tightly linked to kidney injury. Rehydrating with sugary liquids and NSAID intake increased the risk of kidney injury, whereas electrolyte solution consumption was protective. Hypokalemia and hypomagnesemia were associated with kidney injury. DISCUSSION: Heat stress, muscle injury, reduced renal blood flow and fructose metabolism may induce kidney inflammation, the successful resolution of which may be impaired by daily repeating pro-inflammatory triggers. We outline further descriptive, experimental and intervention studies addressing the factors identified in this study.

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.

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.

Intervention to diminish dehydration and kidney damage among sugarcane workers.

Objective The aim of this study was to assess the potential to reduce kidney function damage during the implementation of a water, rest, shade (WRS) and efficiency intervention program among sugarcane workers. Methods A WRS intervention program adapted from the US Occupational Safety and Health Administration (OSHA) coupled with an efficiency program began two months into the 5-month harvest. One of the two groups of workers studied was provided with portable water reservoirs, mobile shaded tents, and scheduled rest periods. Health data (anthropometric and questionnaires), blood, and urine were collected at baseline and at three subsequent times over the course of the harvest. Daily wet bulb globe temperatures (WBGT) were recorded. Results Across a working day there were changes in biomarkers indicating dehydration (urine osmolality) and serum albumin and reduced estimated glomerular filtration rate (eGFR). Cross-shift eGFR decrease was present in both groups; -10.5 mL/min/1.73m2 [95% confidence interval (95% CI) -11.8- -9.1], but smaller for the intervention group after receiving the program. Decreased eGFR over the 5-month harvest was seen in both groups: in the one receiving the intervention -3.4 mL/min/1.73m 2(95% CI -5.5- -1.3) and in the other -5.3 (95% CI -7.9- -2.7). The decrease appeared to halt after the introduction of the intervention in the group receiving the program. Conclusion A WRS and efficiency intervention program was successfully introduced for workers in sugarcane fields and appears to reduce the impact of heat stress on acute and over-harvest biomarkers of kidney function. Further research is needed to determine whether biomarker changes predict reduced risk of chronic kidney disease in this type of work.

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.

Current Density Functional Theory Using Meta-Generalized Gradient Exchange-Correlation Functionals.

We present the self-consistent implementation of current-dependent (hybrid) meta-generalized gradient approximation (mGGA) density functionals using London atomic orbitals. A previously proposed generalized kinetic energy density is utilized to implement mGGAs in the framework of Kohn-Sham current density functional theory (KS-CDFT). A unique feature of the nonperturbative implementation of these functionals is the ability to seamlessly explore a wide range of magnetic fields up to 1 au (∼235 kT) in strength. CDFT functionals based on the TPSS and B98 forms are investigated, and their performance is assessed by comparison with accurate coupled-cluster singles, doubles, and perturbative triples (CCSD(T)) data. In the weak field regime, magnetic properties such as magnetizabilities and nuclear magnetic resonance shielding constants show modest but systematic improvements over generalized gradient approximations (GGA). However, in the strong field regime, the mGGA-based forms lead to a significantly improved description of the recently proposed perpendicular paramagnetic bonding mechanism, comparing well with CCSD(T) data. In contrast to functionals based on the vorticity, these forms are found to be numerically stable, and their accuracy at high field suggests that the extension of mGGAs to CDFT via the generalized kinetic energy density should provide a useful starting point for further development of CDFT approximations.

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.

The Dalton quantum chemistry program system.

Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, Møller-Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.

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.

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 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.