Matching stages of heavy-ion collision models

Heavy-ion reactions and other collective dynamical processes are frequently described by different theoretical approaches for the different stages of the process, like initial equilibration stage, intermediate locally equilibrated fluid dynamical stage, and final freeze-out stage. For the last stage, the best known is the Cooper-Frye description used to generate the phase space distribution of emitted, noninteracting particles from a fluid dynamical expansion or explosion, assuming a final ideal gas distribution, or (less frequently) an out-of-equilibrium distribution. In this work we do not want to replace the Cooper-Frye description, but rather clarify the ways of using it and how to choose the parameters of the distribution and, eventually, how to choose the form of the phase space distribution used in the Cooper-Frye formula. Moreover, the Cooper-Frye formula is used in connection with the freeze-out problem, while the discussion of transition between different stages of the collision is applicable to other transitions also. More recently, hadronization and molecular dynamics models have been matched to the end of a fluid dynamical stage to describe hadronization and freeze-out. The stages of the model description can be matched to each other on space-time hypersurfaces (just like through the frequently used freeze-out hypersurface). This work presents a generalized description of how to match the stages of the description of a reaction to each other, extending the methodology used at freeze-out, in simple covariant form which is easily applicable in its simplest version for most applications.

Influence of the pseudopotential on the properties of liquid rubidium at 315 K

The influence of the pseudopotential on both the structure and the self-diffusion of liquid rubidium at the melting point has been investigated by means of molecular-dynamics calculations. The model potential considered has been computed from the pseudopotential of Ashcroft, the dielectric function of Geldart and Vosko, and a Born-Mayer term. Four different values for the core radius which enters as input in the pseudopotential have been considered. In this way we have been able to observe and interpret the effect of this contribution on the properties of the liquid.

Oxygen adsorption on clean and oxygen precovered Cu(100)

Kuparipinnan hapettumisen alkuvaiheet ovat vielä nykyisin tutkijoille epäselviä. Kuitenkin, jotta hapettumisprosessia voitaisiin säädellä, on sangen tärkeää ymmärtää mistä varsinainen hapettuminen lähtee liikkeelle ja mitkä ovat hapettumisen seuraavat vaiheet. Tähän kysymykseen haetaan vastauksia tässä työssä käyttäen puhtaasti teoreettisia menetelmiä pinnan käsittelyssä. Aikaisempien teoreettisten ja kokeellisten tutkimusten välillä on pieni ristiriita liittyen hapen tarttumistodennäköisyyteen. Teoreettisten tutkimusten mukaan happi ei puhtaalle pinnalle tullessaan näe potentiaalivallia, mutta kokeelliset tutkimukset osoittavat sellaisen kuitenkin olevan. Tuohon ristiriitaan pureudutaan käyttäen aikaisemmista laskuista poikkeavaa kvanttimekaaniseen molekyylidynamiikkaan perustuvaa lähestymistapaa. Työssä havaitaan, että aikaisemmin yleisesti käytetty menetelmä hukkaa huomattavan määrän tietoa ja siten tutkijat eivät voi ainoastaan tyytyä tarkastelemaan kyseisellä menetelmällä saatuja tuloksia. Kuparipinnalle havaittiin, että korkeilla molekyylin kineettisen energian arvolla aikaisemmin suoritetut laskut hajottavista trajektoreista pitävät paikkansa, mutta matalilla kineettisen energian arvoilla molekyyli kohtaa erittäin voimakkaan ``steering'' vaikutuksen ja trajektorit joiden piti olla hajottavia johtavatkin molekulaariseen adsorptioon. Kun hapen konsentraatio pinnalla on suurempi kuin 0.5 ML, pinta rekonstruoituu. Myös rekonstruktion jälkeistä pintaa on tutkittu samanlaisilla menetelmillä kuin puhdasta pintaa. Rekonstruoituneelle pinnalle ei löydetty hajottavia trajektoreita ja havaittiin, että hapelle annetun kineettisen energian matalilla arvoilla myös tässä tapauksessa on erittäin voimakas ``steering'' vaikutus.

The structure of human tRNALys3 anticodon bound to HIV genome is stabilized by modified nucleosides and adjacent mismatch base pairs

Replication of human immunodeficiency virus (HIV) requires base pairing of the reverse transcriptase primer, human tRNA(Lys3), to the viral RNA. Although the major complementary base pairing occurs between the HIV primer binding sequence (PBS) and the tRNA's 3'-terminus, an important discriminatory, secondary contact occurs between the viral A-rich Loop I, 5'-adjacent to the PBS, and the modified, U-rich anticodon domain of tRNA(Lys3). The importance of individual and combined anticodon modifications to the tRNA/HIV-1 Loop I RNA's interaction was determined. The thermal stabilities of variously modified tRNA anticodon region sequences bound to the Loop I of viral sub(sero)types G and B were analyzed and the structure of one duplex containing two modified nucleosides was determined using NMR spectroscopy and restrained molecular dynamics. The modifications 2-thiouridine, s(2)U(34), and pseudouridine, Psi(39), appreciably stabilized the interaction of the anticodon region with the viral subtype G and B RNAs. The structure of the duplex results in two coaxially stacked A-form RNA stems separated by two mismatched base pairs, U(162)*Psi(39) and G(163)*A(38), that maintained a reasonable A-form helix diameter. The tRNA's s(2)U(34) stabilized the interaction between the A-rich HIV Loop I sequence and the U-rich anticodon, whereas the tRNA's Psi(39) stabilized the adjacent mismatched pairs.

Structure-Based, Rational Design of T Cell Receptors.

Adoptive cell transfer using engineered T cells is emerging as a promising treatment for metastatic melanoma. Such an approach allows one to introduce T cell receptor (TCR) modifications that, while maintaining the specificity for the targeted antigen, can enhance the binding and kinetic parameters for the interaction with peptides (p) bound to major histocompatibility complexes (MHC). Using the well-characterized 2C TCR/SIYR/H-2K(b) structure as a model system, we demonstrated that a binding free energy decomposition based on the MM-GBSA approach provides a detailed and reliable description of the TCR/pMHC interactions at the structural and thermodynamic levels. Starting from this result, we developed a new structure-based approach, to rationally design new TCR sequences, and applied it to the BC1 TCR targeting the HLA-A2 restricted NY-ESO-1157-165 cancer-testis epitope. Fifty-four percent of the designed sequence replacements exhibited improved pMHC binding as compared to the native TCR, with up to 150-fold increase in affinity, while preserving specificity. Genetically engineered CD8(+) T cells expressing these modified TCRs showed an improved functional activity compared to those expressing BC1 TCR. We measured maximum levels of activities for TCRs within the upper limit of natural affinity, K D = ∼1 - 5 μM. Beyond the affinity threshold at K D < 1 μM we observed an attenuation in cellular function, in line with the "half-life" model of T cell activation. Our computer-aided protein-engineering approach requires the 3D-structure of the TCR-pMHC complex of interest, which can be obtained from X-ray crystallography. We have also developed a homology modeling-based approach, TCRep 3D, to obtain accurate structural models of any TCR-pMHC complexes when experimental data is not available. Since the accuracy of the models depends on the prediction of the TCR orientation over pMHC, we have complemented the approach with a simplified rigid method to predict this orientation and successfully assessed it using all non-redundant TCR-pMHC crystal structures available. These methods potentially extend the use of our TCR engineering method to entire TCR repertoires for which no X-ray structure is available. We have also performed a steered molecular dynamics study of the unbinding of the TCR-pMHC complex to get a better understanding of how TCRs interact with pMHCs. This entire rational TCR design pipeline is now being used to produce rationally optimized TCRs for adoptive cell therapies of stage IV melanoma.

Técnicas para análise de desempenho de computadores

A series of open source benchmarks for computer performance analysis of personal computers with a focus on computational chemistry calculations is presented. The results returned by these tests are discussed and used to correlate with the actual performance of a set of computers available for research on two computing intensive fields of chemistry, quantum chemical and molecular simulation calculations.

Técnicas de análise do perfil de execução e otimização de programas em química computacional

In this paper we review the basic techniques of performance analysis within the UNIX environment that are relevant in computational chemistry, with particular emphasis on the execution profile using the gprof tool. Two case studies (in ab initio and molecular dynamics calculations) are presented in order to illustrate how execution profiling can be used to effectively identify bottlenecks and to guide source code optimization. Using these profiling and optimization techniques it was possible to obtain significant speedups (of up to 30%) in both cases.

Tetrahydrobenzo[h][1,6]naphthyridine-6-chlorotacrine hybrids as a new family of anti-Alzheimer agents targeting beta-amyloid, tau, and cholinesterase pathologies

Optimization of an essentially inactive 3,4-dihydro-2H-pyrano[3,2-c]quinoline carboxylic ester derivative as acetylcholinesterase (AChE) peripheral anionic site (PAS)-binding motif by double O → NH bioisosteric replacement, combined with molecular hybridization with the AChE catalytic anionic site (CAS) inhibitor 6-chlorotacrine and molecular dynamics-driven optimization of the length of the linker has resulted in the development of the trimethylene-linked 1,2,3,4-tetrahydrobenzo[h][1,6]naphthyridine6-chlorotacrine hybrid 5a as a picomolar inhibitor of human AChE (hAChE). The tetra-, penta-, and octamethylene-linked homologues 5bd have been also synthesized for comparison purposes, and found to retain the nanomolar hAChE inhibitory potency of the parent 6-chlorotacrine. Further biological profiling of hybrids 5ad has shown that they are also potent inhibitors of human butyrylcholinesterase and moderately potent Aβ42 and tau anti-aggregating agents, with IC50 values in the submicromolar and low micromolar range, respectively. Also, in vitro studies using an artificial membrane model have predicted a good brain permeability for hybrids 5ad, and hence, their ability to reach their targets in the central nervous system. The multitarget profile of the novel hybrids makes them promising leads for developing anti-Alzheimer drug candidates with more balanced biological activities.

Tetrahydrobenzo[h][1,6]naphthyridine-6-chlorotacrine hybrids as a new family of anti-Alzheimer agents targeting beta-amyloid, tau, and cholinesterase pathologies

Optimization of an essentially inactive 3,4-dihydro-2H-pyrano[3,2-c]quinoline carboxylic ester derivative as acetylcholinesterase (AChE) peripheral anionic site (PAS)-binding motif by double O → NH bioisosteric replacement, combined with molecular hybridization with the AChE catalytic anionic site (CAS) inhibitor 6-chlorotacrine and molecular dynamics-driven optimization of the length of the linker has resulted in the development of the trimethylene-linked 1,2,3,4-tetrahydrobenzo[h][1,6]naphthyridine6-chlorotacrine hybrid 5a as a picomolar inhibitor of human AChE (hAChE). The tetra-, penta-, and octamethylene-linked homologues 5bd have been also synthesized for comparison purposes, and found to retain the nanomolar hAChE inhibitory potency of the parent 6-chlorotacrine. Further biological profiling of hybrids 5ad has shown that they are also potent inhibitors of human butyrylcholinesterase and moderately potent Aβ42 and tau anti-aggregating agents, with IC50 values in the submicromolar and low micromolar range, respectively. Also, in vitro studies using an artificial membrane model have predicted a good brain permeability for hybrids 5ad, and hence, their ability to reach their targets in the central nervous system. The multitarget profile of the novel hybrids makes them promising leads for developing anti-Alzheimer drug candidates with more balanced biological activities.

Tetrahydrobenzo[h][1,6]naphthyridine-6-chlorotacrine hybrids as a new family of anti-Alzheimer agents targeting beta-amyloid, tau, and cholinesterase pathologies

Optimization of an essentially inactive 3,4-dihydro-2H-pyrano[3,2-c]quinoline carboxylic ester derivative as acetylcholinesterase (AChE) peripheral anionic site (PAS)-binding motif by double O → NH bioisosteric replacement, combined with molecular hybridization with the AChE catalytic anionic site (CAS) inhibitor 6-chlorotacrine and molecular dynamics-driven optimization of the length of the linker has resulted in the development of the trimethylene-linked 1,2,3,4-tetrahydrobenzo[h][1,6]naphthyridine6-chlorotacrine hybrid 5a as a picomolar inhibitor of human AChE (hAChE). The tetra-, penta-, and octamethylene-linked homologues 5bd have been also synthesized for comparison purposes, and found to retain the nanomolar hAChE inhibitory potency of the parent 6-chlorotacrine. Further biological profiling of hybrids 5ad has shown that they are also potent inhibitors of human butyrylcholinesterase and moderately potent Aβ42 and tau anti-aggregating agents, with IC50 values in the submicromolar and low micromolar range, respectively. Also, in vitro studies using an artificial membrane model have predicted a good brain permeability for hybrids 5ad, and hence, their ability to reach their targets in the central nervous system. The multitarget profile of the novel hybrids makes them promising leads for developing anti-Alzheimer drug candidates with more balanced biological activities.

Tetrahydrobenzo[h][1,6]naphthyridine-6-chlorotacrine hybrids as a new family of anti-Alzheimer agents targeting beta-amyloid, tau, and cholinesterase pathologies

Optimization of an essentially inactive 3,4-dihydro-2H-pyrano[3,2-c]quinoline carboxylic ester derivative as acetylcholinesterase (AChE) peripheral anionic site (PAS)-binding motif by double O → NH bioisosteric replacement, combined with molecular hybridization with the AChE catalytic anionic site (CAS) inhibitor 6-chlorotacrine and molecular dynamics-driven optimization of the length of the linker has resulted in the development of the trimethylene-linked 1,2,3,4-tetrahydrobenzo[h][1,6]naphthyridine6-chlorotacrine hybrid 5a as a picomolar inhibitor of human AChE (hAChE). The tetra-, penta-, and octamethylene-linked homologues 5bd have been also synthesized for comparison purposes, and found to retain the nanomolar hAChE inhibitory potency of the parent 6-chlorotacrine. Further biological profiling of hybrids 5ad has shown that they are also potent inhibitors of human butyrylcholinesterase and moderately potent Aβ42 and tau anti-aggregating agents, with IC50 values in the submicromolar and low micromolar range, respectively. Also, in vitro studies using an artificial membrane model have predicted a good brain permeability for hybrids 5ad, and hence, their ability to reach their targets in the central nervous system. The multitarget profile of the novel hybrids makes them promising leads for developing anti-Alzheimer drug candidates with more balanced biological activities.

Control and femtosecond time-resolved imaging of torsion in a chiral molecule

We study how the combination of long and short laser pulses can be used to induce torsion in an axially chiral biphenyl derivative (3,5-difluoro-3 ,5 -dibromo-4 -cyanobiphenyl). A long, with respect to the molecular rotational periods, elliptically polarized laser pulse produces 3D alignment of the molecules, and a linearly polarized short pulse initiates torsion about the stereogenic axis. The torsional motion is monitored in real-time by measuring the dihedral angle using femtosecond time-resolved Coulomb explosion imaging. Within the first 4 picoseconds (ps), torsion occurs with a period of 1.25 ps and an amplitude of 3◦ in excellent agreement with theoretical calculations. At larger times, the quantum states of the molecules describing the torsional motion dephase and an almost isotropic distribution of the dihedral angle is measured.We demonstrate an original application of covariance analysis of two-dimensional ion images to reveal strong correlations between specific ejected ionic fragments from Coulomb explosion. This technique strengthens our interpretation of the experimental data

Um paradigma da química medicinal: a flexibilidade dos ligantes e receptores

In general, molecular modeling techniques applied in medicinal chemistry have been static and drug based. However the active site geometry and the intrinsic flexibility of both receptor and ligand are fundamental properties for molecular recognition and drug action. As a consequence, the use of dynamic models to describe the ligand-receptor complex is becoming a more common procedure. In this work we discuss the relevance of considering the receptor structure in medicinal chemistry studies as well as the flexibility of the ligand-receptor complex.

Computational and Experimental Investigation of the Grooved Roll in Paper Machine Environment

In the paper machine, it is not a desired feature for the boundary layer flows in the fabric and the roll surfaces to travel into the closing nips, creating overpressure. In this thesis, the aerodynamic behavior of the grooved roll and smooth rolls is compared in order to understand the nip flow phenomena, which is the main reason why vacuum and grooved roll constructions are designed. A common method to remove the boundary layer flow from the closing nip is to use the vacuum roll construction. The downside of the use of vacuum rolls is high operational costs due to pressure losses in the vacuum roll shell. The deep grooved roll has the same goal, to create a pressure difference over the paper web and keep the paper attached to the roll or fabric surface in the drying pocket of the paper machine. A literature review revealed that the aerodynamic functionality of the grooved roll is not very well known. In this thesis, the aerodynamic functionality of the grooved roll in interaction with a permeable or impermeable wall is studied by varying the groove properties. Computational fluid dynamics simulations are utilized as the research tool. The simulations have been performed with commercial fluid dynamics software, ANSYS Fluent. Simulation results made with 3- and 2-dimensional fluid dynamics models are compared to laboratory scale measurements. The measurements have been made with a grooved roll simulator designed for the research. The variables in the comparison are the paper or fabric wrap angle, surface velocities, groove geometry and wall permeability. Present-day computational and modeling resources limit grooved roll fluid dynamics simulations in the paper machine scale. Based on the analysis of the aerodynamic functionality of the grooved roll, a grooved roll simulation tool is proposed. The smooth roll simulations show that the closing nip pressure does not depend on the length of boundary layer development. The surface velocity increase affects the pressure distribution in the closing and opening nips. The 3D grooved roll model reveals the aerodynamic functionality of the grooved roll. With the optimal groove size it is possible to avoid closing nip overpressure and keep the web attached to the fabric surface in the area of the wrap angle. The groove flow friction and minor losses play a different role when the wrap angle is changed. The proposed 2D grooved roll simulation tool is able to replicate the grooved aerodynamic behavior with reasonable accuracy. A small wrap angle predicts the pressure distribution correctly with the chosen approach for calculating the groove friction losses. With a large wrap angle, the groove friction loss shows too large pressure gradients, and the way of calculating the air flow friction losses in the groove has to be reconsidered. The aerodynamic functionality of the grooved roll is based on minor and viscous losses in the closing and opening nips as well as in the grooves. The proposed 2D grooved roll model is a simplification in order to reduce computational and modeling efforts. The simulation tool makes it possible to simulate complex paper machine constructions in the paper machine scale. In order to use the grooved roll as a replacement for the vacuum roll, the grooved roll properties have to be considered on the basis of the web handling application.

Novel Robot Solutions for Carrying out Field Joint Welding and Machining in the Assembly of the Vacuum Vessel of ITER

It is necessary to use highly specialized robots in ITER (International Thermonuclear Experimental Reactor) both in the manufacturing and maintenance of the reactor due to a demanding environment. The sectors of the ITER vacuum vessel (VV) require more stringent tolerances than normally expected for the size of the structure involved. VV consists of nine sectors that are to be welded together. The vacuum vessel has a toroidal chamber structure. The task of the designed robot is to carry the welding apparatus along a path with a stringent tolerance during the assembly operation. In addition to the initial vacuum vessel assembly, after a limited running period, sectors need to be replaced for repair. Mechanisms with closed-loop kinematic chains are used in the design of robots in this work. One version is a purely parallel manipulator and another is a hybrid manipulator where the parallel and serial structures are combined. Traditional industrial robots that generally have the links actuated in series are inherently not very rigid and have poor dynamic performance in high speed and high dynamic loading conditions. Compared with open chain manipulators, parallel manipulators have high stiffness, high accuracy and a high force/torque capacity in a reduced workspace. Parallel manipulators have a mechanical architecture where all of the links are connected to the base and to the end-effector of the robot. The purpose of this thesis is to develop special parallel robots for the assembly, machining and repairing of the VV of the ITER. The process of the assembly and machining of the vacuum vessel needs a special robot. By studying the structure of the vacuum vessel, two novel parallel robots were designed and built; they have six and ten degrees of freedom driven by hydraulic cylinders and electrical servo motors. Kinematic models for the proposed robots were defined and two prototypes built. Experiments for machine cutting and laser welding with the 6-DOF robot were carried out. It was demonstrated that the parallel robots are capable of holding all necessary machining tools and welding end-effectors in all positions accurately and stably inside the vacuum vessel sector. The kinematic models appeared to be complex especially in the case of the 10-DOF robot because of its redundant structure. Multibody dynamics simulations were carried out, ensuring sufficient stiffness during the robot motion. The entire design and testing processes of the robots appeared to be complex tasks due to the high specialization of the manufacturing technology needed in the ITER reactor, while the results demonstrate the applicability of the proposed solutions quite well. The results offer not only devices but also a methodology for the assembly and repair of ITER by means of parallel robots.