The gravitational field of a global monopole

We present an exact solution to the non-linear equation which describes a global monopole in the flat space. We re-examine the metric and the geodesics outside the global monopole. We will see that a global monopole produces a repulsive gravitational field outside the core in addition to a solid angular deficit. The lensing property of the global monopole and the global monopole-antimonopole annihilation mechanism are studied.Comment: 8 pages, no figure

The Carriers of the Interstellar Unidentified Infrared Emission Features: Constraints from the Interstellar C-H Stretching Features at 3.2-3.5 Micrometers

The unidentified infrared emission (UIE) features at 3.3, 6.2, 7.7, 8.6, and 11.3 micrometer, commonly attributed to polycyclic aromatic hydrocarbon (PAH) molecules, have been recently ascribed to mixed aromatic/aliphatic organic nanoparticles. More recently, an upper limit of <9% on the aliphatic fraction (i.e., the fraction of carbon atoms in aliphatic form) of the UIE carriers based on the observed intensities of the 3.4 and 3.3 micrometer emission features by attributing them to aliphatic and aromatic C-H stretching modes, respectively, and assuming A_34./A_3.3~0.68 derived from a small set of aliphatic and aromatic compounds, where A_3.4 and A_3.3 are respectively the band strengths of the 3.4 micrometer aliphatic and 3.3 micrometer aromatic C-H bonds. To improve the estimate of the aliphatic fraction of the UIE carriers, here we analyze 35 UIE sources which exhibit both the 3.3 and 3.4 micrometer C-H features and determine I_3.4/I_3.3, the ratio of the power emitted from the 3.4 micrometer feature to that from the 3.3 micrometer feature. We derive the median ratio to be ~ 0.12. We employ density functional theory and second-order perturbation theory to compute A_3.4/A_3.3 for a range of methyl-substituted PAHs. The resulting A_3.4/A_3.3 ratio well exceeds 1.4, with an average ratio of ~1.76. By attributing the 3.4 micrometer feature exclusively to aliphatic C-H stretch (i.e., neglecting anharmonicity and superhydrogenation), we derive the fraction of C atoms in aliphatic form to be ~2%. We therefore conclude that the UIE emitters are predominantly aromatic.Comment: 14 pages, 5 figures, 1 table; accepted for publication in The Astrophysical Journa

The Hamiltonian index of a graph and its branch-bonds

Let $G$ be an undirected and loopless finite graph that is not a path. The minimum $m$ such that the iterated line graph $L^m(G)$ is hamiltonian is called the hamiltonian index of $G,$ denoted by $h(G).$ A reduction method to determine the hamiltonian index of a graph $G$ with $h(G)\geq 2$ is given here. With it we will establish a sharp lower bound and a sharp upper bound for $h(G)$, respectively, which improves some known results of P.A. Catlin et al. [J. Graph Theory 14 (1990)] and H.-J. Lai [Discrete Mathematics 69 (1988)]. Examples show that $h(G)$ may reach all integers between the lower bound and the upper bound. \u

Flexible protein folding by ant colony optimization

Protein structure prediction is one of the most challenging topics in bioinformatics. As the protein structure is found to be closely related to its functions, predicting the folding structure of a protein to judge its functions is meaningful to the humanity. This chapter proposes a flexible ant colony (FAC) algorithm for solving protein folding problems (PFPs) based on the hydrophobic-polar (HP) square lattice model. Different from the previous ant algorithms for PFPs, the pheromones in the proposed algorithm are placed on the arcs connecting adjacent squares in the lattice. Such pheromone placement model is similar to the one used in the traveling salesmen problems (TSPs), where pheromones are released on the arcs connecting the cities. Moreover, the collaboration of effective heuristic and pheromone strategies greatly enhances the performance of the algorithm so that the algorithm can achieve good results without local search methods. By testing some benchmark two-dimensional hydrophobic-polar (2D-HP) protein sequences, the performance shows that the proposed algorithm is quite competitive compared with some other well-known methods for solving the same protein folding problems

Modeling the AgInSbTe Memristor

The AgInSbTe memristor shows gradual resistance tuning characteristics, which makes it a potential candidate to emulate biological plastic synapses. The working mechanism of the device is complex, and both intrinsic charge-trapping mechanism and extrinsic electrochemical metallization effect are confirmed in the AgInSbTe memristor. Mathematical model of the AgInSbTe memristor has not been given before. We propose the flux-voltage controlled memristor model. With piecewise linear approximation technique, we deliver the flux-voltage controlled memristor model of the AgInSbTe memristor based on the experiment data. Our model fits the data well. The flux-voltage controlled memristor model and the piecewise linear approximation method are also suitable for modeling other kinds of memristor devices based on experiment data

Lattice Boltzmann method for relativistic hydrodynamics: Issues on conservation law of particle number and discontinuities

In this paper, we aim to address several important issues about the recently developed lattice Boltzmann (LB) model for relativistic hydrodynamics [M. Mendoza et al., Phys. Rev. Lett. 105, 014502 (2010); Phys. Rev. D 82, 105008 (2010)]. First, we study the conservation law of particle number in the relativistic LB model. Through the Chapman-Enskog analysis, it is shown that in the relativistic LB model the conservation equation of particle number is a convection-diffusion equation rather than a continuity equation, which makes the evolution of particle number dependent on the relaxation time. Furthermore, we investigate the origin of the discontinuities appeared in the relativistic problems with high viscosities, which were reported in a recent study [D. Hupp et al., Phys. Rev. D 84, 125015 (2011)]. A multiple-relaxation-time (MRT) relativistic LB model is presented to examine the influences of different relaxation times on the discontinuities. Numerical experiments show the discontinuities can be eliminated by setting the relaxation time $\tau_e$ (related to the bulk viscosity) to be sufficiently smaller than the relaxation time $\tau_v$ (related to the shear viscosity). Meanwhile, it is found that the relaxation time $\tau_\varepsilon$, which has no effect on the conservation equations at the Navier-Stokes level, will affect the numerical accuracy of the relativistic LB model. Moreover, the accuracy of the relativistic LB model for simulating moderately relativistic problems is also investigated.Comment: 7 figure

Lattice Boltzmann modeling of multiphase flows at large density ratio with an improved pseudopotential model

Owing to its conceptual simplicity and computational efficiency, the pseudopotential multiphase lattice Boltzmann (LB) model has attracted significant attention since its emergence. In this work, we aim to extend the pseudopotential LB model to simulate multiphase flows at large density ratio and relatively high Reynolds number. First, based on our recent work [Li et al., Phys. Rev. E. 86, 016709 (2012)], an improved forcing scheme is proposed for the multiple-relaxation-time pseudopotential LB model in order to achieve thermodynamic consistency and large density ratio in the model. Next, through investigating the effects of the parameter a in the Carnahan-Starling equation of state, we find that the interface thickness is approximately proportional to 1/sqrt(a). Using a smaller a will lead to a wider interface thickness, which can reduce the spurious currents and enhance the numerical stability of the pseudopotential model at large density ratio. Furthermore, it is found that a lower liquid viscosity can be gained in the pseudopotential model by increasing the kinematic viscosity ratio between the vapor and liquid phases. The improved pseudopotential LB model is numerically validated via the simulations of stationary droplet and droplet oscillation. Using the improved model as well as the above treatments, numerical simulations of droplet splashing on a thin liquid film are conducted at a density ratio in excess of 500 with Reynolds numbers ranging from 40 to 1000. The dynamics of droplet splashing is correctly reproduced and the predicted spread radius is found to obey the power law reported in the literature.Comment: 9 figures, 2 tables, accepted by Physical Review E (in press

Model performance of downscaling 1999-2004 hydrometeorological fields to the upper Rio Grande basin using different forcing datasets

This study downscaled more than five years of data (1999-2004) for hydrometeorological fields over the upper Rio Grande basin (URGB) to a 4-km resolution using a regional model [fifth-generation Pennsylvania State University-National Center for Atmospheric Research (NCAR) Mesoscale Model (MM5, version 3)] and two forcing datasets that include National Centers for Environmental Prediction (NCEP)-NCAR reanalysis-1 (R1) and North America Regional Reanalysis (NARR) data. The long-term high-resolution simulation results show detailed patterns of hydroclimatological fields that are highly related to the characteristics of the regional terrain; the most important of these patterns are precipitation localization features caused by the complex topography. In comparison with station observational data, the downscaling processing, on whichever forcing field is used, generated more accurate surface temperature and humidity fields than the Eta Model and NARR data, although it still included marked errors, such as a negative (positive) bias toward the daily maximum (minimum) temperature and overestimated precipitation, especially in the cold season. Comparing the downscaling results forced by the NARR and R1 with both the gridded and station observational data shows that under the NARR forcing, the MM5 model produced generally better results for precipitation, temperature, and humidity than it did under the R1 forcing. These improvements were more apparent in winter and spring. During the warm season, although the use of NARR improved the precipitation estimates statistically at the regional (basin) scale, it substantially underestimated them over the southern upper Rio Grande basin, partly because the NARR forcing data exhibited warm and dry biases in the monsoon-active region during the simulation period and improper domain selection. Analyses also indicate that over mountainous regions, both the Climate Prediction Center's (CPC's) gridded (0.25°) and NARR forcings underestimate precipitation in comparison with station gauge data. © 2008 American Meteorological Society