Winners vs. losers: Schistosoma mansoni intestinal and liver eggs exhibit striking differences in gene expression and immunogenicity.

The eggs of the blood fluke Schistosoma mansoni are the main cause of the clinical manifestations of chronic schistosomiasis. After laying, the egg "winners" attach to the endothelium of the mesenteric vein and, after a period of development, induce the growth of a small granuloma, which facilitates their passage to the intestinal lumen. Egg "losers" carried by the bloodstream to non-specific tissues also undergo full development and induce large granuloma formation, but their life ends there. Although these trapped eggs represent a dead end in the parasite life cycle, the vast majority of studies attempting to describe the biology of the S. mansoni eggs have studied these liver-trapped "losers" instead of migrating intestinal "winners". This raises the fundamental question of how these eggs differ. With robust comparative transcriptomic analysis performed on S. mansoni eggs isolated 7 weeks post infection, we show that gene expression is critically dependent on tissue localization, both in the early and late stages of development. While mitochondrial genes and venom allergen-like proteins are significantly upregulated in mature intestinal eggs, well-described egg immunomodulators IPSE/alpha-1 and omega-1, together with micro-exon genes, are predominantly expressed in liver eggs. In addition, several proteases and protease inhibitors previously implicated in egg-host interactions display clear tissue-specific gene expression patterns. These major differences in gene expression could be then reflected in the observed different ability of liver and intestinal soluble egg antigens to elicit host immune responses and in the shorter viability of miracidia hatched from liver eggs. Our comparative analysis provides a new perspective on the biology of parasite's eggs in the context of their development and tissue localization. These findings could contribute to a broader and more accurate understanding of parasite eggs interactions with the host, which have historically been often restricted to liver eggs and sometimes inaccurately generalized.

Correlations among different platelet aggregation pathways in a group of healthy volunteers.

Platelet aggregation is a complicated process mediated by different signaling pathways. As the process is highly complex and apparently redundant, the relationships between these pathways are not yet fully known. The aim of this project was to study the interconnections among seven different aggregation pathways in a group of 53 generally healthy volunteers aged 20 to 66 years. Platelet aggregation was induced with thrombin receptor activating peptide 6 (TRAP), arachidonic acid (AA), platelet activating factor 16 (PAF), ADP, collagen, thromboxane A2 analogue U46619 or ristocetin (platelet agglutination) ex vivo in fasting blood samples according to standardized timetable protocol. Additionally, some samples were pre-treated with known clinically used antiplatelet drugs (vorapaxar, ticagrelor or acetylsalicylic acid (ASA)). Significant correlations among all used inducers were detected (Pearson correlation coefficients (rP): 0.3 to 0.85). Of all the triggers, AA showed to be the best predictor of the response to other inducers with rP ranging from 0.66 to 0.85. Interestingly, the antiplatelet response to ticagrelor strongly predicted the response to unrelated drug vorapaxar (rP = 0.71). Our results indicate that a response to one inducer can predict the response for other triggers or even to an antiplatelet drug. These data are useful for future testing but should be also confirmed in patients.

What is the context?• Platelet activation is a complicated process with multiple signaling cascades involved.• A total of seven common platelet triggers (ADP, collagen, TRAP-6, PAF, arachidonic acid/AA/, ristocetin and U46619) were tested.• The process is dependent on many factors including sex, age, concomitant disease(s), pharmacotherapy.What is new?• There were significant correlations between all tested aggregatory cascades.• AA has the highest rate of response predictability in our heterogeneous generally healthy volunteer group.• There was no correlation between impedance aggregometry in whole blood and turbidimetric measurement with platelet-rich plasma.What is the impact?• The effect of antiplatelet drugs can be assessed from the reaction to different trigger(s) at least in this group of healthy patients.• Future studies must test these relationships in patients with different diseases.

Uncovering the essential roles of glutamate carboxypeptidase 2 orthologs in Caenorhabditis elegans.

Human glutamate carboxypeptidase 2 (GCP2) from the M28B metalloprotease group is an important target for therapy in neurological disorders and an established tumor marker. However, its physiological functions remain unclear. To better understand general roles, we used the model organism Caenorhabditis elegans to genetically manipulate its three existing orthologous genes and evaluate the impact on worm physiology. The results of gene knockout studies showed that C. elegans GCP2 orthologs affect the pharyngeal physiology, reproduction, and structural integrity of the organism. Promoter-driven GFP expression revealed distinct localization for each of the three gene paralogs, with gcp-2.1 being most abundant in muscles, intestine, and pharyngeal interneurons, gcp-2.2 restricted to the phasmid neurons, and gcp-2.3 located in the excretory cell. The present study provides new insight into the unique phenotypic effects of GCP2 gene knockouts in C. elegans, and the specific tissue localizations. We believe that elucidation of particular roles in a non-mammalian organism can help to explain important questions linked to physiology of this protease group and in extension to human GCP2 involvement in pathophysiological processes.

Cost-Efficient High-Resolution Linear Absorption Spectra through Extrapolating the Dipole Moment from Real-Time Time-Dependent Electronic-Structure Theory.

We present a novel function fitting method for approximating the propagation of the time-dependent electric dipole moment from real-time electronic structure calculations. Real-time calculations of the electronic absorption spectrum require discrete Fourier transforms of the electric dipole moment. The spectral resolution is determined by the total propagation time, i.e., the trajectory length of the dipole moment, causing a high computational cost. Our developed method uses function fitting on shorter trajectories of the dipole moment, achieving arbitrary spectral resolution through extrapolation. Numerical testing shows that the fitting method can reproduce high-resolution spectra by using short dipole trajectories. The method converges with as little as 100 a.u. dipole trajectories for some systems, though the difficulty converging increases with the spectral density. We also introduce an error estimate of the fit, reliably assessing its convergence and hence the quality of the approximated spectrum.

The Impact of Convertase Subtilisin/Kexin Type 9 Monoclonal Antibodies with and without Apheresis on Platelet Aggregation in Familial Hypercholesterolemia.

BACKGROUND AND AIMS: It is well known that elevated cholesterol is associated with enhanced platelet aggregation and patients suffering from familial hypercholesterolemia (FH) have a high risk of thrombotic cardiovascular events. Although decreasing cholesterol level is associated with attenuation of platelet hyperactivity, there are currently no data on the effect of convertase subtilisin/kexin type 9 monoclonal antibodies (PCSK9ab) on platelet reactivity in FH. The aim of the study was to analyse the impact of different therapies including PCSK9ab on platelet aggregation in FH. METHODS: This study enrolled all 15 patients treated in the University Hospital Hradec Králové for FH. PCSK9ab have been administered in 12 of 15 patients while 8 patients were also undergoing lipid apheresis. Blood samples from all patients including pre- and post-apheresis period were tested for platelet aggregation triggered by 7 inducers, and the effect of 3 clinically used drugs (acetylsalicylic acid, ticagrelor and vorapaxar) was compared as well. RESULTS: Although apheresis decreased the reactivity of platelets in general, platelet responses were not different between non-apheresis patients treated with PCSK9ab and apheresis patients (post-apheresis values) with the exception of ristocetin. However, when compared to age-matched healthy population, FH patients had significantly lower platelet aggregation responses to 4 out of 7 used inducers and higher profit from 2 out of 3 used antiplatelet drugs even after exclusion of FH patients regularly receiving conventional antiplatelet treatment. CONCLUSION: This study showed for the first time the suitability of PCSK9ab treatment for reduction of platelet reactivity in FH patients.

The Effect of 4-Methylcatechol on Platelets in Familial Hypercholesterolemic Patients Treated with Lipid Apheresis and/or Proprotein Convertase Subtilisin Kexin 9 Monoclonal Antibodies.

Elevated low-density lipoprotein (LDL) cholesterol levels lead to atherosclerosis and platelet hyperaggregability, both of which are known culprits of arterial thrombosis. Normalization of LDL cholesterol in familial hypercholesterolemia (FH) is not an easy task and frequently requires specific treatment, such as regularly performed lipid apheresis and/or novel drugs such as proprotein convertase subtilisin kexin 9 monoclonal antibodies (PCSK9Ab). Moreover, a high resistance rate to the first-line antiplatelet drug acetylsalicylic acid (ASA) stimulated research of novel antiplatelet drugs. 4-methylcatechol (4-MC), a known metabolite of several dietary flavonoids, may be a suitable candidate. The aim of this study was to analyse the antiplatelet effect of 4-MC in FH patients and to compare its impact on two FH treatment modalities via whole-blood impedance aggregometry. When compared to age-matched, generally healthy controls, the antiplatelet effect of 4-MC against collagen-induced aggregation was higher in FH patients. Apheresis itself improved the effect of 4-MC on platelet aggregation and blood from patients treated with this procedure and pretreated with 4-MC had lower platelet aggregability when compared to those solely treated with PCKS9Ab. Although this study had some inherent limitations, e.g., a low number of patients and possible impact of administered drugs, it confirmed the suitability of 4-MC as a promising antiplatelet agent and also demonstrated the effect of 4-MC in patients with a genetic metabolic disease for the first time.

Exact Two-Component TDDFT with Simple Two-Electron Picture-Change Corrections: X-ray Absorption Spectra Near L- and M-Edges of Four-Component Quality at Two-Component Cost.

X-ray absorption spectroscopy (XAS) has gained popularity in recent years as it probes matter with high spatial and elemental sensitivities. However, the theoretical modeling of XAS is a challenging task since XAS spectra feature a fine structure due to scalar (SC) and spin-orbit (SO) relativistic effects, in particular near L and M absorption edges. While full four-component (4c) calculations of XAS are nowadays feasible, there is still interest in developing approximate relativistic methods that enable XAS calculations at the two-component (2c) level while maintaining the accuracy of the parent 4c approach. In this article we present theoretical and numerical insights into two simple yet accurate 2c approaches based on an (extended) atomic mean-field exact two-component Hamiltonian framework, (e)amfX2C, for the calculation of XAS using linear eigenvalue and damped response time-dependent density functional theory (TDDFT). In contrast to the commonly used one-electron X2C (1eX2C) Hamiltonian, both amfX2C and eamfX2C account for the SC and SO two-electron and exchange-correlation picture-change (PC) effects that arise from the X2C transformation. As we demonstrate on L- and M-edge XAS spectra of transition metal and actinide compounds, the absence of PC corrections in the 1eX2C approximation results in a substantial overestimation of SO splittings, whereas (e)amfX2C Hamiltonians reproduce all essential spectral features such as shape, position, and SO splitting of the 4c references in excellent agreement, while offering significant computational savings. Therefore, the (e)amfX2C PC correction models presented here constitute reliable relativistic 2c quantum-chemical approaches for modeling XAS.

SC-560 and mofezolac isosteres as new potent COX-1 selective inhibitors with antiplatelet effect.

Selective cyclooxygenase (COX)-1 inhibitors can be employed as potential cardioprotective drugs. Moreover, COX-1 plays a key role in inflammatory processes and its activity is associated with some types of cancer. In this work, we designed and synthesized a set of compounds that structurally mimic the selective COX-1 inhibitors, SC-560 and mofezolac, the central cores of which were replaced either with triazole or benzene rings. The advantage of this approach is a relatively simple synthesis in comparison with the syntheses of parent compounds. The newly synthesized compounds exhibited remarkable activity and selectivity toward COX-1 in the enzymatic in vitro assay. The most potent compound, 10a (IC50 = 3 nM for COX-1 and 850 nM for COX-2), was as active as SC-560 (IC50  = 2.4 nM for COX-1 and 470 nM for COX-2) toward COX-1 and it was even more selective. The in vitro COX-1 enzymatic activity was further confirmed in the cell-based whole-blood antiplatelet assay, where three out of four selected compounds (10a,c,d, and 3b) exerted outstanding IC50 values in the nanomolar range (9-252 nM). Moreover, docking simulations were performed to reveal key interactions within the COX-1 binding pocket. Furthermore, the toxicity of the selected compounds was tested using the normal human kidney HK-2 cell line.

Accurate Relativistic Real-Time Time-Dependent Density Functional Theory for Valence and Core Attosecond Transient Absorption Spectroscopy.

First principles theoretical modeling of out-of-equilibrium processes observed in attosecond pump-probe transient absorption spectroscopy (TAS) triggering pure electron dynamics remains a challenging task, especially for heavy elements and/or core excitations containing fingerprints of scalar and spin-orbit relativistic effects. To address this, we formulate a methodology for simulating TAS within the relativistic real-time, time-dependent density functional theory (RT-TDDFT) framework, for both the valence and core energy regimes. Especially for TAS, full four-component (4c) RT simulations are feasible but computationally demanding. Therefore, in addition to the 4c approach, we also introduce the atomic mean-field exact two-component (amfX2C) Hamiltonian accounting for one- and two-electron picture-change corrections within RT-TDDFT. amfX2C preserves the accuracy of the parent 4c method at a fraction of its computational cost. Finally, we apply the methodology to study valence and near-L2,3-edge TAS processes of experimentally relevant systems and provide additional physical insights using relativistic nonequilibrium response theory.

Proteases and their inhibitors involved in Schistosoma mansoni egg-host interaction revealed by comparative transcriptomics with Fasciola hepatica eggs.

Schistosoma mansoni eggs are the main causative agents of the pathological manifestations of schistosomiasis. The eggs are laid in the host bloodstream, then they migrate through the intestinal wall into the lumen. However, a significant proportion of the eggs become lodged in the liver, where they cause inflammation and fibrosis. In this study, we focus on a specific group of proteins expressed by the egg, namely proteases and their inhibitors. These molecules are often involved in schistosome-host interactions, but are still unexplored in the egg stage. Using RNA-seq and comparative transcriptomics of immature and mature S. mansoni eggs, we mapped the portfolio of proteases and their inhibitors, and determined their gene expression levels. In addition, we compared these data with gene expression of proteases and their inhibitors in Fasciola hepatica eggs. Fasciola hepatica eggs served as a useful comparative model, as they do not migrate through tissues and inflict pathology. We detected transcription of 135 and 117 proteases in S. mansoni and F. hepatica eggs, respectively, with 87 identified as orthologous between the two species. In contrast, we observed only four orthologous inhibitors out of 21 and 16 identified in S. mansoni and F. hepatica eggs, respectively. Among others, we measured high and developmentally regulated levels of expression of metalloproteases in S. mansoni eggs, specifically aminopeptidase N1, endothelin-converting enzyme 1, and several leishmanolysin-like peptidases. We identified highly transcribed protease inhibitors serpin and alpha-2-macroglobulin that are unique to S. mansoni eggs, and antistasin-like inhibitor in F. hepatica eggs. This study provides new insights into the portfolio of proteases and inhibitors expressed by S. mansoni with potential roles in egg tissue migration, stimulation of angiogenesis, and interaction with host blood and immunity.

Sex-Related Differences in Platelet Aggregation: A Literature Review Supplemented with Local Data from a Group of Generally Healthy Individuals.

The process of platelet aggregation is often influenced by several factors including sex and age. A literature review confirmed the existence of sex-related differences in platelet aggregation. Although 68 out of 78 papers found such differences, there are still some controversies regarding these differences, which can be due to multiple factors (age, trigger, concomitant disease, sample handling, etc.). These outcomes are discussed in line with novel results obtained from a local study, in which blood samples from a total of 53 overall healthy women and men with ages ranging from 20 to 66 years were collected. Aggregation was induced with seven different triggers (ristocetin, thrombin receptor activating peptide 6 [TRAP-6], arachidonic acid [AA], platelet-activating factor 16 [PAF-16], ADP, collagen, or thromboxane A2 analog U-46619) ex vivo. In addition, three FDA-approved antiplatelet drugs (vorapaxar, ticagrelor, or acetylsalicylic acid [ASA]) were also tested. In general, women had higher aggregation responses to some agonists (ADP, TRAP), as well as lower benefit from inhibitors (ASA, vorapaxar). The aggregatory responses to AA and TRAP decreased with age in both sexes, while responses to ADP, U-46619, and PAF were affected by age only in women. In conclusion, more studies are needed to decipher the biological importance of sex-related differences in platelet aggregation in part to enable personalized antiplatelet treatment.

Neurovascular responses to neuronal activity during sensory development.

Understanding the development of intercellular communication in sensory regions is relevant to elucidate mechanisms of physiological and pathological responses to oxygen shortage in the newborn brain. Decades of studies in laboratory rodents show that neuronal activity impacts sensory maturation during two periods of postnatal development distinguished by the maturation of accessory structures at the sensory periphery. During the first of these developmental periods, angiogenesis is modulated by neuronal activity, and physiological levels of neuronal activity cause local tissue hypoxic events. This correlation suggests that neuronal activity is upstream of the production of angiogenic factors, a process that is mediated by intermittent hypoxia caused by neuronal oxygen consumption. In this perspective article we address three theoretical implications based on this hypothesis: first, that spontaneous activity of sensory neurons has properties that favor the generation of intermittent tissue hypoxia in neonate rodents; second, that intermittent hypoxia promotes the expression of hypoxia inducible transcription factors (HIFs) in sensory neurons and astrocytes; and third, that activity-dependent production of angiogenic factors is involved in pathological oxygen contexts.

The Antiplatelet Effect of 4-Methylcatechol in a Real Population Sample and Determination of the Mechanism of Action.

A polyphenol-rich diet has beneficial effects on cardiovascular health. However, dietary polyphenols generally have low bioavailability and reach low plasma concentrations. Small phenolic metabolites of these compounds formed by human microbiota are much more easily absorbable and could be responsible for this effect. One of these metabolites, 4-methylcatechol (4-MC), was suggested to be a potent anti-platelet compound. The effect of 4-MC was tested ex vivo in a group of 53 generally healthy donors using impedance blood aggregometry. The mechanism of action of this compound was also investigated by employing various aggregation inducers/inhibitors and a combination of aggregometry and enzyme linked immunosorbent assay (ELISA) methods. 4-MC was confirmed to be more potent than acetylsalicylic acid on both arachidonic acid and collagen-triggered platelet aggregation. Its clinically relevant effect was found even at a concentration of 10 µM. Mechanistic studies showed that 4-MC is able to block platelet aggregation caused by the stimulation of different pathways (receptors for the von Willebrand factor and platelet-activating factor, glycoprotein IIb/IIIa, protein kinase C, intracellular calcium elevation). The major mechanism was defined as interference with cyclooxygenase-thromboxane synthase coupling. This study confirmed the strong antiplatelet potential of 4-MC in a group of healthy donors and defined its mechanism of action.

Comparison of Antiplatelet Effects of Phenol Derivatives in Humans.

Flavonoids are associated with positive cardiovascular effects. However, due to their low bioavailability, metabolites are likely responsible for these properties. Recently, one of these metabolites, 4-methylcatechol, was described to be a very potent antiplatelet compound. This study aimed to compare its activity with its 22 close derivatives both of natural or synthetic origin in order to elucidate a potential structure-antiplatelet activity relationship. Blood from human volunteers was induced to aggregate by arachidonic acid (AA), collagen or thrombin, and plasma coagulation was also studied. Potential toxicity was tested on human erythrocytes as well as on a cancer cell line. Our results indicated that 17 out of the 22 compounds were very active at a concentration of 40 µM and, importantly, seven of them had an IC50 on AA-triggered aggregation below 3 µM. The effects of the most active compounds were confirmed on collagen-triggered aggregation too. None of the tested compounds was toxic toward erythrocytes at 50 µM and four compounds partly inhibited proliferation of breast cancer cell line at 100 µM but not at 10 µM. Additionally, none of the compounds had a significant effect on blood coagulation or thrombin-triggered aggregation. This study hence reports four phenol derivatives (4-ethylcatechol, 4-fluorocatechol, 2-methoxy-4-ethylphenol and 3-methylcatechol) suitable for future in vivo testing.

Accurate X-ray Absorption Spectra near L- and M-Edges from Relativistic Four-Component Damped Response Time-Dependent Density Functional Theory.

The simulation of X-ray absorption spectra requires both scalar and spin-orbit (SO) relativistic effects to be taken into account, particularly near L- and M-edges where the SO splitting of core p and d orbitals dominates. Four-component Dirac-Coulomb Hamiltonian-based linear damped response time-dependent density functional theory (4c-DR-TDDFT) calculates spectra directly for a selected frequency region while including the relativistic effects variationally, making the method well suited for X-ray applications. In this work, we show that accurate X-ray absorption spectra near L2,3- and M4,5-edges of closed-shell transition metal and actinide compounds with different central atoms, ligands, and oxidation states can be obtained by means of 4c-DR-TDDFT. While the main absorption lines do not change noticeably with the basis set and geometry, the exchange-correlation functional has a strong influence with hybrid functionals performing the best. The energy shift compared to the experiment is shown to depend linearly on the amount of Hartee-Fock exchange with the optimal value being 60% for spectral regions above 1000 eV, providing relative errors below 0.2% and 2% for edge energies and SO splittings, respectively. Finally, the methodology calibrated in this work is used to reproduce the experimental L2,3-edge X-ray absorption spectra of [RuCl2(DMSO)2(Im)2] and [WCl4(PMePh2)2], and resolve the broad bands into separated lines, allowing an interpretation based on ligand field theory and double point groups. These results support 4c-DR-TDDFT as a reliable method for calculating and analyzing X-ray absorption spectra of chemically interesting systems, advance the accuracy of state-of-the art relativistic DFT approaches, and provide a reference for benchmarking more approximate techniques.

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.

Relativistic four-component linear damped response TDDFT for electronic absorption and circular dichroism calculations.

We present a detailed theory, implementation, and a benchmark study of a linear damped response time-dependent density functional theory (TDDFT) based on the relativistic four-component (4c) Dirac-Kohn-Sham formalism using the restricted kinetic balance condition for the small-component basis and a noncollinear exchange-correlation kernel. The damped response equations are solved by means of a multifrequency iterative subspace solver utilizing decomposition of the equations according to Hermitian and time-reversal symmetry. This partitioning leads to robust convergence, and the detailed algorithm of the solver for relativistic multicomponent wavefunctions is also presented. The solutions are then used to calculate the linear electric- and magnetic-dipole responses of molecular systems to an electric perturbation, leading to frequency-dependent dipole polarizabilities, electronic absorption, circular dichroism (ECD), and optical rotatory dispersion (ORD) spectra. The methodology has been implemented in the relativistic spectroscopy DFT program ReSpect, and its performance was assessed on a model series of dimethylchalcogeniranes, C4H8X (X = O, S, Se, Te, Po, Lv), and on larger transition metal complexes that had been studied experimentally, [M(phen)3]3+ (M = Fe, Ru, Os). These are the first 4c damped linear response TDDFT calculations of ECD and ORD presented in the literature.

Resolution-of-identity accelerated relativistic two- and four-component electron dynamics approach to chiroptical spectroscopies.

We present an implementation and application of electron dynamics based on real-time time-dependent density functional theory (RT-TDDFT) and relativistic 2-component X2C and 4-component Dirac-Coulomb (4c) Hamiltonians to the calculation of electron circular dichroism and optical rotatory dispersion spectra. In addition, the resolution-of-identity approximation for the Coulomb term (RI-J) is introduced into RT-TDDFT and formulated entirely in terms of complex quaternion algebra. The proposed methodology was assessed on the dimethylchalcogenirane series, C4H8X (X = O, S, Se, Te, Po, Lv), and the spectra obtained by non-relativistic and relativistic methods start to disagree for Se and Te, while dramatic differences are observed for Po and Lv. The X2C approach, even in its simplest one-particle form, reproduces the reference 4c results surprisingly well across the entire series while offering an 8-fold speed-up of the simulations. An overall acceleration of RT-TDDFT by means of X2C and RI-J increases with system size and approaches a factor of almost 25 when compared to the full 4c treatment, without compromising the accuracy of the final spectra. These results suggest that one-particle X2C electron dynamics with RI-J acceleration is an attractive method for the calculation of chiroptical spectra in the valence region.

Acceleration of Relativistic Electron Dynamics by Means of X2C Transformation: Application to the Calculation of Nonlinear Optical Properties.

The Liouville-von Neumann equation based on the four-component matrix Dirac-Kohn-Sham Hamiltonian is transformed to a quasirelativistic exact two-component (X2C) form and then used to solve the time evolution of the electronic states only. By this means, a significant acceleration by a factor of 7 or more has been achieved. The transformation of the original four-component equation of motion is formulated entirely in matrix algebra, following closely the X2C decoupling procedure of Ilias and Saue [ J. Chem. Phys. 2007 , 126 , 064102 ] proposed earlier for a static (time-independent) case. In a dynamic (time-dependent) regime, however, an adiabatic approximation must in addition be introduced in order to preserve the block-diagonal form of the time-dependent Dirac-Fock operator during the time evolution. The resulting X2C Liouville-von Neumann electron dynamics (X2C-LvNED) is easy to implement as it does not require an explicit form of the picture-change transformed operators responsible for the (higher-order) relativistic corrections and/or interactions with external fields. To illustrate the accuracy and performance of the method, numerical results and computational timings for nonlinear optical properties are presented. All of the time domain X2C-LvNED results show excellent agreement with the reference four-component calculations as well as with the results obtained from frequency domain response theory.

Excitation Energies from Real-Time Propagation of the Four-Component Dirac-Kohn-Sham Equation.

We report the first implementation of real-time time-dependent density functional theory (RT-TDDFT) at the relativistic four-component level of theory. In contrast to the perturbative linear-response TDDFT approach (LR-TDDFT), the RT-TDDFT approach performs an explicit time propagation of the Dirac-Kohn-Sham density matrix, offering the possibility to simulate molecular spectroscopies involving strong electromagnetic fields while, at the same time, treating relativistic scalar and spin-orbit corrections variationally. The implementation is based on the matrix representation of the Dirac-Coulomb Hamiltonian in the basis of restricted kinetically balanced Gaussian-type functions, exploiting the noncollinear Kramers unrestricted formalism implemented in the program ReSpect. We also present an analytic form for the delta-type impulse commonly used in RT-TDDFT calculations, as well as a dipole-weighted transition matrix analysis, facilitating the interpretation of spectral transitions in terms of ground-state molecular orbitals. The possibilities offered by the methodology are illustrated by investigating vertical excitation energies and oscillator strengths for ground-state to excited-state transitions in the Group 12 atoms and in heavy-element hydrides. The accuracy of the method is assessed by comparing the excitation energies obtained with earlier relativistic linear response TDDFT results and available experimental data.