Temperature-dependent vortex motion in a square mesoscopic superconducting cylinder: Ginzburg-Landau calculations

In this work we investigate the dynamics of vortices in a square mesoscopic superconductor. As time evolves we show how the vortices are nucleated into the sample to form a multivortex, single vortex, and giant vortex states. We illustrate how the vortices move around at the transition fields before they accommodate into an equilibrium configuration. We also calculate the magnetization and the free energy as functions of the applied magnetic field for several values of temperature. In addition, we evaluate the upper critical field.

Monte Carlo and modified Tanford-Kirkwood results for macromolecular electrostatics calculations

The understanding of electrostatic interactions is an essential aspect of the complex correlation between structure and function of biological macromolecules. It is also important in protein engineering and design. Theoretical studies of such interactions are predominantly done within the framework of Debye-Huckel theory. A classical example is the Tanford-Kirkwood (TK) model. Besides other limitations, this model assumes an infinitesimally small macromolecule concentration. By comparison to Monte Carlo (MC) simulations, it is shown that TK predictions for the shifts in ion binding constants upon addition of salt become less reliable even at moderately macromolecular concentrations. A simple modification based on colloidal literature is suggested to the TK scheme. The modified TK models suggested here satisfactorily predict MC and experimental shifts in the calcium binding constant as a function of protein concentration for the calbindin D-9k mutant and calmodulin.

The STATFLUX code: a statistical method for calculation of flow and set of parameters, based on the Multiple-Compartment Biokinetical Model

The code STATFLUX, implementing a new and simple statistical procedure for the calculation of transfer coefficients in radionuclide transport to animals and plants, is proposed. The method is based on the general multiple-compartment model, which uses a system of linear equations involving geometrical volume considerations. Flow parameters were estimated by employing two different least-squares procedures: Derivative and Gauss-Marquardt methods, with the available experimental data of radionuclide concentrations as the input functions of time. The solution of the inverse problem, which relates a given set of flow parameter with the time evolution of concentration functions, is achieved via a Monte Carlo Simulation procedure.Program summaryTitle of program: STATFLUXCatalogue identifier: ADYS_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADYS_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneComputer for which the program is designed and others on which it has been tested: Micro-computer with Intel Pentium III, 3.0 GHzInstallation: Laboratory of Linear Accelerator, Department of Experimental Physics, University of São Paulo, BrazilOperating system: Windows 2000 and Windows XPProgramming language used: Fortran-77 as implemented in Microsoft Fortran 4.0. NOTE: Microsoft Fortran includes non-standard features which are used in this program. Standard Fortran compilers such as, g77, f77, ifort and NAG95, are not able to compile the code and therefore it has not been possible for the CPC Program Library to test the program.Memory, required to execute with typical data: 8 Mbytes of RAM memory and 100 MB of Hard disk memoryNo. of bits in a word: 16No. of lines in distributed program, including test data, etc.: 6912No. of bytes in distributed Program, including test data, etc.: 229 541Distribution format: tar.gzNature of the physical problem: the investigation of transport mechanisms for radioactive substances, through environmental pathways, is very important for radiological protection of populations. One such pathway, associated with the food chain, is the grass-animal-man sequence. The distribution of trace elements in humans and laboratory animals has been intensively studied over the past 60 years [R.C. Pendlenton, C.W. Mays, R.D. Lloyd, A.L. Brooks, Differential accumulation of iodine-131 from local fallout in people and milk, Health Phys. 9 (1963) 1253-1262]. In addition, investigations on the incidence of cancer in humans, and a possible causal relationship to radioactive fallout, have been undertaken [E.S. Weiss, M.L. Rallison, W.T. London, W.T. Carlyle Thompson, Thyroid nodularity in southwestern Utah school children exposed to fallout radiation, Amer. J. Public Health 61 (1971) 241-249; M.L. Rallison, B.M. Dobyns, F.R. Keating, J.E. Rall, F.H. Tyler, Thyroid diseases in children, Amer. J. Med. 56 (1974) 457-463; J.L. Lyon, M.R. Klauber, J.W. Gardner, K.S. Udall, Childhood leukemia associated with fallout from nuclear testing, N. Engl. J. Med. 300 (1979) 397-402]. From the pathways of entry of radionuclides in the human (or animal) body, ingestion is the most important because it is closely related to life-long alimentary (or dietary) habits. Those radionuclides which are able to enter the living cells by either metabolic or other processes give rise to localized doses which can be very high. The evaluation of these internally localized doses is of paramount importance for the assessment of radiobiological risks and radiological protection. The time behavior of trace concentration in organs is the principal input for prediction of internal doses after acute or chronic exposure. The General Multiple-Compartment Model (GMCM) is the powerful and more accepted method for biokinetical studies, which allows the calculation of concentration of trace elements in organs as a function of time, when the flow parameters of the model are known. However, few biokinetics data exist in the literature, and the determination of flow and transfer parameters by statistical fitting for each system is an open problem.Restriction on the complexity of the problem: This version of the code works with the constant volume approximation, which is valid for many situations where the biological half-live of a trace is lower than the volume rise time. Another restriction is related to the central flux model. The model considered in the code assumes that exist one central compartment (e.g., blood), that connect the flow with all compartments, and the flow between other compartments is not included.Typical running time: Depends on the choice for calculations. Using the Derivative Method the time is very short (a few minutes) for any number of compartments considered. When the Gauss-Marquardt iterative method is used the calculation time can be approximately 5-6 hours when similar to 15 compartments are considered. (C) 2006 Elsevier B.V. All rights reserved.

THEORY OF ELECTRICAL CHARACTERISTICS OF A SCHOTTKY-BARRIER HAVING EXPONENTIALLY DISTRIBUTED IMPURITY STATES AND METAL-INSULATOR METAL STRUCTURES

The metal-insulator or metal-amorphous semiconductor blocking contact is still not well understood. Here, we discuss the steady state characteristics of a non-intimate metal-insulator Schottky barrier. We consider an exponential distribution (in energy) of impurity states in addition to impurity states at a single energy level within the depletion region. We present analytical expressions for the electrical potential, field, thickness of depletion region, capacitance, and charge accumulated in the depletion region. We also discuss ln I versus V(ap) data. Finally, we compare the characteristics in three cases: (i) impurity states at only a single energy level; (ii) uniform energy distribution of impurity states; and (iii) exponential energy distribution of impurity states.In general, the electrical characteristics of Schottky barriers and metal-insulator-metal structures with Schottky barriers depend strongly on the energy distribution of impurity states.

THEORY OF NON-STEADY-STATE ELECTRICAL CHARACTERISTICS OF METAL-INSULATOR-METAL STRUCTURES WITH SCHOTTKY BARRIERS AND UNIFORMLY DISTRIBUTED INTERFACE IMPURITY STATES

The metal-insulator (or amorphous semiconductor) blocking contact is still not well understood. In the present paper, we discuss the non steady state characteristics of Metal-lnsulator-Metal Structure with non-intimate blocking contacts (i.e. Metal-Oxide-Insulator-Metal Structure). We consider a uniform distribution (in energy) of impurity states in addition to impurity states at a single energy level within the depletion region. We discuss thermal as well as isothermal characteristics and present expressions for the temperature of maximum current (T-m) and a method to calculate the density of uniformly distributed impurity states. The variation of mobility with electrical field has also been considered. Finally we plot the theoretical curves under different conditions. The present results are closing into available experimental results.

Simulation of free surface flows in a distributed memory environment

A parallel technique, for a distributed memory machine, based on domain decomposition for solving the Navier-Stokes equations in cartesian and cylindrical coordinates in two dimensions with free surfaces is described. It is based on the code by Tome and McKee (J. Comp. Phys. 110 (1994) 171-186) and Tome (Ph.D. Thesis, University of Strathclyde, Glasgow, 1993) which in turn is based on the SMAC method by Amsden and Harlow (Report LA-4370, Los Alamos Scientific Laboratory, 1971), which solves the Navier-Stokes equations in three steps: the momentum and Poisson equations and particle movement, These equations are discretized by explicit and 5-point finite differences. The parallelization is performed by splitting the computation domain into vertical panels and assigning each of these panels to a processor. All the computation can then be performed using nearest neighbour communication. Test runs comparing the performance of the parallel with the serial code, and a discussion of the load balancing question are presented. PVM is used for communication between processes. (C) 1999 Elsevier B.V. B.V. All rights reserved.

A combined wavelet-element free Galerkin method for numerical calculations of electromagnetic fields

A combined wavelet-element free Galerkin (EFG) method is proposed for solving electromagnetic EM) field problems. The bridging scales are used to preserve the consistency and linear independence properties of the entire bases. A detailed description of the development of the discrete model and its numerical implementations is given to facilitate the reader to. understand the proposed algorithm. A numerical example to validate the proposed method is also reported.

Crystal and molecular structures of 4-substituted 3,4-dihydropyrimidin-2(1H)-ones studied by X-ray and AM1 and B3LYP calculations

(1) C13H13N3O5, Mr = 291.26, P (1) over bar, a = 7.4629(9), b = 7.9203(9), c = 12.126(2) angstrom, alpha = 86.804(5), beta = 78.471(7), gamma = 69.401(8)degrees, V = 657.3(2)angstrom(3), Z = 2, R-1 = 0.0454; (2) C11H12N2O4, Mr=236.23, Pbca, a=7.2713(9), b=14.234(1), c=20.848(3)angstrom, V= 2157.8(4) angstrom(3), Z=8, R-1=0.0504; (3) C13H13N2O3Cl, Mr = 280.70, P2/n, a = 17.344(2), b = 9.237(1), c = 18.398(2) angstrom; beta = 92.61(2)degrees, V = 2944.4(6) angstrom(3), Z = 8, R-1 = 0.0714. The conformational features of three 4-substituted-3-4-dihydropyrimidin-2(1H)-ones were investigated by computational and single crystal X-ray crystallographic studies. The geometries were optimized using semiempirical (AM1) and first principle calculations (B3LYP/6-31G**) methods, the rotational barriers for important functional groups were studied. In all structures the pyrimidinone rings are in a more or less distorted boat conformation. The phenyl and the furane rings are almost perpendicular to the best least-squares plane through the dihydropyrimidinone ring.

Distributed feedback multipeak laser emission in Rhodamine 6G doped organic-inorganic hybrids

Organic-inorganic hybrid materials were prepared from an ureasil precursor (ureapropyltriethoxysilane designated as UPTES) and acrylic acid modified zirconium (IV) n-propoxide. Thin films containing rhodamine 6G (Rh6G) were prepared by spin-coating on glass substrates with different Zr:Si molar ratios (Zr:Si = 75:25, 50:50 and 25:75). Refractive index, thickness, number of propagating modes and attenuation coefficient were measured at 543.5, 632.8 and 1550 nm wavelengths by the prism coupling technique. Distributed feedback (DFB) laser effect was observed and studied as a function of films thickness and refractive index.

Nonadiabatic calculations of the oscillator strengths for the helium atom in the hyperspherical adiabatic approach

Energies and wavefunctions are calculated for the bound states of the helium atom in the hyperspherical adiabatic approach by the full inclusion of nonadiabatic couplings. We show that the use of appropriate asymptotic radial boundary conditions not only allows the efficient calculation of energies accurate up to a few ppm for the ground state but also gives increasingly precise results for high-lying excited states with a unique set of equations. The accuracy of the wavefunctions is demonstrated by the calculation of oscillator strengths in the length form for transitions between stares ii S-1(e) and (n + 1) P-1(0) up to n = 29, in agreement with variational calculations.

Structure characterization of molecular complexes for non-linear optical materials. II. 1 : 1 complexes of 4-methyl-morpholine-N-oxide (1) and 3-picoline-N-oxide (2) with 2,4,6-trinitrophenol, studied by X-ray, AM1 and DFT calculations

(1) C6H2N3O7- center dot C5H12NO2+, Mr = 346.26, P2(1)/c, a = 7.2356(6), b = 10.5765(9), c = 19.593(2) angstrom, 3 beta=95.101(6)degrees, V = 1493.5(2) angstrom(3), Z = 4, R-1 = 0.0414; (2) C6H2N3O7- center dot C6H8NO+, Mr = 38.24, P2(1)/n, a = 7.8713(5), b = 6.1979(7), c = 28.697(3) angstrom, beta = 90.028(7)degrees, V = 1400.0(2) angstrom(3), Z = 4, R-1 = 0.0416. The packing units in both compounds consist of hydrogen bonded cation-anion pairs. The (hyper)polarizabilities have been calculated for the crystallographic and optimized molecules, by AM1 and at the DFT/B3LYP(6-31G**) level.

Development of basis sets to calculations of the electronic structure of YMnO3

The generator coordinate Hartree-Fock method was used to develop 20s17p, 30s20p14d, and 30s21p16d Gaussian basis sets for the O ((3)p), Mn (S-6), and Y (D-2) atoms, respectively. The Gaussian basis sets were contracted to 20s17p/9s7p, 30s20p14d/11s7p7d, and 30s21p16d/14s7p7d and utilized in calculations of total energy and orbital energies of the (MnO1+)-Mn-5 and (YO1+)-Y-3 fragments to evaluate its quality in molecular studies. Finally, the contracted basis set for O atom was supplemented with one polarization function of d symmetry and used along with the other contracted basis sets (for Mn and Y) to calculate dipole moments, total energy, and total atomic charges in YMnO3 in space group D-6h. The analysis of those properties showed that is reasonable to believe that YMnO3 present behavior of piezoelectric material. (C) 2003 Elsevier B.V. All rights reserved.

A revised taxonomy for time warp based distributed synchronization protocols

This paper presents an improved and updated taxonomy for Time Warp based distributed synchronization protocols. This taxonomy aims to allow the grouping of several optmistic distributed simulation synchronization protocols, with the objective to facilitate the task to decide which protocol is better for a specific simulation.

Evaluating distributed time-varying generation through a multiobjective index

In the last decade, distributed generation, with its various technologies, has increased its presence in the energy mix presenting distribution networks with challenges in terms of evaluating the technical impacts that require a wide range of network operational effects to be qualified and quantified. The inherent time-varying behavior of demand and distributed generation (particularly when renewable sources are used), need to be taken into account since considering critical scenarios of loading and generation may mask the impacts. One means of dealing with such complexity is through the use of indices that indicate the benefit or otherwise of connections at a given location and for a given horizon. This paper presents a multiobjective performance index for distribution networks with time-varying distributed generation which consider a number of technical issues. The approach has been applied to a medium voltage distribution network considering hourly demand and wind speeds. Results show that this proposal has a better response to the natural behavior of loads and generation than solely considering a single operation scenario.

Gaussian basis sets by generator coordinate Hartree-Fock method to ab initio calculations of electron affinities of enolates

The Generator Coordinate Hartree-Fock (GCHF) method is employed to generate uncontracted 15s and 18s11p gaussian basis sets for the H, C and O atoms, respectively. These basis sets are then contracted to 3s and 4s H atom and 6s5p, for C and O atoms by a standard procedure. For quality evaluation of contracted basis sets in molecular calculations, we have accomplished calculations of total and orbital energies in the Hartree-Fock-Roothaaii (HFR) approach for CH, C(2) and CO molecules. The results obtained with the uncontracted basis sets are compared with values obtained with the standard D95, 6-311G basis sets and with values reported in the literature. The 4s and 6s5p basis sets are enriched with polarization and diffuse functions for atoms of the parent neutral systems and of the enolates anions (cycloheptanone enolate, 2,5-dimethyleyelopentanone enolate, 4-heptanone enolate, and di-isopropyl ketone enolate) from the literature, in order to assess their performance in ab initio molecular calculations, and applied for calculations of electron affinities of the enolates. The calculations were performed at the DFT (BLYP and B3LYP) and HF levels and compared with the corresponding experimental values and with those obtained by using other 6-3 1 + +G((*)) and 6-311 + +G((*)) basis sets from literature. For the enolates studied, the differences between the electron affinities obtained with GCHF basis sets, at the B3LYP level, and the experimental values are -0.001, -0,014, -0.001, and -0.001 eV. (C) 2002 Elsevier B.V. B.V. All rights reserved.