Numerical calculation of the left and right fractional derivatives

The calculation of fractional derivatives is an important topic in scientific research. While formal definitions are clear from the mathematical point of view, they pose limitations in applied sciences that have not been yet tackled. This paper addresses the problem of obtaining left and right side derivatives when adopting numerical approximations. The results reveal the relationship between the resulting distinct values for different fractional orders and types of signals.

Development of a method for numerical calculation of wave refraction /

"October 1964."

Numerical calculation of diffraction coefficients from building edges using FDTD method

Finite Difference Time Domain (FDTD) Method and software are applied to obtain diffraction waves from modulated Gaussian plane wave illumination for right angle wedges and Fast Fourier Transform (FFT) is used to get diffraction coefficients in a wideband in the illuminated lit region. Theta and Phi polarization in 3-dimensional, TM and TE polarization in 2-dimensional cases are considered respectively for soft and hard diffraction coefficients. Results using FDTD method of perfect electric conductor (PEC) wedge are compared with asymptotic expressions from Uniform Theory of Diffraction (UTD). Extend the PEC wedges to some homogenous conducting and dielectric building materials for diffraction coefficients that are not available analytically in practical conditions. ^

Quantum Monte Carlo : some theoretical and numerical studies

In Part I, theoretical derivations for Variational Monte Carlo calculations are compared with results from a numerical calculation of He; both indicate that minimization of the ratio estimate of Evar , denoted EMC ' provides different optimal variational parameters than does minimization of the variance of E MC • Similar derivations for Diffusion Monte Carlo calculations provide a theoretical justification for empirical observations made by other workers. In Part II, Importance sampling in prolate spheroidal coordinates allows Monte Carlo calculations to be made of E for the vdW molecule var He2' using a simplifying partitioning of the Hamiltonian and both an HF-SCF and an explicitly correlated wavefunction. Improvements are suggested which would permit the extension of the computational precision to the point where an estimate of the interaction energy could be made~

Numerical loop calculations with contour deformation

Gegenstand dieser Arbeit ist die nummerische Berechnung von Schleifenintegralen welche in höheren Ordnungen der Störungstheorie auftreten.rnAnalog zur reellen Emission kann man auch in den virtuellen Beiträgen Subtraktionsterme einführen, welche die kollinearen und soften Divergenzen des Schleifenintegrals entfernen. Die Phasenraumintegration und die Schleifenintegration können dann in einer einzigen Monte Carlo Integration durchgeführt werden. In dieser Arbeit zeigen wir wie eine solche numerische Integration unter zu Hilfenahme einer Kontourdeformation durchgeführt werden kann. Ausserdem zeigen wir wie man die benötigeten Integranden mit Rekursionsformeln berechnen kann.

Analytical and numerical study of polarons

The aim of this thesis is to introduce the polaron concept and to perform a DFT numerical calculation of a small polaron in the rutile phase of TiO2. In the first chapters, we present an analytical study of small and large polarons, based on the Holstein and Fröhlich Hamiltonians. The necessary mathematical formalism and physics fundamentals are briefly reviewed in the first chapter. In the second part of the thesis, Density Functional Theory (DFT) is introduced together with the DFT+U correction and its implementation in the Vienna Ab-Initio Simulation Package (VASP). The calculation of a small polaron in rutile is then described and discussed at a qualitative level. The polaronic solution is compared with the one of a delocalized electron. The calculation showed how the polaron creates a new energy level 0.70 eV below the conduction band. The energy level is visible both in the band structure diagram and in the density of states diagram. The electron is localized on a titanium atom, distorting the surrounding lattice. In particular, the four oxygen atoms closer to the titanium atom are displaced by 0.085 Å outwards, whereas the two further oxygen atoms by 0.023 Å. The results are compatible, at a qualitative level, with the literature. Further developments of this work may try to improve the precision of the results and to quantitatively compare them with the literature.

On the anharmonic, multi-phonon, Debye-Waller contributions to the phonon-limited resistivity of metals : applications to Na and K

The anharmonic, multi-phonon (MP), and Oebye-Waller factor (OW) contributions to the phonon limited resistivity (;0) of metals derived by Shukla and Muller (1979) by the doubletime temperature dependent Green function method have been numerically evaluated for Na and K in the high temperature limit. The anharmonic contributions arise from the cubic and quartic shift of phonons (CS, QS), and phonon width (W) and the interference term (1). The QS, MP and OW contributions to I' are also derived by the matrix element method and the results are in agreement with those of Shukla and Muller (1979). In the high temperature limit, the contributions to;O from each of the above mentioned terms are of the type BT2 For numerical calculations suitable expressions are derived for the anharmonic contributions to ~ in terms of the third and fourth rank tensors obtained by the Ewald procedure. The numerical calculation of the contributions to;O from the OW, MP term and the QS have been done exactly and from the CS, Wand I terms only approximately in the partial and total Einstein approximations (PEA, TEA), using a first principle approach (Shukla and Taylor (1976)). The results obtained indicate that there is a strong pairwise cancellation between the: OW and MP terms, the QS and CS and the Wand I terms. The sum total of these contributions to;O for Na and K amounts to 4 to 11% and 2 to 7%, respectively, in the PEA while in the TEA they amount to 3 to 7% and 1 to 4%, respectively, in the temperature range.

Énergie des vortex dans un modèle abélien de Higgs en 2+1 dimensions avec un potentiel d’ordre six

Dans ce travail, j’étudierai principalement un modèle abélien de Higgs en 2+1 dimensions, dans lequel un champ scalaire interagit avec un champ de jauge. Des défauts topologiques, nommés vortex, sont créés lorsque le potentiel possède un minimum brisant spontanément la symétrie U(1). En 3+1 dimensions, ces vortex deviennent des défauts à une dimension. Ils ap- paraissent par exemple en matière condensée dans les supraconducteurs de type II comme des lignes de flux magnétique. J’analyserai comment l’énergie des solutions statiques dépend des paramètres du modèle et en particulier du nombre d’enroulement du vortex. Pour le choix habituel de potentiel (un poly- nôme quartique dit « BPS »), la relation entre les masses des deux champs mène à deux types de comportements : type I si la masse du champ de jauge est plus grande que celle du champ sca- laire et type II inversement. Selon le cas, la dépendance de l’énergie au nombre d’enroulement, n, indiquera si les vortex auront tendance à s’attirer ou à se repousser, respectivement. Lorsque le flux emprisonné est grand, les vortex présentent un profil où la paroi est mince, permettant certaines simplifications dans l’analyse. Le potentiel, un polynôme d’ordre six (« non-BPS »), est choisi tel que le centre du vortex se trouve dans le vrai vide (minimum absolu du potentiel) alors qu’à l’infini le champ scalaire se retrouve dans le faux vide (minimum relatif du potentiel). Le taux de désintégration a déjà été estimé par une approximation semi-classique pour montrer l’impact des défauts topologiques sur la stabilité du faux vide. Le projet consiste d’abord à établir l’existence de vortex classi- quement stables de façon numérique. Puis, ma contribution fut une analyse des paramètres du modèle révélant le comportement énergétique de ceux-ci en fonction du nombre d’enroulement. Ce comportement s’avèrera être différent du cas « BPS » : le ratio des masses ne réussit pas à décrire le comportement observé numériquement.

Water wave scattering by finite arrays of circular structures

The scattering of small amplitude water waves by a finite array of locally axisymmetric structures is considered. Regions of varying quiescent depth are included and their axisymmetric nature, together with a mild-slope approximation, permits an adaptation of well-known interaction theory which ultimately reduces the problem to a simple numerical calculation. Numerical results are given and effects due to regions of varying depth on wave loading and free-surface elevation are presented.

Comment on “A modified leapfrog scheme for shallow water equations” by Wen-Yih Sun and Oliver M.T. Sun

A recent paper published in this journal considers the numerical integration of the shallow-water equations using the leapfrog time-stepping scheme [Sun Wen-Yih, Sun Oliver MT. A modified leapfrog scheme for shallow water equations. Comput Fluids 2011;52:69–72]. The authors of that paper propose using the time-averaged height in the numerical calculation of the pressure-gradient force, instead of the instantaneous height at the middle time step. The authors show that this modification doubles the maximum Courant number (and hence the maximum time step) at which the integrations are stable, doubling the computational efficiency. Unfortunately, the pressure-averaging technique proposed by the authors is not original. It was devised and published by Shuman [5] and has been widely used in the atmosphere and ocean modelling community for over 40 years.

Thermoluminescence, electron paramagnetic resonance and optical absorption in natural and synthetic rhodonite crystals

Thermoluminescence, electron paramagnetic resonance and optical absorption properties of rhodonite, a natural silicate mineral, have been investigated and compared to those of synthetic crystal, pure and doped. The TL peaks grow linearly for radiation dose up to 4 kGy, and then saturate. In all the synthetic samples, 140 and 340 degrees C TL peaks are observed; the difference occurs in their relative intensities, but only 340 degrees C peak grows strongly for high doses. Al(2)O(3) and Al(2)O(3) + CaO-doped synthetic samples presented several decades intenser TL compared to that of synthetic samples doped with other impurities. A heating rate of 4 degrees C/s has been used in all the TL readings. The EPR spectrum of natural rhodonite mineral has only one huge signal around g = 2.0 with width extending from 1,000 to 6,000 G. This is due to Mn dipolar interaction, a fact proved by numerical calculation based on Van Vleck dipolar broadening expression. The optical absorption spectrum is rich in absorption bands in near-UV, visible and near-IR intervals. Several bands in the region from 540 to 340 nm are interpreted as being due to Mn(3+) in distorted octahedral environment. A broad and intense band around 1,040 nm is due to Fe(2+). It decays under heating up to 900 degrees C. At this temperature it is reduced by 80% of its original intensity. The pink, natural rhodonite, heated in air starts becoming black at approximately 600 degrees C.

Understanding the tsunami with a simple model

In this paper, we use the approximation of shallow water waves (Margaritondo G 2005 Eur. J. Phys. 26 401) to understand the behaviour of a tsunami in a variable depth. We deduce the shallow water wave equation and the continuity equation that must be satisfied when a wave encounters a discontinuity in the sea depth. A short explanation about how the tsunami hit the west coast of India is given based on the refraction phenomenon. Our procedure also includes a simple numerical calculation suitable for undergraduate students in physics and engineering.

Simulação numérica do fluxo hídrico subterrâneo na bacia hidrográfica do Rio Pitimbu RN

Urban centers in Pitimbu Watershed use significant groundwater sources for public supply. Therefore, studies in Dunas Barreiras aquifer are relevant to expand knowledge about it and help manage water resources in the region. An essential tool for this management is the numerical modeling of groundwater flow. In this work, we developed a groundwater flow model for Pitimbu Watershed, using the Visual Modflow, version 2.7.1., which uses finite difference method for solving the govern equation of the dynamics of groundwater flow. We carried out the numerical simulation of steady-state model for the entire region of the basin. The model was built in the geographical, geomorphological and hydrogeological study of the area, which defined the boundary conditions and the parameters required for the numerical calculation. Owing to unavailability of current data based on monitoring of the aquifer it was not possible to calibrate the model. However, the simulation results showed that the overall water balance approached zero, therefore satisfying the equation for the three-dimensional behavior of the head water in steady state. Variations in aquifer recharge data were made to verify the impact of this contribution on the water balance of the system, especially in the scenario in which recharge due to drains and sinks was removed. According to the results generated by Visual Modflow occurred significantly hydraulic head lowering, ranging from 16,4 to 82 feet of drawdown. With the results obtained, it can be said that modeling is performed as a valid tool for the management of water resources in Pitimbu River Basin, and to support new studies

Propagação de danos em sistemas cooperativos: propriedades termodinâmicas

High-precision calculations of the correlation functions and order parameters were performed in order to investigate the critical properties of several two-dimensional ferro- magnetic systems: (i) the q-state Potts model; (ii) the Ashkin-Teller isotropic model; (iii) the spin-1 Ising model. We deduced exact relations connecting specific damages (the difference between two microscopic configurations of a model) and the above mentioned thermodynamic quanti- ties which permit its numerical calculation, by computer simulation and using any ergodic dynamics. The results obtained (critical temperature and exponents) reproduced all the known values, with an agreement up to several significant figures; of particular relevance were the estimates along the Baxter critical line (Ashkin-Teller model) where the exponents have a continuous variation. We also showed that this approach is less sensitive to the finite-size effects than the standard Monte-Carlo method. This analysis shows that the present approach produces equal or more accurate results, as compared to the usual Monte Carlo simulation, and can be useful to investigate these models in circumstances for which their behavior is not yet fully understood

Espectros fractais em sistemas nanoestruturados e cristais fotônicos

The study of the elementary excitations such as photons, phonons, plasmons, polaritons, polarons, excitons and magnons, in crystalline solids and nanostructures systems are nowdays important active field for research works in solid state physics as well as in statistical physics. With this aim in mind, this work has two distinct parts. In the first one, we investigate the propagation of excitons polaritons in nanostructured periodic and quasiperiodic multilayers, from the description of the behavior for bulk and surface modes in their individual constituents. Through analytical, as well as computational numerical calculation, we obtain the spectra for both surface and bulk exciton-polaritons modes in the superstructures. Besides, we investigate also how the quasiperiodicity modifies the band structure related to the periodic case, stressing their amazing self-similar behavior leaving to their fractal/multifractal aspects. Afterwards, we present our results related to the so-called photonic crystals, the eletromagnetic analogue of the electronic crystalline structure. We consider periodic and quasiperiodic structures, in which one of their component presents a negative refractive index. This unusual optic characteristic is obtained when the electric permissivity and the magnetic permeability µ are both negatives for the same range of angular frequency ω of the incident wave. The given curves show how the transmission of the photon waves is modified, with a striking self-similar profile. Moreover, we analyze the modification of the usual Planck´s thermal spectrum when we use a quasiperiodic fotonic superlattice as a filter.