Quantum noise in optical fibers. I. Stochastic equations

We analyze the quantum dynamics of radiation propagating in a single-mode optical fiber with dispersion, nonlinearity, and Raman coupling to thermal phonons. We start from a fundamental Hamiltonian that includes the principal known nonlinear effects and quantum-noise sources, including linear gain and loss. Both Markovian and frequency-dependent, non-Markovian reservoirs are treated. This treatment allows quantum Langevin equations, which have a classical form except for additional quantum-noise terms, to be calculated. In practical calculations, it is more useful to transform to Wigner or 1P quasi-probability operator representations. These transformations result in stochastic equations that can be analyzed by use of perturbation theory or exact numerical techniques. The results have applications to fiber-optics communications, networking, and sensor technology.

Shift of the subharmonic resonances and suppression of fluorescence in a two-level atom driven by a bichromatic field

We study the behavior of a two-level atom that is driven by a bichromatic field consisting of a strong resonant component and a weaker tunable component. In addition to the splitting of the energy levels (the multiphoton AC Stark effect), we find that the weaker component also shifts the subharmonic resonances, an effect we attribute to a dynamic Stark shift. When the weaker component is tuned to a shifted resonance, no fluorescence occurs at either the frequency of the strong component or the three-photon mixing frequency. Results are obtained with numerical techniques and explained in terms of the dressed-atom model of the system. (C) 1998 Optical Society of America [S0740-3224(98)01508-2] OCIS codes: 270.4180, 270.6620, 270.0270.

Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

Maximum Likelihood Estimation of Mixture Densities for Binned and Truncated Multivariate Data

Binning and truncation of data are common in data analysis and machine learning. This paper addresses the problem of fitting mixture densities to multivariate binned and truncated data. The EM approach proposed by McLachlan and Jones (Biometrics, 44: 2, 571-578, 1988) for the univariate case is generalized to multivariate measurements. The multivariate solution requires the evaluation of multidimensional integrals over each bin at each iteration of the EM procedure. Naive implementation of the procedure can lead to computationally inefficient results. To reduce the computational cost a number of straightforward numerical techniques are proposed. Results on simulated data indicate that the proposed methods can achieve significant computational gains with no loss in the accuracy of the final parameter estimates. Furthermore, experimental results suggest that with a sufficient number of bins and data points it is possible to estimate the true underlying density almost as well as if the data were not binned. The paper concludes with a brief description of an application of this approach to diagnosis of iron deficiency anemia, in the context of binned and truncated bivariate measurements of volume and hemoglobin concentration from an individual's red blood cells.

Computation of the linear elastic properties of random porous materials with a wide variety of microstructure

A finite-element method is used to study the elastic properties of random three-dimensional porous materials with highly interconnected pores. We show that Young's modulus, E, is practically independent of Poisson's ratio of the solid phase, nu(s), over the entire solid fraction range, and Poisson's ratio, nu, becomes independent of nu(s) as the percolation threshold is approached. We represent this behaviour of nu in a flow diagram. This interesting but approximate behaviour is very similar to the exactly known behaviour in two-dimensional porous materials. In addition, the behaviour of nu versus nu(s) appears to imply that information in the dilute porosity limit can affect behaviour in the percolation threshold limit. We summarize the finite-element results in terms of simple structure-property relations, instead of tables of data, to make it easier to apply the computational results. Without using accurate numerical computations, one is limited to various effective medium theories and rigorous approximations like bounds and expansions. The accuracy of these equations is unknown for general porous media. To verify a particular theory it is important to check that it predicts both isotropic elastic moduli, i.e. prediction of Young's modulus alone is necessary but not sufficient. The subtleties of Poisson's ratio behaviour actually provide a very effective method for showing differences between the theories and demonstrating their ranges of validity. We find that for moderate- to high-porosity materials, none of the analytical theories is accurate and, at present, numerical techniques must be relied upon.

Caracterização dinâmica de um edifício existente de betão armado através de ensaios de vibração ambiental

Dissertação de natureza Científica para obtenção do grau de Mestre em Engenharia Civil

Comparação de técnicas analíticas e numéricas para previsão de resistência de juntas adesivas de sobreposição simples

As juntas adesivas têm vindo a ser usadas em diversas áreas e contam com inúmeras aplicações práticas. Devido ao fácil e rápido fabrico, as juntas de sobreposição simples (JSS) são um tipo de configuração bastante comum. O aumento da resistência, a redução de peso e a resistência à corrosão são algumas das vantagens que este tipo de junta oferece relativamente aos processos de ligação tradicionais. Contudo, a concentração de tensões nas extremidades do comprimento da ligação é uma das principais desvantagens. Existem poucas técnicas de dimensionamento precisas para a diversidade de ligações que podem ser encontradas em situações reais, o que constitui um obstáculo à utilização de juntas adesivas em aplicações estruturais. O presente trabalho visa comparar diferentes métodos analíticos e numéricos na previsão da resistência de JSS com diferentes comprimentos de sobreposição (LO). O objectivo fundamental é avaliar qual o melhor método para prever a resistência das JSS. Foram produzidas juntas adesivas entre substratos de alumínio utilizando um adesivo époxido frágil (Araldite® AV138), um adesivo epóxido moderadamente dúctil (Araldite® 2015), e um adesivo poliuretano dúctil (SikaForce® 7888). Consideraram-se diferentes métodos analíticos e dois métodos numéricos: os Modelos de Dano Coesivo (MDC) e o Método de Elementos Finitos Extendido (MEFE), permitindo a análise comparativa. O estudo possibilitou uma percepção crítica das capacidades de cada método consoante as características do adesivo utilizado. Os métodos analíticos funcionam apenas relativamente bem em condições muito específicas. A análise por MDC com lei triangular revelou ser um método bastante preciso, com excepção de adesivos que sejam bastante dúcteis. Por outro lado, a análise por MEFE demonstrou ser uma técnica pouco adequada, especialmente para o crescimento de dano em modo misto.

Design of Novel Frequency Generation Circuit with Application to Ultra-Wideband Distributed Oscillators

Given the urgence of a new paradigm in wireless digital trasmission which should allow for higher bit rate, lower latency and tigher delay constaints, it has been proposed to investigate the fundamental building blocks that at the circuital/device level, will boost the change towards a more efficient network architecture, with high capacity, higher bandwidth and a more satisfactory end user experience. At the core of each transciever, there are inherently analog devices capable of providing the carrier signal, the oscillators. It is strongly believed that many limitations in today's communication protocols, could be relieved by permitting high carrier frequency radio transmission, and having some degree of reconfigurability. This led us to studying distributed oscillator architectures which work in the microwave range and possess wideband tuning capability. As microvave oscillators are essentially nonlinear devices, a full nonlinear analyis, synthesis, and optimization had to be considered for their implementation. Consequently, all the most used nonlinear numerical techniques in commercial EDA software had been reviewed. An application of all the aforementioned techniques has been shown, considering a systems of three coupled oscillator ("triple push" oscillator) in which the stability of the various oscillating modes has been studied. Provided that a certain phase distribution is maintained among the oscillating elements, this topology permits a rise in the output power of the third harmonic; nevertheless due to circuit simmetry, "unwanted" oscillating modes coexist with the intenteded one. Starting with the necessary background on distributed amplification and distributed oscillator theory, the design of a four stage reverse mode distributed voltage controlled oscillator (DVCO) using lumped elments has been presented. All the design steps have been reported and for the first time a method for an optimized design with reduced variations in the output power has been presented. Ongoing work is devoted to model a wideband DVCO and to implement a frequency divider.

Can extreme black Holes have (long) Abelian Higgs hair?

It has been argued that a black hole horizon can support the long-range fields of a Nielsen-Olesen string and that one can think of such a vortex as black hole "hair." In this paper, we examine the properties of an Abelian Higgs vortex in the presence of a charged black hole as we allow the hole to approach extremality. Using both analytical and numerical techniques, we show that the magnetic field lines (as well as the scalar field) of the vortex are completely expelled from the black hole in the extreme limit. This was to be expected, since extreme black holes in Einstein-Maxwell theory are known to exhibit such a "Meissner effect" in general. This would seem to imply that a vortex does not want to be attached to an extreme black hole. We calculate the total energy of the vortex fields in the presence of an extreme black hole. When the hole is small relative to the size of the vortex, it is energetically favored for the hole to remain inside the vortex region, contrary to the intuition that the hole should be expelled. However, as we allow the extreme horizon radius to become very large compared to the radius of the vortex, we do find evidence of an instability. This proves that it is energetically unfavorable for a thin vortex to interact with a large extreme black hole. This would seem to dispel the notion that a black hole can support "long" Abelian Higgs hair in the extreme limit. We show that these considerations do not go through in the near-extreme limit. Finally, we discuss the implications for strings that end at black holes, as in the processes where a string snaps by nucleating black holes.

The seasonal evolution of high vertical–mode internal waves in a deep reservoir

The causal mechanism and seasonal evolution of the internal wave field in a deep, warm, monomictic reservoirare examined through the analysis of field observations and numerical techniques. The study period extends fromthe onset of thermal stratification in the spring until midsummer in 2005. During this time, wind forcing wasperiodic, with a period of 24 h (typical of land–sea breezes), and the thermal structure in the lake wascharacterized by the presence of a shallow surface layer overlying a thick metalimnion, typical of small to mediumsized reservoirs with deep outtakes. Basin-scale internal seiches of high vertical mode (ranging from mode V3 toV5) were observed in the metalimnion. The structure of the dominant modes of oscillation changed asstratification evolved on seasonal timescales, but in all cases, their periods were close to that of the local windforcing (i.e., 24 h), suggesting a resonant response. Nonresonant oscillatory modes of type V1 and V2 becamedominant after large frontal events, which disrupted the diurnal periodicity of the wind forcing

Quatro alternativas para resolver a equação de Schrödinger para o átomo de hidrogênio

Quantum chemistry describes the hydrogen atom as one of the few systems that permits an exact solution of the Schrödinger equation. Students tend to consider that little can be learned from the hydrogen atom and forget that it can be used as a standard to test numerical procedures used to calculate properties of multielectronic systems. In this paper, four different numerical procedures are described in order to solve the Schrödinger equation for the hydrogen atom. The basic motivation is to identify new insights and methods that can be obtained from the application of powerful numerical techniques in a well-known system.

Multidata inverse problems and Bayesian solution methods in astronomy

Statistical analyses of measurements that can be described by statistical models are of essence in astronomy and in scientific inquiry in general. The sensitivity of such analyses, modelling approaches, and the consequent predictions, is sometimes highly dependent on the exact techniques applied, and improvements therein can result in significantly better understanding of the observed system of interest. Particularly, optimising the sensitivity of statistical techniques in detecting the faint signatures of low-mass planets orbiting the nearby stars is, together with improvements in instrumentation, essential in estimating the properties of the population of such planets, and in the race to detect Earth-analogs, i.e. planets that could support liquid water and, perhaps, life on their surfaces. We review the developments in Bayesian statistical techniques applicable to detections planets orbiting nearby stars and astronomical data analysis problems in general. We also discuss these techniques and demonstrate their usefulness by using various examples and detailed descriptions of the respective mathematics involved. We demonstrate the practical aspects of Bayesian statistical techniques by describing several algorithms and numerical techniques, as well as theoretical constructions, in the estimation of model parameters and in hypothesis testing. We also apply these algorithms to Doppler measurements of nearby stars to show how they can be used in practice to obtain as much information from the noisy data as possible. Bayesian statistical techniques are powerful tools in analysing and interpreting noisy data and should be preferred in practice whenever computational limitations are not too restrictive.

Projection-Based Statistical Inference in Linear Structural Models with Possibly Weak Instruments

It is well known that standard asymptotic theory is not valid or is extremely unreliable in models with identification problems or weak instruments [Dufour (1997, Econometrica), Staiger and Stock (1997, Econometrica), Wang and Zivot (1998, Econometrica), Stock and Wright (2000, Econometrica), Dufour and Jasiak (2001, International Economic Review)]. One possible way out consists here in using a variant of the Anderson-Rubin (1949, Ann. Math. Stat.) procedure. The latter, however, allows one to build exact tests and confidence sets only for the full vector of the coefficients of the endogenous explanatory variables in a structural equation, which in general does not allow for individual coefficients. This problem may in principle be overcome by using projection techniques [Dufour (1997, Econometrica), Dufour and Jasiak (2001, International Economic Review)]. AR-types are emphasized because they are robust to both weak instruments and instrument exclusion. However, these techniques can be implemented only by using costly numerical techniques. In this paper, we provide a complete analytic solution to the problem of building projection-based confidence sets from Anderson-Rubin-type confidence sets. The latter involves the geometric properties of “quadrics” and can be viewed as an extension of usual confidence intervals and ellipsoids. Only least squares techniques are required for building the confidence intervals. We also study by simulation how “conservative” projection-based confidence sets are. Finally, we illustrate the methods proposed by applying them to three different examples: the relationship between trade and growth in a cross-section of countries, returns to education, and a study of production functions in the U.S. economy.

Sur un système de deux oscillateurs FitzHugh-Nagumo couplés

Ce mémoire consiste en l’étude du comportement dynamique de deux oscillateurs FitzHugh-Nagumo identiques couplés. Les paramètres considérés sont l’intensité du courant injecté et la force du couplage. Juqu’à cinq solutions stationnaires, dont on analyse la stabilité asymptotique, peuvent co-exister selon les valeurs de ces paramètres. Une analyse de bifurcation, effectuée grâce à des méthodes tant analytiques que numériques, a permis de détecter différents types de bifurcations (point de selle, Hopf, doublement de période, hétéroclinique) émergeant surtout de la variation du paramètre de couplage. Une attention particulière est portée aux conséquences de la symétrie présente dans le système.

Ultrasonic Study Of The Elastic Properties And Phase Transitions In Selected Mixed Sulphate Crystals

The thesis investigated the elastic properties and phase transitions in selected mixed sulphate crystals – Lithium Hydrazinium Sulphate [LiN2H2SO4], Lithium Ammonium Sulphate [LiNH4SO4] and Lithium Potassium Sulphate [LiKSO4] – using ultrasonic technique. The pulse echo overlap technique has been used for measuring ultrasonic velocity and its dependence on temperature along different directions with waves of longitudinal and transverse polarizations. Two major numerical techniques and the corresponding computer programs developed as part of present work are presented in this thesis. All the 9 elastic constants of LHS are determined accurately from ultrasonic measurements and applying misorientation correction refines the constants. Ultrasonic measurements are performed in LAS to determine the elastic constants and to study the low temperature phase transitions. Temperature variation studies of elastic constant of LAS are performed for 6 different modes of propagation for heating and cooling at low temperatures. All the 5 independent elastic constants of LPS is determined using ultrasonic measurements. It is concluded that LPS crystal does not undergo a phase transition near this temperature. A comparison of the three crystals studied shows that LPS has maximum number of phase transitions and LHS has the least number. It is interesting to note that LPS has the simplest formula unit among the three. There is considerable scope for the future work on these crystals and others belonging to the sulphate family.