Approximate analytic expression for the eigenenergies of the anharmonic oscillator V(x)=Ax6+Bx2

An approximate analytic expression for the eigenenergies of the anharmonic oscillator V(x)=Ax6+Bx2 is introduced, starting from particular analytic solutions which are valid when certain relations between the parameters A and B are held. © 1995 The American Physical Society.

Algebraic solution of an anisotropic ring-shaped oscillator

Using an algebraic technique related to the SO (2, 1) group we construct the Green function for the potential ar2 + b(r sin θ)-2 + c(r cos θ)-2 + dr2 sin2θ + er2 cos2θ. The energy spectrum and the normalized wave functions are also obtained. © 1990.

Accounting for uncertainty in ecological analysis: the strengths and limitations of hierarchical statistical modeling

Analyses of ecological data should account for the uncertainty in the process(es) that generated the data. However, accounting for these uncertainties is a difficult task, since ecology is known for its complexity. Measurement and/or process errors are often the only sources of uncertainty modeled when addressing complex ecological problems, yet analyses should also account for uncertainty in sampling design, in model specification, in parameters governing the specified model, and in initial and boundary conditions. Only then can we be confident in the scientific inferences and forecasts made from an analysis. Probability and statistics provide a framework that accounts for multiple sources of uncertainty. Given the complexities of ecological studies, the hierarchical statistical model is an invaluable tool. This approach is not new in ecology, and there are many examples (both Bayesian and non-Bayesian) in the literature illustrating the benefits of this approach. In this article, we provide a baseline for concepts, notation, and methods, from which discussion on hierarchical statistical modeling in ecology can proceed. We have also planted some seeds for discussion and tried to show where the practical difficulties lie. Our thesis is that hierarchical statistical modeling is a powerful way of approaching ecological analysis in the presence of inevitable but quantifiable uncertainties, even if practical issues sometimes require pragmatic compromises.

Experimental Assessment of the Sensitiveness of an Electrochemical Oscillator towards Chemical Perturbations

In this study we address the problem of the response of a (electro)chemical oscillator towards chemical perturbations of different magnitudes. The chemical perturbation was achieved by addition of distinct amounts of trifluoromethanesulfonate (TFMSA), a rather stable and non-specifically adsorbing anion, and the system under investigation was the methanol electro-oxidation reaction under both stationary and oscillatory regimes. Increasing the anion concentration resulted in a decrease in the reaction rates of methanol oxidation and a general decrease in the parameter window where oscillations occurred. Furthermore, the addition of TFMSA was found to decrease the induction period and the total duration of oscillations. The mechanism underlying these observations was derived mathematically and revealed that inhibition in the methanol oxidation through blockage of active sites was found to further accelerate the intrinsic non-stationarity of the unperturbed system. Altogether, the presented results are among the few concerning the experimental assessment of the sensitiveness of an oscillator towards chemical perturbations. The universal nature of the complex chemical oscillator investigated here might be used for reference when studying the dynamics of other less accessible perturbed networks of (bio)chemical reactions.

Interactions between 1-butyl-3-methylimidazolium tetrafluoroborate ionic liquid and gamma-Al2O3 (100) surface calculated by density functional theory

The main aim of this work is to investigate the 1-butyl-3-methylimidazolium tetrafluoroborate ([C4C1Im]+[BF4]-) ionic liquid (IL) adsorption on the gamma-Al2O3 (100) by density functional theory calculations to try to rationalize the adsorption as an electrostatic phenomenon. Optimized geometries and interaction energies of IL one-monolayer on the gamma-Al2O3 were obtained on high surface coverage (one cationanion pair per 94.96 nm2). A study of dispersion force was made to estimate its contribution to the adsorption. Overall, the process is ruled by electrostatic interaction between ions and surface. Adsorption of the anion [BF4]- and cation [C4C1Im]+ was also studied by Bader charge analysis and charge density difference for supported and unsupported situations. It is suggested that the IL ions have their charges maintained with significant anion cloud polarization inward to the acid aluminum sites. (c) 2012 Wiley Periodicals, Inc.

Genotype x location x feeding interactions in Nellore cattle calculated through three-mode principal component analysis

The main objective of this study was to apply three-mode principal component analysis to assess the triple interaction (genotype x location x feeding) on direct genetic value for weight at 205 days of age. We used 60 sires with offspring in three regions of northeastern Brazil (Maranhao, Mata and Agreste, and Reconcavo Baiano) and raised on a pasture regime or with supplementation. There was no interaction between genotype and location, but there was a correlation between genotype and direct effect of feeding. The use of sires should be dictated according to the system of rearing of their offspring.

CdSe/CdTe interface band gaps and band offsets calculated using spin-orbit and self-energy corrections

We performed ab initio calculations of the electronic structures of bulk CdSe and CdTe, and their interface band alignments on the CdSe in-plane lattice parameters. For this, we employed the LDA-1/2 self-energy correction scheme [L.G. Ferreira, M. Marques, L.K. Teles, Phys. Rev. B 78 (2008) 125116] to obtain corrected band gaps and band offsets. Our calculations include the spin-orbit effects for the bulk cases, which have shown to be of importance for the equilibrium systems and are possibly degraded in these strained semiconductors. Therefore, the SO showed reduced importance for the band alignment of this particular system. Moreover, the electronic structure calculated along the transition region across the CdSe/CdTe interface shows an interesting non-monotonic variation of the band gap in the range 0.8-1.8 eV, which may enhance the absorption of light for corresponding frequencies at the interface between these two materials in photovoltaic applications. (C) 2012 Elsevier B.V. All rights reserved.

Period determination in the food-entrainable and methamphetamine-sensitive circadian oscillator(s)

Daily rhythmic processes are coordinated by circadian clocks, which are present in numerous central and peripheral tissues. In mammals, two circadian clocks, the food-entrainable oscillator (FEO) and methamphetamine-sensitive circadian oscillator (MASCO), are "black box" mysteries because their anatomical loci are unknown and their outputs are not expressed under normal physiological conditions. In the current study, the investigation of the timekeeping mechanisms of the FEO and MASCO in mice with disruption of all three paralogs of the canonical clock gene, Period, revealed unique and convergent findings. We found that both the MASCO and FEO in Per1(-/-)/Per2(-/-)/Per3(-/-) mice are circadian oscillators with unusually short (similar to 21 h) periods. These data demonstrate that the canonical Period genes are involved in period determination in the FEO and MASCO, and computational modeling supports the hypothesis that the FEO and MASCO use the same timekeeping mechanism or are the same circadian oscillator. Finally, these studies identify Per1(-/-)/Per2(-/-)/Per3(-/-) mice as a unique tool critical to the search for the elusive anatomical location(s) of the FEO and MASCO.

Dynamical changes from harmonic vibrations of a limited power supply driving a Duffing oscillator

The harmonic oscillations of a Duffing oscillator driven by a limited power supply are investigated as a function of the alternative strength of the rotor. The semi-trivial and non-trivial solutions are derived. We examine the stability of these solutions and then explore the complex behaviors associated with the bifurcations sequences. Interestingly, a 3D diagram provides a global view of the effects of alternate strength on the appearance of chaos and hyperchaos on the system.

Push-out bond strengths of different dental cements used to cement glass fiber posts

Statement of problem: Since the introduction of glass fiber posts, irreversible vertical root fractures have become a rare occurrence; however, adhesive failure has become the primary failure mode. Purpose: The purpose of this study was to evaluate the push-out bond strength of glass fiber posts cemented with different luting agents on 3 segments of the root. Material and methods: Eighty human maxillary canines with similar root lengths were randomly divided into 8 groups (n=10) according to the cement assessed (Rely X luting, Luting and Lining, Ketac Cem, Rely X ARC, Biscem, Duo-link, Rely X U100, and Variolink II). After standardized post space preparation, the root dentin was pretreated for dualpolymerizing resin cements and untreated for the other cements. The mixed luting cement paste was inserted into post spaces with a spiral file and applied to the post surface that was seated into the canal. After 7 days, the teeth were sectioned perpendicular to their long axis into 1-mm-thick sections. The push-out test was performed at a speed of 0.5 mm/min until extrusion of the post occurred. The results were evaluated by 2-way ANOVA and the all pairwise multiple comparison procedures (Tukey test) (?=.05). Results: ANOVA showed that the type of interaction between cement and root location significantly influenced the push-out strength (P<.05). The highest push-out strength results with root location were obtained with Luting and Lining (S3) (19.5 ±4.9 MPa), Ketac Cem (S2) (18.6 ±5.5 MPa), and Luting and Lining (S1) (18.0 ±7.6 MPa). The lowest mean values were recorded with Variolink II (S1) (4.6 ±4.0 MPa), Variolink II (S2) (1.6 ±1.5 MPa), and Rely X ARC (S3) (0.9 ±1.1 MPa). Conclusions: Self-adhesive cements and glass ionomer cements showed significantly higher values compared to dual-polymerizing resin cements. In all root segments, dual-polymerizing resin cements provided significantly lower bond strength. Significant differences among root segments were found only for Duo-link cement.

Autoionizing states and their relevance in electron-ion recombination

Atomic physics plays an important role in determining the evolution stages in a wide range of laboratory and cosmic plasmas. Therefore, the main contribution to our ability to model, infer and control plasma sources is the knowledge of underlying atomic processes. Of particular importance are reliable low temperature dielectronic recombination (DR) rate coefficients. This thesis provides systematically calculated DR rate coefficients of lithium-like beryllium and sodium ions via ∆n = 0 doubly excited resonant states. The calculations are based on complex-scaled relativistic many-body perturbation theory in an all-order formulation within the single- and double-excitation coupled-cluster scheme, including radiative corrections. Comparison of DR resonance parameters (energy levels, autoionization widths, radiative transition probabilities and strengths) between our theoretical predictions and the heavy-ion storage rings experiments (CRYRING-Stockholm and TSRHeidelberg) shows good agreement. The intruder state problem is a principal obstacle for general application of the coupled-cluster formalism on doubly excited states. Thus, we have developed a technique designed to avoid the intruder state problem. It is based on a convenient partitioning of the Hilbert space and reformulation of the conventional set of pairequations. The general aspects of this development are discussed, and the effectiveness of its numerical implementation (within the non-relativistic framework) is selectively illustrated on autoionizing doubly excited states of helium.

Statistical properties of Radio Halos and the re-acceleration model

Galaxy clusters occupy a special position in the cosmic hierarchy as they are the largest bound structures in the Universe. There is now general agreement on a hierarchical picture for the formation of cosmic structures, in which galaxy clusters are supposed to form by accretion of matter and merging between smaller units. During merger events, shocks are driven by the gravity of the dark matter in the diffuse barionic component, which is heated up to the observed temperature. Radio and hard-X ray observations have discovered non-thermal components mixed with the thermal Intra Cluster Medium (ICM) and this is of great importance as it calls for a “revision” of the physics of the ICM. The bulk of present information comes from the radio observations which discovered an increasing number of Mpcsized emissions from the ICM, Radio Halos (at the cluster center) and Radio Relics (at the cluster periphery). These sources are due to synchrotron emission from ultra relativistic electrons diffusing through µG turbulent magnetic fields. Radio Halos are the most spectacular evidence of non-thermal components in the ICM and understanding the origin and evolution of these sources represents one of the most challenging goal of the theory of the ICM. Cluster mergers are the most energetic events in the Universe and a fraction of the energy dissipated during these mergers could be channelled into the amplification of the magnetic fields and into the acceleration of high energy particles via shocks and turbulence driven by these mergers. Present observations of Radio Halos (and possibly of hard X-rays) can be best interpreted in terms of the reacceleration scenario in which MHD turbulence injected during these cluster mergers re-accelerates high energy particles in the ICM. The physics involved in this scenario is very complex and model details are difficult to test, however this model clearly predicts some simple properties of Radio Halos (and resulting IC emission in the hard X-ray band) which are almost independent of the details of the adopted physics. In particular in the re-acceleration scenario MHD turbulence is injected and dissipated during cluster mergers and thus Radio Halos (and also the resulting hard X-ray IC emission) should be transient phenomena (with a typical lifetime <» 1 Gyr) associated with dynamically disturbed clusters. The physics of the re-acceleration scenario should produce an unavoidable cut-off in the spectrum of the re-accelerated electrons, which is due to the balance between turbulent acceleration and radiative losses. The energy at which this cut-off occurs, and thus the maximum frequency at which synchrotron radiation is produced, depends essentially on the efficiency of the acceleration mechanism so that observations at high frequencies are expected to catch only the most efficient phenomena while, in principle, low frequency radio surveys may found these phenomena much common in the Universe. These basic properties should leave an important imprint in the statistical properties of Radio Halos (and of non-thermal phenomena in general) which, however, have not been addressed yet by present modellings. The main focus of this PhD thesis is to calculate, for the first time, the expected statistics of Radio Halos in the context of the re-acceleration scenario. In particular, we shall address the following main questions: • Is it possible to model “self-consistently” the evolution of these sources together with that of the parent clusters? • How the occurrence of Radio Halos is expected to change with cluster mass and to evolve with redshift? How the efficiency to catch Radio Halos in galaxy clusters changes with the observing radio frequency? • How many Radio Halos are expected to form in the Universe? At which redshift is expected the bulk of these sources? • Is it possible to reproduce in the re-acceleration scenario the observed occurrence and number of Radio Halos in the Universe and the observed correlations between thermal and non-thermal properties of galaxy clusters? • Is it possible to constrain the magnetic field intensity and profile in galaxy clusters and the energetic of turbulence in the ICM from the comparison between model expectations and observations? Several astrophysical ingredients are necessary to model the evolution and statistical properties of Radio Halos in the context of re-acceleration model and to address the points given above. For these reason we deserve some space in this PhD thesis to review the important aspects of the physics of the ICM which are of interest to catch our goals. In Chapt. 1 we discuss the physics of galaxy clusters, and in particular, the clusters formation process; in Chapt. 2 we review the main observational properties of non-thermal components in the ICM; and in Chapt. 3 we focus on the physics of magnetic field and of particle acceleration in galaxy clusters. As a relevant application, the theory of Alfv´enic particle acceleration is applied in Chapt. 4 where we report the most important results from calculations we have done in the framework of the re-acceleration scenario. In this Chapter we show that a fraction of the energy of fluid turbulence driven in the ICM by the cluster mergers can be channelled into the injection of Alfv´en waves at small scales and that these waves can efficiently re-accelerate particles and trigger Radio Halos and hard X-ray emission. The main part of this PhD work, the calculation of the statistical properties of Radio Halos and non-thermal phenomena as expected in the context of the re-acceleration model and their comparison with observations, is presented in Chapts.5, 6, 7 and 8. In Chapt.5 we present a first approach to semi-analytical calculations of statistical properties of giant Radio Halos. The main goal of this Chapter is to model cluster formation, the injection of turbulence in the ICM and the resulting particle acceleration process. We adopt the semi–analytic extended Press & Schechter (PS) theory to follow the formation of a large synthetic population of galaxy clusters and assume that during a merger a fraction of the PdV work done by the infalling subclusters in passing through the most massive one is injected in the form of magnetosonic waves. Then the processes of stochastic acceleration of the relativistic electrons by these waves and the properties of the ensuing synchrotron (Radio Halos) and inverse Compton (IC, hard X-ray) emission of merging clusters are computed under the assumption of a constant rms average magnetic field strength in emitting volume. The main finding of these calculations is that giant Radio Halos are naturally expected only in the more massive clusters, and that the expected fraction of clusters with Radio Halos is consistent with the observed one. In Chapt. 6 we extend the previous calculations by including a scaling of the magnetic field strength with cluster mass. The inclusion of this scaling allows us to derive the expected correlations between the synchrotron radio power of Radio Halos and the X-ray properties (T, LX) and mass of the hosting clusters. For the first time, we show that these correlations, calculated in the context of the re-acceleration model, are consistent with the observed ones for typical µG strengths of the average B intensity in massive clusters. The calculations presented in this Chapter allow us to derive the evolution of the probability to form Radio Halos as a function of the cluster mass and redshift. The most relevant finding presented in this Chapter is that the luminosity functions of giant Radio Halos at 1.4 GHz are expected to peak around a radio power » 1024 W/Hz and to flatten (or cut-off) at lower radio powers because of the decrease of the electron re-acceleration efficiency in smaller galaxy clusters. In Chapt. 6 we also derive the expected number counts of Radio Halos and compare them with available observations: we claim that » 100 Radio Halos in the Universe can be observed at 1.4 GHz with deep surveys, while more than 1000 Radio Halos are expected to be discovered in the next future by LOFAR at 150 MHz. This is the first (and so far unique) model expectation for the number counts of Radio Halos at lower frequency and allows to design future radio surveys. Based on the results of Chapt. 6, in Chapt.7 we present a work in progress on a “revision” of the occurrence of Radio Halos. We combine past results from the NVSS radio survey (z » 0.05 − 0.2) with our ongoing GMRT Radio Halos Pointed Observations of 50 X-ray luminous galaxy clusters (at z » 0.2−0.4) and discuss the possibility to test our model expectations with the number counts of Radio Halos at z » 0.05 − 0.4. The most relevant limitation in the calculations presented in Chapt. 5 and 6 is the assumption of an “averaged” size of Radio Halos independently of their radio luminosity and of the mass of the parent clusters. This assumption cannot be released in the context of the PS formalism used to describe the formation process of clusters, while a more detailed analysis of the physics of cluster mergers and of the injection process of turbulence in the ICM would require an approach based on numerical (possible MHD) simulations of a very large volume of the Universe which is however well beyond the aim of this PhD thesis. On the other hand, in Chapt.8 we report our discovery of novel correlations between the size (RH) of Radio Halos and their radio power and between RH and the cluster mass within the Radio Halo region, MH. In particular this last “geometrical” MH − RH correlation allows us to “observationally” overcome the limitation of the “average” size of Radio Halos. Thus in this Chapter, by making use of this “geometrical” correlation and of a simplified form of the re-acceleration model based on the results of Chapt. 5 and 6 we are able to discuss expected correlations between the synchrotron power and the thermal cluster quantities relative to the radio emitting region. This is a new powerful tool of investigation and we show that all the observed correlations (PR − RH, PR − MH, PR − T, PR − LX, . . . ) now become well understood in the context of the re-acceleration model. In addition, we find that observationally the size of Radio Halos scales non-linearly with the virial radius of the parent cluster, and this immediately means that the fraction of the cluster volume which is radio emitting increases with cluster mass and thus that the non-thermal component in clusters is not self-similar.

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.

Computer simulation of rod-coil block copolymers and tetrapod,polymer mixtures

In this thesis I present a new coarse-grained model suitable to investigate the phase behavior of rod-coil block copolymers on mesoscopic length scales. In this model the rods are represented by hard spherocylinders, whereas the coil block consists of interconnected beads. The interactions between the constituents are based on local densities. This facilitates an efficient Monte-Carlo sampling of the phase space. I verify the applicability of the model and the simulation approach by means of several examples. I treat pure rod systems and mixtures of rod and coil polymers. Then I append coils to the rods and investigate the role of the different model parameters. Furthermore, I compare different implementations of the model. I prove the capability of the rod-coil block copolymers in our model to exhibit typical micro-phase separated configurations as well as extraordinary phases, such as the wavy lamellar state, percolating structuresrnand clusters. Additionally, I demonstrate the metastability of the observed zigzag phase in our model. A central point of this thesis is the examination of the phase behavior of the rod-coil block copolymers in dependence of different chain lengths and interaction strengths between rods and coil. The observations of these studies are summarized in a phase diagram for rod-coil block copolymers. Furthermore, I validate a stabilization of the smectic phase with increasing coil fraction.rnIn the second part of this work I present a side project in which I derive a model permitting the simulation of tetrapods with and without grafted semiconducting block copolymers. The effect of these polymers is added in an implicit manner by effective interactions between the tetrapods. While the depletion interaction is described in an approximate manner within the Asakura-Oosawa model, the free energy penalty for the brush compression is calculated within the Alexander-de Gennes model. Recent experiments with CdSe tetrapods show that grafted tetrapods are clearly much better dispersed in the polymer matrix than bare tetrapods. My simulations confirm that bare tetrapods tend to aggregate in the matrix of excess polymers, while clustering is significantly reduced after grafting polymer chains to the tetrapods. Finally, I propose a possible extension enabling the simulation of a system with fluctuating volume and demonstrate its basic functionality. This study is originated in a cooperation with an experimental group with the goal to analyze the morphology of these systems in order to find the ideal morphology for hybrid solar cells.

Feasibility of texture analysis for the assessment of biochemical changes in meniscal tissue on T1 maps calculated from delayed gadolinium-enhanced magnetic resonance imaging of cartilage data: comparison with conventional relaxation time measurements

To (1) establish the feasibility of texture analysis for the in vivo assessment of biochemical changes in meniscal tissue on delayed gadolinium-enhanced magnetic resonance imaging of cartilage (dGEMRIC), and (2) compare textural with conventional T1 relaxation time measurements calculated from dGEMRIC data ("T1(Gd) relaxation times").