Higher-order molecular properties and excitation energies in single-reference and multireference coupled-cluster theory

Coupled-cluster (CC) theory is one of the most successful approaches in high-accuracy quantum chemistry. The present thesis makes a number of contributions to the determination of molecular properties and excitation energies within the CC framework. The multireference CC (MRCC) method proposed by Mukherjee and coworkers (Mk-MRCC) has been benchmarked within the singles and doubles approximation (Mk-MRCCSD) for molecular equilibrium structures. It is demonstrated that Mk-MRCCSD yields reliable results for multireference cases where single-reference CC methods fail. At the same time, the present work also illustrates that Mk-MRCC still suffers from a number of theoretical problems and sometimes gives rise to results of unsatisfactory accuracy. To determine polarizability tensors and excitation spectra in the MRCC framework, the Mk-MRCC linear-response function has been derived together with the corresponding linear-response equations. Pilot applications show that Mk-MRCC linear-response theory suffers from a severe problem when applied to the calculation of dynamic properties and excitation energies: The Mk-MRCC sufficiency conditions give rise to a redundancy in the Mk-MRCC Jacobian matrix, which entails an artificial splitting of certain excited states. This finding has established a new paradigm in MRCC theory, namely that a convincing method should not only yield accurate energies, but ought to allow for the reliable calculation of dynamic properties as well. In the context of single-reference CC theory, an analytic expression for the dipole Hessian matrix, a third-order quantity relevant to infrared spectroscopy, has been derived and implemented within the CC singles and doubles approximation. The advantages of analytic derivatives over numerical differentiation schemes are demonstrated in some pilot applications.

Method of Detecting System Function by Measuring Frequency Response

Real-time battery impedance spectrum is acquired using a one-time record. Fast Summation Transformation (FST) is a parallel method of acquiring a real-time battery impedance spectrum using a one-time record that enables battery diagnostics. An excitation current to a battery is a sum of equal amplitude sine waves of frequencies that are octave harmonics spread over a range of interest. A sample frequency is also octave and harmonically related to all frequencies in the sum. The time profile of this signal has a duration that is a few periods of the lowest frequency. The voltage response of the battery, average deleted, is the impedance of the battery in the time domain. Since the excitation frequencies are known and octave and harmonically related, a simple algorithm, FST, processes the time record by rectifying relative to the sine and cosine of each frequency. Another algorithm yields real and imaginary components for each frequency.

Method of Detecting System Function by Measuring Frequency Response

Methods of rapidly measuring an impedance spectrum of an energy storage device in-situ over a limited number of logarithmically distributed frequencies are described. An energy storage device is excited with a known input signal, and aresponse is measured to ascertain the impedance spectrum. An excitation signal is a limited time duration sum-of-sines consisting of a select number offrequencies. In one embodiment, magnitude and phase of each frequency ofinterest within the sum-of-sines is identified when the selected frequencies and sample rate are logarithmic integer steps greater than two. This technique requires a measurement with a duration of one period of the lowest frequency. In another embodiment, where selected frequencies are distributed in octave steps, the impedance spectrum can be determined using a captured time record that is reduced to a half-period of the lowest frequency.

Method of Detecting System Function by Measuring Frequency Response

Real time battery impedance spectrum is acquired using one time record, Compensated Synchronous Detection (CSD). This parallel method enables battery diagnostics. The excitation current to a test battery is a sum of equal amplitude sin waves of a few frequencies spread over range of interest. The time profile of this signal has duration that is a few periods of the lowest frequency. The voltage response of the battery, average deleted, is the impedance of the battery in the time domain. Since the excitation frequencies are known, synchronous detection processes the time record and each component, both magnitude and phase, is obtained. For compensation, the components, except the one of interest, are reassembled in the time domain. The resulting signal is subtracted from the original signal and the component of interest is synchronously detected. This process is repeated for each component.

Measurement of the vascular input function in mice for DCE-MRI

DCE-MRI is an important technique in the study of small animal cancer models because its sensitivity to vascular changes opens the possibility of quantitative assessment of early therapeutic response. However, extraction of physiologically descriptive parameters from DCE-MRI data relies upon measurement of the vascular input function (VIF), which represents the contrast agent concentration time course in the blood plasma. This is difficult in small animal models due to artifacts associated with partial volume, inflow enhancement, and the limited temporal resolution achievable with MR imaging. In this work, the development of a suite of techniques for high temporal resolution, artifact resistant measurement of the VIF in mice is described. One obstacle in VIF measurement is inflow enhancement, which decreases the sensitivity of the MR signal to the presence of contrast agent. Because the traditional techniques used to suppress inflow enhancement degrade the achievable spatiotemporal resolution of the pulse sequence, improvements can be achieved by reducing the time required for the suppression. Thus, a novel RF pulse which provides spatial presaturation contemporaneously with the RF excitation was implemented and evaluated. This maximizes the achievable temporal resolution by removing the additional RF and gradient pulses typically required for suppression of inflow enhancement. A second challenge is achieving the temporal resolution required for accurate characterization of the VIF, which exceeds what can be achieved with conventional imaging techniques while maintaining adequate spatial resolution and tumor coverage. Thus, an anatomically constrained reconstruction strategy was developed that allows for sampling of the VIF at extremely high acceleration factors, permitting capture of the initial pass of the contrast agent in mice. Simulation, phantom, and in vivo validation of all components were performed. Finally, the two components were used to perform VIF measurement in the murine heart. An in vivo study of the VIF reproducibility was performed, and an improvement in the measured injection-to-injection variation was observed. This will lead to improvements in the reliability of quantitative DCE-MRI measurements and increase their sensitivity.

Different dark conformations function in color-sensitive photosignaling by the sensory rhodopsin I-HtrI complex.

The haloarchaeal phototaxis receptor sensory rhodopsin I (SRI) in complex with its transducer HtrI delivers an attractant signal from excitation with an orange photon and a repellent signal from a second near-UV photon excitation. Using a proteoliposome system with purified SRI in complex with its transducer HtrI, we identified by site-directed fluorescence labeling a site (Ser(155)) on SRI that is conformationally active in signal relay to HtrI. Using site-directed spin labeling of Ser(155)Cys with a nitroxide side chain, we detected a change in conformation following one-photon excitation such that the spin probe exhibits a splitting of the outer hyperfine extrema (2A'(zz)) significantly smaller than that of the electron paramagnetic resonance spectrum in the dark state. The dark conformations of five mutant complexes that do not discriminate between orange and near-UV excitation show shifts to lower or higher 2A'(zz) values correlated with the alterations in their motility behavior to one- and two-photon stimuli. These data are interpreted in terms of a model in which the dark complex is populated by two conformers in the wild type, one that inhibits the CheA kinase (A) and the other that activates it (R), shifted in the dark by mutations and shifted in the wild-type SRI-HtrI complex in opposite directions by one-photon and two-photon reactions.

Relation between Excitation Power Density and Er3+ Doping Yielding the Highest Absolute Upconversion Quantum Yield

The upconversion quantum yield (UCQY) is one of the most significant parameters for upconverter materials. A high UCQY is essential for a succesful integration of upconversion in many applications, such as harvesting of the solar radiation. However, little is known about which doping level of the rare-earth ions yields the highest UCQY in the different host lattices and what are the underlying causes. Here, we investigate which Er3+ doping yields the highest UCQY in the host lattices β-NaYF4 and Gd2O2S under 4I15/2 → 4I13/2 excitation. We show for both host lattices that the optimum Er3+ doping is not fixed and it actually decreases as the irradiance of the excitation increases. To find the optimum Er3+ doping for a given irradiance, we determined the peak position of the internal UCQY as a function of the average Er−Er distance. For this purpose, we used a fit on experimental data, where the average Er−Er distance was calculated from the Er3+ doping of the upconverter samples and the lattice parameters of the host materials. We observe optimum average Er−Er distances for the host lattices β-NaYF4 and Gd2O2S with differences <14% at the same irradiance levels, whereas the optimum Er3+ doping are around 2× higher for β-NaYF4 than for Gd2O2S. Estimations by extrapolation to higher irradiances indicate that the optimum average Er−Er distance converges to values around 0.88 and 0.83 nm for β-NaYF4 and Gd2O2S, respectively. Our findings point to a fundamental relationship and focusing on the average distance between the active rare-earth ions might be a very efficient way to optimize the doping of rare-earth ions with regard to the highest achievable UCQY.

Na+ channel function, regulation, structure, trafficking and sequestration.

This paper is the second of a series of three reviews published in this issue resulting from the University of California Davis Cardiovascular Symposium 2014: Systems approach to understanding cardiac excitation-contraction coupling and arrhythmias: Na(+) channel and Na(+) transport. The goal of the symposium was to bring together experts in the field to discuss points of consensus and controversy on the topic of sodium in the heart. The present review focuses on Na(+) channel function and regulation, Na(+) channel structure and function, and Na(+) channel trafficking, sequestration and complexing.

Mechanical characterisation of tungsten-1wt.% yttrium oxide as a function of temperature and atmosphere

This study evaluates the mechanical behaviour of an Y2O3-dispersed tungsten (W) alloy and compares it to a pure W reference material. Both materials were processed via mechanical alloying (MA) and subsequent hot isostatic pressing (HIP). We performed non-standard three-point bending (TPB) tests in both an oxidising atmosphere and vacuum across a temperature range from 77 K, obtained via immersion in liquid nitrogen, to 1473 K to determine the mechanical strength, yield strength and fracture toughness. This research aims to evaluate how the mechanical behaviour of the alloy is affected by oxides formed within the material at high temperatures, primarily from 873 K, when the materials undergo a massive thermal degradation. The results indicate that the alloy is brittle to a high temperature (1473 K) under both atmospheres and that the mechanical properties degrade significantly above 873 K. We also used Vickers microhardness tests and the dynamic modulus by impulse excitation technique (IET) to determine the elastic modulus at room temperature. Moreover, we performed nanoindentation tests to determine the effect of size on the hardness and elastic modulus; however, no significant differences were found. Additionally, we calculated the relative density of the samples to assess the porosity of the alloy. Finally, we analysed the microstructure and fracture surfaces of the tested materials via field emission scanning electron microscopy (FE-SEM) and transmission electron microscopy (TEM). In this way, the relationship between the macroscopic mechanical properties and micromechanisms of failure could be determined based on the temperature and oxides formed

An unusual pathway of excitation energy deactivation in carotenoids: Singlet-to-triplet conversion on an ultrafast timescale in a photosynthetic antenna

Carotenoids are important biomolecules that are ubiquitous in nature and find widespread application in medicine. In photosynthesis, they have a large role in light harvesting (LH) and photoprotection. They exert their LH function by donating their excited singlet state to nearby (bacterio)chlorophyll molecules. In photosynthetic bacteria, the efficiency of this energy transfer process can be as low as 30%. Here, we present evidence that an unusual pathway of excited state relaxation in carotenoids underlies this poor LH function, by which carotenoid triplet states are generated directly from carotenoid singlet states. This pathway, operative on a femtosecond and picosecond timescale, involves an intermediate state, which we identify as a new, hitherto uncharacterized carotenoid singlet excited state. In LH complex-bound carotenoids, this state is the precursor on the reaction pathway to the triplet state, whereas in extracted carotenoids in solution, this state returns to the singlet ground state without forming any triplets. We discuss the possible identity of this excited state and argue that fission of the singlet state into a pair of triplet states on individual carotenoid molecules constitutes the mechanism by which the triplets are generated. This is, to our knowledge, the first ever direct observation of a singlet-to-triplet conversion process on an ultrafast timescale in a photosynthetic antenna.

Mechanical bases of frequency tuning and neural excitation at the base of the cochlea: Comparison of basilar-membrane vibrations and auditory-nerve-fiber responses in chinchilla

We review the mechanical origin of auditory-nerve excitation, focusing on comparisons of the magnitudes and phases of basilar-membrane (BM) vibrations and auditory-nerve fiber responses to tones at a basal site of the chinchilla cochlea with characteristic frequency ≈ 9 kHz located 3.5 mm from the oval window. At this location, characteristic frequency thresholds of fibers with high spontaneous activity correspond to magnitudes of BM displacement or velocity in the order of 1 nm or 50 μm/s. Over a wide range of stimulus frequencies, neural thresholds are not determined solely by BM displacement but rather by a function of both displacement and velocity. Near-threshold, auditory-nerve responses to low-frequency tones are synchronous with peak BM velocity toward scala tympani but at 80–90 dB sound pressure level (in decibels relative to 20 microPascals) and at 100–110 dB sound pressure level responses undergo two large phase shifts approaching 180°. These drastic phase changes have no counterparts in BM vibrations. Thus, although at threshold levels the encoding of BM vibrations into spike trains appears to involve only relatively minor signal transformations, the polarity of auditory-nerve responses does not conform with traditional views of how BM vibrations are transmitted to the inner hair cells. The response polarity at threshold levels, as well as the intensity-dependent phase changes, apparently reflect micromechanical interactions between the organ of Corti, the tectorial membrane and the subtectorial fluid, and/or electrical and synaptic processes at the inner hair cells.

Intramembrane charge movements and excitation– contraction coupling expressed by two-domain fragments of the Ca2+ channel

To investigate the molecular basis of the voltage sensor that triggers excitation–contraction (EC) coupling, the four-domain pore subunit of the dihydropyridine receptor (DHPR) was cut in the cytoplasmic linker between domains II and III. cDNAs for the I-II domain (α1S 1–670) and the III-IV domain (α1S 701-1873) were expressed in dysgenic α1S-null myotubes. Coexpression of the two fragments resulted in complete recovery of DHPR intramembrane charge movement and voltage-evoked Ca2+ transients. When fragments were expressed separately, EC coupling was not recovered. However, charge movement was detected in the I-II domain expressed alone. Compared with I-II and III-IV together, the charge movement in the I-II domain accounted for about half of the total charge (Qmax = 3 ± 0.23 vs. 5.4 ± 0.76 fC/pF, respectively), and the half-activation potential for charge movement was significantly more negative (V1/2 = 0.2 ± 3.5 vs. 22 ± 3.4 mV, respectively). Thus, interactions between the four internal domains of the pore subunit in the assembled DHPR profoundly affect the voltage dependence of intramembrane charge movement. We also tested a two-domain I-II construct of the neuronal α1A Ca2+ channel. The neuronal I-II domain recovered charge movements like those of the skeletal I-II domain but could not assist the skeletal III-IV domain in the recovery of EC coupling. The results demonstrate that a functional voltage sensor capable of triggering EC coupling in skeletal myotubes can be recovered by the expression of complementary fragments of the DHPR pore subunit. Furthermore, the intrinsic voltage-sensing properties of the α1A I-II domain suggest that this hemi-Ca2+ channel could be relevant to neuronal function.

Formation of junctions involved in excitation-contraction coupling in skeletal and cardiac muscle.

During excitation-contraction (e-c) coupling of striated muscle, depolarization of the surface membrane is converted into Ca2+ release from internal stores. This process occurs at intracellular junctions characterized by a specialized composition and structural organization of membrane proteins. The coordinated arrangement of the two key junctional components--the dihydropyridine receptor (DHPR) in the surface membrane and the ryanodine receptor (RyR) in the sarcoplasmic reticulum--is essential for their normal, tissue-specific function in e-c coupling. The mechanisms involved in the formation of the junctions and a potential participation of DHPRs and RyRs in this process have been subject of intensive studies over the past 5 years. In this review we discuss recent advances in understanding the organization of these molecules in skeletal and cardiac muscle, as well as their concurrent and independent assembly during development of normal and mutant muscle. From this information we derive a model for the assembly of the junctions and the establishment of the precise structural relationship between DHPRs and RyRs that underlies their interaction in e-c coupling.

Modulation of Sarco/Endoplasmic Reticulum Ca2+-ATPase 2 (SERCA2) function by acetylation following the treatment with histone deacetylase inhibitor suberoylanilide hydroxamic acid (SAHA)

Background: Acetylation and deacetylation at specific lysine (K) residues is mediated by histone acetylases (HATs) and deacetylases (HDACs), respectively. HATs and HDACs act on both histone and non-histone proteins, regulating various processes, including cardiac impulse propagation. Aim of the present work was to establish whether the function of the Ca2+ ATPase SERCA2, one of the major players in Ca2+ reuptake during excitation-contraction coupling in cardiac myocytes (CMs), could be modulated by direct K acetylation. Materials and methods: HL-1 atrial mouse cells (donated by Prof. Claycomb), zebrafish and Streptozotocin-induced diabetic rat CMs were treated with the pan-inhibitor of class I and II HDACs suberanilohydroxamic acid (SAHA) for 1.5 hour. Evaluation of SERCA2 acetylation was analyzed by co-immunoprecipitation. SERCA2 activity was measured on microsomes by pyruvate/NADH coupled reaction assay. SERCA2 mutants were obtained after cloning wild-type and mutated sequences into the pCDNA3 vector and transfected into HEK cells. Ca2+ transients in CMs (loading with Fluo3-AM, field stimulation, 0.5 Hz) and in transfected HEK cells (loading with FLUO-4, caffeine pulse) were recorded. Results: Co-Immunoprecipitation experiments performed on HL-1 cells demonstrated a significant increase in the acetylation of SERCA2 after SAHA-treatment (2.5 µM, n=3). This was associated with an increase in SERCA2 activity in microsomes obtained from HL-1 cells, after SAHA exposure (n=5). Accordingly, SAHA-treatment significantly shortened the Ca2+ reuptake time of adult zebrafish CMs. Further, SAHA 2.5 nM restored to control values the recovery time of Ca2+ transients decay in diabetic rat CMs. HDAC inhibition also improved contraction parameters, such as fraction of shortening, and increased pump activity in microsomes isolated from diabetic CMs (n=4). Notably, the K464, identified by bioinformatic tools as the most probable acetylation site on human SERCA2a, was mutated into Glutamine (Q) or Arginine (R) mimicking acetylation and deacetylation respectively. Measurements of Ca2+ transients in HEK cells revealed that the substitution of K464 with R significantly delayed the transient recovery time, thus indicating that deacetylation has a negative impact on SERCA2 function. Conclusions: Our results indicate that SERCA2 function can be improved by pro-acetylation interventions and that this mechanism of regulation is conserved among species. Therefore, the present work provides the basis to open the search for novel pharmacological tools able to specifically improve SERCA2 activity in diseases where its expression and/or function is impaired, such as diabetic cardiomyopathy.

Simulation of the secondary electrons energy deposition produced by proton beams in PMMA: influence of the target electronic excitation description

We have studied the radial dependence of the energy deposition of the secondary electron generated by swift proton beams incident with energies T = 50 keV–5 MeV on poly(methylmethacrylate) (PMMA). Two different approaches have been used to model the electronic excitation spectrum of PMMA through its energy loss function (ELF), namely the extended-Drude ELF and the Mermin ELF. The singly differential cross section and the total cross section for ionization, as well as the average energy of the generated secondary electrons, show sizeable differences at T ⩽ 0.1 MeV when evaluated with these two ELF models. In order to know the radial distribution around the proton track of the energy deposited by the cascade of secondary electrons, a simulation has been performed that follows the motion of the electrons through the target taking into account both the inelastic interactions (via electronic ionizations and excitations as well as electron-phonon and electron trapping by polaron creation) and the elastic interactions. The radial distribution of the energy deposited by the secondary electrons around the proton track shows notable differences between the simulations performed with the extended-Drude ELF or the Mermin ELF, being the former more spread out (and, therefore, less peaked) than the latter. The highest intensity and sharpness of the deposited energy distributions takes place for proton beams incident with T ~ 0.1–1 MeV. We have also studied the influence in the radial distribution of deposited energy of using a full energy distribution of secondary electrons generated by proton impact or using a single value (namely, the average value of the distribution); our results show that differences between both simulations become important for proton energies larger than ~0.1 MeV. The results presented in this work have potential applications in materials science, as well as hadron therapy (due to the use of PMMA as a tissue phantom) in order to properly consider the generation of electrons by proton beams and their subsequent transport and energy deposition through the target in nanometric scales.