Variational calculation of static and dynamic vibrational nonlinear optical properties

The vibrational configuration interaction method used to obtain static vibrational (hyper)polarizabilities is extended to dynamic nonlinear optical properties in the infinite optical frequency approximation. Illustrative calculations are carried out on H2 O and N H3. The former molecule is weakly anharmonic while the latter contains a strongly anharmonic umbrella mode. The effect on vibrational (hyper)polarizabilities due to various truncations of the potential energy and property surfaces involved in the calculation are examined

Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA

New Radicals and new applications of Dynamic Nuclear Polarization of biological systems

Dynamic Nuclear Polarization (DNP) is an emerging technique that could revolutionize the NMR study of small molecules at very low concentrations by the increase in sensitivity that results from transfer of polarization between electronic and nuclear spins. Although the underlying physics has been known for a long time, in the last few years there has been a lot of excitement on the chemistry and biology NMR community caused by the demonstration that the highly polarized nuclei that are prepared in solid state at very low temperatures (1-2 K) could be rapidly transferred to liquid samples at room temperature and studied in solution by conventional NMR techniques. In favorable cases several order of magnitude increases in sensitivity have been achieved. The technique is now mature enough that a commercial instrument is available. The efficiency of DNP depends on two crucial aspects: i) the efficiency of the nuclear polarization process and ii) the efficiency of the transfer from the initial solid state to the fluid state in which NMR is measured. The preferred areas of application (iii) will be dictated by situations in which the low concentration of the sample or its intrinsic low receptivity are the limiting factors .

Glutamatergic and GABAergic energy metabolism measured in the rat brain by (13) C NMR spectroscopy at 14.1 T.

Energy metabolism supports both inhibitory and excitatory neurotransmission processes. This study investigated the specific contribution of astrocytic metabolism to γ-aminobutyric acid (GABA) synthesis and inhibitory GABAergic neurotransmission that remained to be ilucidated in vivo. Therefore, we measured (13) C incorporation into brain metabolites by dynamic (13) C nuclear magnetic resonance spectroscopy at 14.1 T in rats under α-chloralose anaesthesia during infusion of [1,6-(13) C]glucose. The enhanced sensitivity at 14.1 T allowed to quantify incorporation of (13) C into the three aliphatic carbons of GABA non-invasively. Metabolic fluxes were determined with a mathematical model of brain metabolism comprising glial, glutamatergic and GABAergic compartments. GABA synthesis rate was 0.11 ± 0.01 μmol/g/min. GABA-glutamine cycle was 0.053 ± 0.003 μmol/g/min and accounted for 22 ± 1% of total neurotransmitter cycling between neurons and glia. Cerebral glucose oxidation was 0.47 ± 0.02 μmol/g/min, of which 35 ± 1% and 7 ± 1% was diverted to the glutamatergic and GABAergic tricarboxylic acid cycles, respectively. The remaining fraction of glucose oxidation was in glia, where 12 ± 1% of the TCA cycle flux was dedicated to oxidation of GABA. 16 ± 2% of glutamine synthesis was provided to GABAergic neurons. We conclude that substantial metabolic activity occurs in GABAergic neurons and that glial metabolism supports both glutamatergic and GABAergic neurons in the living rat brain. We performed (13) C NMR spectroscopy in vivo at high magnetic field (14.1 T) upon administration of [1,6-(13) C]glucose. This allowed to measure (13) C incorporation into the three aliphatic carbons of GABA in the rat brain, in addition to those of glutamate, glutamine and aspartate. These data were then modelled to determine fluxes of energy metabolism in GABAergic and glutamatergic neurons and glial cells.

Coherent tunnelling in Cu2+- and Ag2+-doped MgO and CaO: Cu2+ explored through ab initio calculations

The observation of coherent tunnelling in Cu2+ - and Ag2+ -doped MgO and CaO:Cu2+ was a crucial discovery in the realm of the Jahn-Teller (JT) effect. The main reasons favoring this dynamic behavior are now clarified through ab initio calculations on Cu2+ - and Ag2+ -doped cubic oxides. Small JT distortions and an unexpected low anharmonicity of the eg JT mode are behind energy barriers smaller than 25 cm-1 derived through CASPT2 calculations for Cu2+ - and Ag2+ -doped MgO and CaO:Cu2+ . The low anharmonicity is shown to come from a strong vibrational coupling of MO610- units (M=Cu,Ag) to the host lattice. The average distance between the d9 impurity and ligands is found to vary significantly on passing from MgO to SrO following to a good extent the lattice parameter.

In vivo 13C NMR spectroscopy and metabolic modeling in the brain: a practical perspective.

In vivo 13C NMR spectroscopy has the unique capability to measure metabolic fluxes noninvasively in the brain. Quantitative measurements of metabolic fluxes require analysis of the 13C labeling time courses obtained experimentally with a metabolic model. The present work reviews the ingredients necessary for a dynamic metabolic modeling study, with particular emphasis on practical issues.

Dynamic and quasistatic trajectories in quasifission reactions and particle emission

We show that the quasifission paths predicted by the one-body dissipation dynamics, in the slowest phase of a binary reaction, follow a quasistatic path, which represents a sequence of states of thermal equilibrium at a fixed value of the deformation coordinate. This establishes the use of the statistical particle-evaporation model in the case of dynamical time-evolving systems. Pre- and post-scission multiplicities of neutrons and total multiplicities of protons and α particles in fission reactions of 63Cu+92Mo, 60Ni+100Mo, 63Cu+100Mo at 10 MeV/u and 20Ne+144,148,154Sm at 20 MeV/u are reproduced reasonably well with statistical model calculations performed along dynamic trajectories whose slow stage (from the most compact configuration up to the point where the neck starts to develop) lasts some 35×10−21 s.

The Dynamic Customer Portfolio Method in Customer Relationship Management

Tutkielman tavoitteena on määritellä keskeiset ja sopivat asiakasportfoliomallit ja asiakasmatriisit asiakassuhteen määrittämiseen. Tutkimus keskittyy asiakassuhteen arvottamiseen ja avainasiakkaiden määrittämiseen kohdeyrityksessä. Keskeisimmät ja sopivimmat asiakasportfliomallit huomioidaan asiakkaiden arvioinnissa. Tutkielman teoriaosassa esitellään tunnetuimmat ja käytetyimmät asiakasportfoliomallit ja matriisit alan kirjallisuuden perusteella. Tämän lisäksi asiakasportfoliomalleihin yhdistetään näkökulmia suhdemarkkinoinnin, asiakkuuksien johtamisen ja tuoteportfolioiden teorioista. Keskeisimmät kirjallisuuden lähteet ovat johtamisen ja markkinoinnin alalta. Tutkielman empiriaosassa esitellään kohdeyritys ja sen tämän hetkinen asiakassuhteiden johtamiskäytäntö. Lisäksi tehdään parannusehdotuksia kohdeyrityksen nykyiseen asiakassuhteiden arvottamismenetelmään jotta asiakassuhteiden arvon laskeminen vastaisi mahdollisimman hyvin kohdeyrityksen nykyisiä tarpeita. Asiakassuhteen arvon määrittämiseksi käytetään myös fokusryhmähaastattelua. Avainasiakkaat määritellään ja tilannetta havainnollistetaan sijoittamalla avainasiakkaat asiakasportfolioon.

NMR chemical shifts of the rhodopsin chromophore in the dark state and in bathorhodopsin: a hybrid QM/MM molecular dynamics study.

We investigate nuclear magnetic resonance (NMR) parameters of the rhodopsin chromophore in the dark state of the protein and in the early photointermediate bathorhodopsin via first-principles molecular dynamics simulations and NMR chemical shift calculations in a hybrid quantum/classical (QM/MM) framework. NMR parameters are particularly sensitive to structural properties and to the chemical environment, which allows us to address different questions about the retinal chromophore in situ. Our calculations show that both the 13C and the 1H NMR chemical shifts are rather insensitive to the protonation state of Glu181, an ionizable amino acid side chain located in the vicinity of the isomerizing 11-cis bond. Thus, other techniques should be better suited to establish its protonation state. The calculated chemical shifts for bathorhodopsin further support our previously published theoretical structure, which is in very good agreement with more recent X-ray data.

Dynamic electrostatic force microscopy in liquid media

We present the implementation of dynamic electrostatic force microscopy in liquid media. This implementation enables the quantitative imaging of local dielectric properties of materials in electrolyte solutions with nanoscale spatial resolution. Local imaging capabilities are obtained by probing the frequency-dependent and ionic concentration-dependent electrostatic forces at high frequency (>1 MHz), while quantification of the interaction forces is obtained with finite-element numerical calculations. The results presented open a wide range of possibilities in a number of fields where the dielectric properties of materials need to be probed at the nanoscale and in a liquid environment.

Dynamic electrostatic force microscopy in liquid media

We present the implementation of dynamic electrostatic force microscopy in liquid media. This implementation enables the quantitative imaging of local dielectric properties of materials in electrolyte solutions with nanoscale spatial resolution. Local imaging capabilities are obtained by probing the frequency-dependent and ionic concentration-dependent electrostatic forces at high frequency (>1 MHz), while quantification of the interaction forces is obtained with finite-element numerical calculations. The results presented open a wide range of possibilities in a number of fields where the dielectric properties of materials need to be probed at the nanoscale and in a liquid environment.

Comparative analysis of the trypanocidal activity and chemical properties of E-lychnophoric acid and its derivatives using theoretical calculations

E-Lychnophoric acid 1, its derivative ester 2 and alcohol 3 killed 100% of trypomastigote blood forms of Trypanosoma cruzi at the concentrations of 13.86, 5.68, and 6.48 µg/mL, respectively. Conformational distribution calculations (AM1) of 1, 2 and 3 gave minimum energies for the conformers a, b, c, and d, which differ from each other only in the cyclononene ring geometry. Calculations (DFT/BLYP/6-31G*) of geometry optimization and chemical properties were performed for conformers of 1, 2, and 3. The theoretical results were numerically compared to the trypanocidal activity. Calculated values of atomic charge, orbital population, and vibrational frequencies showed that the C-4-C-5 pi-endocyclic bond does not affect the trypanocidal activity of the studied compounds. Nevertheless, the structure of the group at C-4 strongly influences the activity. However, the theoretical results indicated that the intra-ring (C-1 and C-9) and pi-exocycle (C-8 and C-14) carbons of caryophyllene-type structures promote the trypanocidal activity of these compounds.

Numerical Investigation of Flow Past a Circular Cylinder and in a Staggered Tube Bundle Using Various Turbulence Models

Transitional flow past a three-dimensional circular cylinder is a widely studied phenomenon since this problem is of interest with respect to many technical applications. In the present work, the numerical simulation of flow past a circular cylinder, performed by using a commercial CFD code (ANSYS Fluent 12.1) with large eddy simulation (LES) and RANS (κ - ε and Shear-Stress Transport (SST) κ - ω! model) approaches. The turbulent flow for ReD = 1000 & 3900 is simulated to investigate the force coefficient, Strouhal number, flow separation angle, pressure distribution on cylinder and the complex three dimensional vortex shedding of the cylinder wake region. The numerical results extracted from these simulations have good agreement with the experimental data (Zdravkovich, 1997). Moreover, grid refinement and time-step influence have been examined. Numerical calculations of turbulent cross-flow in a staggered tube bundle continues to attract interest due to its importance in the engineering application as well as the fact that this complex flow represents a challenging problem for CFD. In the present work a time dependent simulation using κ – ε, κ - ω! and SST models are performed in two dimensional for a subcritical flow through a staggered tube bundle. The predicted turbulence statistics (mean and r.m.s velocities) have good agreement with the experimental data (S. Balabani, 1996). Turbulent quantities such as turbulent kinetic energy and dissipation rate are predicted using RANS models and compared with each other. The sensitivity of grid and time-step size have been analyzed. Model constants sensitivity study have been carried out by adopting κ – ε model. It has been observed that model constants are very sensitive to turbulence statistics and turbulent quantities.

Scheduling dynamic dataflow graphs with model checking

With the shift towards many-core computer architectures, dataflow programming has been proposed as one potential solution for producing software that scales to a varying number of processor cores. Programming for parallel architectures is considered difficult as the current popular programming languages are inherently sequential and introducing parallelism is typically up to the programmer. Dataflow, however, is inherently parallel, describing an application as a directed graph, where nodes represent calculations and edges represent a data dependency in form of a queue. These queues are the only allowed communication between the nodes, making the dependencies between the nodes explicit and thereby also the parallelism. Once a node have the su cient inputs available, the node can, independently of any other node, perform calculations, consume inputs, and produce outputs. Data ow models have existed for several decades and have become popular for describing signal processing applications as the graph representation is a very natural representation within this eld. Digital lters are typically described with boxes and arrows also in textbooks. Data ow is also becoming more interesting in other domains, and in principle, any application working on an information stream ts the dataflow paradigm. Such applications are, among others, network protocols, cryptography, and multimedia applications. As an example, the MPEG group standardized a dataflow language called RVC-CAL to be use within reconfigurable video coding. Describing a video coder as a data ow network instead of with conventional programming languages, makes the coder more readable as it describes how the video dataflows through the different coding tools. While dataflow provides an intuitive representation for many applications, it also introduces some new problems that need to be solved in order for data ow to be more widely used. The explicit parallelism of a dataflow program is descriptive and enables an improved utilization of available processing units, however, the independent nodes also implies that some kind of scheduling is required. The need for efficient scheduling becomes even more evident when the number of nodes is larger than the number of processing units and several nodes are running concurrently on one processor core. There exist several data ow models of computation, with different trade-offs between expressiveness and analyzability. These vary from rather restricted but statically schedulable, with minimal scheduling overhead, to dynamic where each ring requires a ring rule to evaluated. The model used in this work, namely RVC-CAL, is a very expressive language, and in the general case it requires dynamic scheduling, however, the strong encapsulation of dataflow nodes enables analysis and the scheduling overhead can be reduced by using quasi-static, or piecewise static, scheduling techniques. The scheduling problem is concerned with nding the few scheduling decisions that must be run-time, while most decisions are pre-calculated. The result is then an, as small as possible, set of static schedules that are dynamically scheduled. To identify these dynamic decisions and to find the concrete schedules, this thesis shows how quasi-static scheduling can be represented as a model checking problem. This involves identifying the relevant information to generate a minimal but complete model to be used for model checking. The model must describe everything that may affect scheduling of the application while omitting everything else in order to avoid state space explosion. This kind of simplification is necessary to make the state space analysis feasible. For the model checker to nd the actual schedules, a set of scheduling strategies are de ned which are able to produce quasi-static schedulers for a wide range of applications. The results of this work show that actor composition with quasi-static scheduling can be used to transform data ow programs to t many different computer architecture with different type and number of cores. This in turn, enables dataflow to provide a more platform independent representation as one application can be fitted to a specific processor architecture without changing the actual program representation. Instead, the program representation is in the context of design space exploration optimized by the development tools to fit the target platform. This work focuses on representing the dataflow scheduling problem as a model checking problem and is implemented as part of a compiler infrastructure. The thesis also presents experimental results as evidence of the usefulness of the approach.