Limit on the pion distribution amplitude

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Scattering and bound states of fermions in the modified Hulthen potential

In the present work the scattering of a fermion in the modified Hulthen potential is considered with a general vector and scalar and we solved the Dirac equation in the one-dimensional space. The transmission and reflection coefficients are reported. The bound-state solution is also given. The study shows the asymptotic behavior of the wave function in bound-state and scattering states solutions.

Oscilador harmônico unidimensional via transformada unilateral de Fourier e mapeamento do potencial de morse unidimensional nos potenciais pseudo-harmônico e de kratzer tridimensionais

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estudo de sistemas não-lineares em mecânica quântica

The main goal of this work is to investigate the effects of a nonlinear cubic term inserted in the Schrödinger equation for one-dimensional potentials studied in Quantum Mechanics textbooks. Being the main tool the numerical analysis in a large number of works, the analysis of this effect by this term in the potential itself, in order to work with an analytical solution, can be considered something new. For the harmonic oscillator potential, the analysis was made from a numerical method, comparing the result with the known results in the literature. In the case of the infinite well potential and the step potential, hoping to work with an analytical solution, by construction we started with the known wavefunction for the linear case noting the effects in the other physical quantities. The coupling of the physical quantities involved in this work has yielded, besides many complications in the calculations, a series of conditions on the existence and validity of the solutions in regard to the system possible configurations

Oscilador harmônico unidimensional via transformada unilateral de Fourier e mapeamento do potencial de morse unidimensional nos potenciais pseudo-harmônico e de kratzer tridimensionais

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estudo de sistemas não-lineares em mecânica quântica

The main goal of this work is to investigate the effects of a nonlinear cubic term inserted in the Schrödinger equation for one-dimensional potentials studied in Quantum Mechanics textbooks. Being the main tool the numerical analysis in a large number of works, the analysis of this effect by this term in the potential itself, in order to work with an analytical solution, can be considered something new. For the harmonic oscillator potential, the analysis was made from a numerical method, comparing the result with the known results in the literature. In the case of the infinite well potential and the step potential, hoping to work with an analytical solution, by construction we started with the known wavefunction for the linear case noting the effects in the other physical quantities. The coupling of the physical quantities involved in this work has yielded, besides many complications in the calculations, a series of conditions on the existence and validity of the solutions in regard to the system possible configurations

One-dimensional point interaction with Griffiths' boundary conditions

Griffiths proposed a pair of boundary conditions that define a point interaction in one dimensional quantum mechanics. The conditions involve the nth derivative of the wave function where n is a non-negative integer. We re-examine the interaction so defined and explicitly confirm that it is self-adjoint for any even value of n and for n = 1. The interaction is not self-adjoint for odd n > 1. We then propose a similar but different pair of boundary conditions with the nth derivative of the wave function such that the ensuing point interaction is self-adjoint for any value of n.

On the Deactivation Mechanisms of Adenine-Thymine Base Pair

In this contribution, the multiconfigurational second-order perturbation theory method based on a complete active space reference wave function (CASSCF/CASPT2) is applied to study all possible single and double proton/hydrogen transfers between the nucleobases in the adenine-thymine (AT) base pair, analyzing the role of excited states with different nature [localized (LE) and charge transfer (CT)] and considering concerted as well as step-wise mechanisms. According to the findings, once the lowest excited states, localized in adenine, are populated during UV irradiation of the Watson-Crick base pair, the proton transfer in the N-O bridge does not require high energy in order to populate a CT state. The latter state will immediately relax toward a crossing with the ground state, which will funnel the system to either the canonical structure or the imino-enol tautomer. The base pair is also capable of repairing itself easily since the imino-enol species is unstable to thermal conversion.

Stopping power and depth dose profile of H+ and He+ ion beams in hydroxyapatite thin films.

Hadron therapy is a promising technique to treat deep-seated tumors. For an accurate treatment planning, the energy deposition in the soft and hard human tissue must be well known. Water has been usually employed as a phantom of soft tissues, but other biomaterials, such as hydroxyapatite (HAp), used as bone substitute, are also relevant as a phantom for hard tissues. The stopping power of HAp for H+ and He+ beams has been studied experimentally and theoretically. The measurements have been done using the Rutherford backscattering technique in an energy range of 450-2000 keV for H+ and of 400-5000 keV for He+ projectiles. The theoretical calculations are based in the dielectric formulation together with the MELF-GOS (Mermin Energy-Loss Function – Generalized Oscillator Strengths) method [1] to describe the target excitation spectrum. A quite good agreement between the experimental data and the theoretical results has been found. The depth dose profile of H+ and He+ ion beams in HAp has been simulated by the SEICS (Simulation of Energetic Ions and Clusters through Solids) code [2], which incorporates the electronic stopping force due to the energy loss by collisions with the target electrons, including fluctuations due to the energy-loss straggling, the multiple elastic scattering with the target nuclei, with their corresponding nuclear energy loss, and the dynamical charge-exchange processes in the projectile charge state. The energy deposition by H+ and He+ as a function of the depth are compared, at several projectile energies, for HAp and liquid water, showing important differences.

Strahlungsbegleitete Pion-Photoproduktion am Proton

Im Jahre 1997 wurden von Tatischeff et al. bei der Reaktion p p -> X p pi+ resonanzartige Zustände im Spektrum der invarianten Masse des fehlenden Nukleons X bei M = 1004, 1044 und 1094 MeV gefunden. In einem zweiten Experiment von Filkov et al. beobachtete man bei der Reaktion p d -> p p X Resonanzstrukturen bei M = 966, 986 und 1003 MeV. Solche exotischen Resonanzen widersprechen etablierten Nukleonenmodellen, die die Delta(1232)-Resonanz als ersten Anregungszustand beschreiben. Zur Deutung der beobachteten Strukturen wurden Quarkcluster-Modelle mit und ohne Farb-Magnet-Wechselwirkungen entwickelt. Lvov et al. zweifelten die experimentellen Ergebnisse an, da keine Strukturen in den Daten zur reellen Comptonstreuung gefunden wurden. Als Gegenargument wurde von Kobushkin vorgeschlagen, dass diese Resonanzen eine total-antisymmetrische Spin-Flavour-Wellenfunktion haben und nur der N-2Gamma-Zerfall erlaubt wäre. In dieser Arbeit wurde die Reaktion g p -> X pi+ -> n g g pi+ zur Suche nach diesen exotischen Resonanzen verwendet. Die Daten wurden parallel zur Messung der Pion-Polarisierbarkeiten am Mainzer Beschleuniger MAMI genommen. Durch Bremsstrahlung der Elektronen an einer Radiatorfolie wurden reelle Photonen erzeugt, deren Energie von der A2-Photonenmarkierungsanlage (Glasgow-Tagger) bestimmt wurde. Als Protontarget wurde ein 10 cm langes Flüssigwasserstoff-Target verwendet. Geladene Reaktionsprodukte wurden unter Vorwärtswinkeln Theta < 20 Grad bezüglich der Strahlachse in einer Vieldraht-Proportionalkammer nachgewiesen, während Photonen im Spektrometer TAPS mit 526 BaF2-Kristallen unter Polarwinkeln Theta > 60 Grad detektiert wurden. Zum Nachweis von Neutronen stand ein Flugzeitdetektor mit insgesamt 111 Einzelmodulen zur Verfügung. Zum Test der Analysesoftware und des experimentellen Aufbaus wurden zusätzlich die Reaktionskanäle g p -> p pi0 und g p -> n pi0 pi+ ausgewertet. Für die Ein-Pion-Produktion wurden differentielle Wirkungsquerschnitte unter Rückwärtswinkeln bestimmt und mit theoretischen Modellen und experimentellen Werten verglichen. Für den Kanal g p -> n pi0 pi+ wurden Spektren invarianter Massen für verschiedene Teilchenkombinationen ermittelt und mit einer Simulation verglichen. Die Daten legen nahe, dass die Reaktion hauptsächlich über eine Anregung der Delta0(1232)-Resonanz verläuft. Bei der Suche nach exotischen Resonanzen wurden keine statistisch signifikanten Strukturen gefunden. Es wurden Obergrenzen für den differentiellen Wirkungsquerschnitt ermittelt.

Renormalization of fermion mixing

Precision measurements of phenomena related to fermion mixing require the inclusion of higher order corrections in the calculation of corresponding theoretical predictions. For this, a complete renormalization scheme for models that allow for fermion mixing is highly required. The correct treatment of unstable particles makes this task difficult and yet, no satisfactory and general solution can be found in the literature. In the present work, we study the renormalization of the fermion Lagrange density with Dirac and Majorana particles in models that involve mixing. The first part of the thesis provides a general renormalization prescription for the Lagrangian, while the second one is an application to specific models. In a general framework, using the on-shell renormalization scheme, we identify the physical mass and the decay width of a fermion from its full propagator. The so-called wave function renormalization constants are determined such that the subtracted propagator is diagonal on-shell. As a consequence of absorptive parts in the self-energy, the constants that are supposed to renormalize the incoming fermion and the outgoing antifermion are different from the ones that should renormalize the outgoing fermion and the incoming antifermion and not related by hermiticity, as desired. Instead of defining field renormalization constants identical to the wave function renormalization ones, we differentiate the two by a set of finite constants. Using the additional freedom offered by this finite difference, we investigate the possibility of defining field renormalization constants related by hermiticity. We show that for Dirac fermions, unless the model has very special features, the hermiticity condition leads to ill-defined matrix elements due to self-energy corrections of external legs. In the case of Majorana fermions, the constraints for the model are less restrictive. Here one might have a better chance to define field renormalization constants related by hermiticity. After analysing the complete renormalized Lagrangian in a general theory including vector and scalar bosons with arbitrary renormalizable interactions, we consider two specific models: quark mixing in the electroweak Standard Model and mixing of Majorana neutrinos in the seesaw mechanism. A counter term for fermion mixing matrices can not be fixed by only taking into account self-energy corrections or fermion field renormalization constants. The presence of unstable particles in the theory can lead to a non-unitary renormalized mixing matrix or to a gauge parameter dependence in its counter term. Therefore, we propose to determine the mixing matrix counter term by fixing the complete correction terms for a physical process to experimental measurements. As an example, we calculate the decay rate of a top quark and of a heavy neutrino. We provide in each of the chosen models sample calculations that can be easily extended to other theories.

Simulations of charge transport in organic compounds

To aid the design of organic semiconductors, we study the charge transport properties of organic liquid crystals, i.e. hexabenzocoronene and carbazole macrocycle, and single crystals, i.e. rubrene, indolocarbazole and benzothiophene derivatives (BTBT, BBBT). The aim is to find structure-property relationships linking the chemical structure as well as the morphology with the bulk charge carrier mobility of the compounds. To this end, molecular dynamics (MD) simulations are performed yielding realistic equilibrated morphologies. Partial charges and molecular orbitals are calculated based on single molecules in vacuum using quantum chemical methods. The molecular orbitals are then mapped onto the molecular positions and orientations, which allows calculation of the transfer integrals between nearest neighbors using the molecular orbital overlap method. Thus we obtain realistic transfer integral distributions and their autocorrelations. In case of organic crystals the differences between two descriptions of charge transport, namely semi-classical dynamics (SCD) in the small polaron limit and kinetic Monte Carlo (KMC) based on Marcus rates, are studied. The liquid crystals are investigated solely in the hopping limit. To simulate the charge dynamics using KMC, the centers of mass of the molecules are mapped onto lattice sites and the transfer integrals are used to compute the hopping rates. In the small polaron limit, where the electronic wave function is spread over a limited number of neighboring molecules, the Schroedinger equation is solved numerically using a semi-classical approach. The results are compared for the different compounds and methods and, where available, with experimental data. The carbazole macrocycles form columnar structures arranged on a hexagonal lattice with side chains facing inwards, so columns can closely approach each other allowing inter-columnar and thus three-dimensional transport. When taking only intra-columnar transport into account, the mobility is orders of magnitude lower than in the three-dimensional case. BTBT is a promising material for solution-processed organic field-effect transistors. We are able to show that, on the time-scales of charge transport, static disorder due to slow side chain motions is the main factor determining the mobility. The resulting broad transfer integral distributions modify the connectivity of the system but sufficiently many fast percolation paths remain for the charges. Rubrene, indolocarbazole and BBBT are examples of crystals without significant static disorder. The high mobility of rubrene is explained by two main features: first, the shifted cofacial alignment of its molecules, and second, the high center of mass vibrational frequency. In comparsion to SCD, only KMC based on Marcus rates is capable of describing neighbors with low coupling and of taking static disorder into account three-dimensionally. Thus it is the method of choice for crystalline systems dominated by static disorder. However, it is inappropriate for the case of strong coupling and underestimates the mobility of well-ordered crystals. SCD, despite its one-dimensionality, is valuable for crystals with strong coupling and little disorder. It also allows correct treatment of dynamical effects, such as intermolecular vibrations of the molecules. Rate equations are incapable of this, because simulations are performed on static snapshots. We have thus shown strengths and weaknesses of two state of the art models used to study charge transport in organic compounds, partially developed a program to compute and visualize transfer integral distributions and other charge transport properties, and found structure-mobility relations for several promising organic semiconductors.

Higher-order perturbative relativistic corrections to energies and properties

Relativistic effects need to be considered in quantum-chemical calculations on systems including heavy elements or when aiming at high accuracy for molecules containing only lighter elements. In the latter case, consideration of relativistic effects via perturbation theory is an attractive option. Among the available techniques, Direct Perturbation Theory (DPT) in its lowest order (DPT2) has become a standard tool for the calculation of relativistic corrections to energies and properties.In this work, the DPT treatment is extended to the next order (DPT4). It is demonstrated that the DPT4 correction can be obtained as a second derivative of the energy with respect to the relativistic perturbation parameter. Accordingly, differentiation of a suitable Lagrangian, thereby taking into account all constraints on the wave function, provides analytic expressions for the fourth-order energy corrections. The latter have been implemented at the Hartree-Fock level and within second-order Møller-Plesset perturbaton theory using standard analytic second-derivative techniques into the CFOUR program package. For closed-shell systems, the DPT4 corrections consist of higher-order scalar-relativistic effects as well as spin-orbit corrections with the latter appearing here for the first time in the DPT series.Relativistic corrections are reported for energies as well as for first-order electrical properties and compared to results from rigorous four-component benchmark calculations in order to judge the accuracy and convergence of the DPT expansion for both the scalar-relativistic as well as the spin-orbit contributions. Additionally, the importance of relativistic effects to the bromine and iodine quadrupole-coupling tensors is investigated in a joint experimental and theoretical study concerning the rotational spectra of CH2BrF, CHBrF2, and CH2FI.

Variational Monte Carlo study of the hydrogen electronic structure at ultra-high pressure

Die causa finalis der vorliegenden Arbeit ist das Verständnis des Phasendiagramms von Wasserstoff bei ultrahohen Drücken, welche von nichtleitendem H2 bis hin zu metallischem H reichen. Da die Voraussetzungen für ultrahohen Druck im Labor schwer zu schaffen sind, bilden Computersimulationen ein wichtiges alternatives Untersuchungsinstrument. Allerdings sind solche Berechnungen eine große Herausforderung. Eines der größten Probleme ist die genaue Auswertung des Born-Oppenheimer Potentials, welches sowohl für die nichtleitende als auch für die metallische Phase geeignet sein muss. Außerdem muss es die starken Korrelationen berücksichtigen, die durch die kovalenten H2 Bindungen und die eventuellen Phasenübergänge hervorgerufen werden. Auf dieses Problem haben unsere Anstrengungen abgezielt. Im Kontext von Variationellem Monte Carlo (VMC) ist die Shadow Wave Function (SWF) eine sehr vielversprechende Option. Aufgrund ihrer Flexibilität sowohl lokalisierte als auch delokalisierte Systeme zu beschreiben sowie ihrer Fähigkeit Korrelationen hoher Ordnung zu berücksichtigen, ist sie ein idealer Kandidat für unsere Zwecke. Unglücklicherweise bringt ihre Formulierung ein Vorzeichenproblem mit sich, was die Anwendbarkeit limitiert. Nichtsdestotrotz ist es möglich diese Schwierigkeit zu umgehen indem man die Knotenstruktur a priori festlegt. Durch diesen Formalismus waren wir in der Lage die Beschreibung der Elektronenstruktur von Wasserstoff signifikant zu verbessern, was eine sehr vielversprechende Perspektive bietet. Während dieser Forschung haben wir also die Natur des Vorzeichenproblems untersucht, das sich auf die SWF auswirkt, und dabei ein tieferes Verständnis seines Ursprungs erlangt. Die vorliegende Arbeit ist in vier Kapitel unterteilt. Das erste Kapitel führt VMC und die SWF mit besonderer Ausrichtung auf fermionische Systeme ein. Kapitel 2 skizziert die Literatur über das Phasendiagramm von Wasserstoff bei ultrahohem Druck. Das dritte Kapitel präsentiert die Implementierungen unseres VMC Programms und die erhaltenen Ergebnisse. Zum Abschluss fasst Kapitel 4 unsere Bestrebungen zur Lösung des zur SWF zugehörigen Vorzeichenproblems zusammen.

Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections

For (H2O)n where n = 1–10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 K.