BELIEF PROPAGATION FOR OPTIMAL EDGE COVER IN THE RANDOM COMPLETE GRAPH

We apply the objective method of Aldous to the problem of finding the minimum-cost edge cover of the complete graph with random independent and identically distributed edge costs. The limit, as the number of vertices goes to infinity, of the expected minimum cost for this problem is known via a combinatorial approach of Hessler and Wastlund. We provide a proof of this result using the machinery of the objective method and local weak convergence, which was used to prove the (2) limit of the random assignment problem. A proof via the objective method is useful because it provides us with more information on the nature of the edge's incident on a typical root in the minimum-cost edge cover. We further show that a belief propagation algorithm converges asymptotically to the optimal solution. This can be applied in a computational linguistics problem of semantic projection. The belief propagation algorithm yields a near optimal solution with lesser complexity than the known best algorithms designed for optimality in worst-case settings.

Minimization of Incompletely Specified Sequential Machines

A simple yet efficient method for the minimization of incompletely specified sequential machines (ISSMs) is proposed. Precise theorems are developed, as a consequence of which several compatibles can be deleted from consideration at the very first stage in the search for a minimal closed cover. Thus, the computational work is significantly reduced. Initial cardinality of the minimal closed cover is further reduced by a consideration of the maximal compatibles (MC's) only; as a result the method converges to the solution faster than the existing procedures. "Rank" of a compatible is defined. It is shown that ordering the compatibles, in accordance with their rank, reduces the number of comparisons to be made in the search for exclusion of compatibles. The new method is simple, systematic, and programmable. It does not involve any heuristics or intuitive procedures. For small- and medium-sized machines, it canle used for hand computation as well. For one of the illustrative examples used in this paper, 30 out of 40 compatibles can be ignored in accordance with the proposed rules and the remaining 10 compatibles only need be considered for obtaining a minimal solution.

SuSeFLAV 1.2: Program for supersymmetric mass spectra with seesaw mechanism and rare lepton flavor violating decays

Accurate supersymmetric spectra are required to confront data from direct and indirect searches of supersymmetry. SuSeFLAV is a numerical tool capable of computing supersymmetric spectra precisely for various supersymmetric breaking scenarios applicable even in the presence of flavor violation. The program solves MSSM RGEs with complete 3 x 3 flavor mixing at 2-loop level and one loop finite threshold corrections to all MSSM parameters by incorporating radiative electroweak symmetry breaking conditions. The program also incorporates the Type-I seesaw mechanism with three massive right handed neutrinos at user defined mass scales and mixing. It also computes branching ratios of flavor violating processes such as l(j) -> l(i)gamma, l(j) -> 3 l(i), b -> s gamma and supersymmetric contributions to flavor conserving quantities such as (g(mu) - 2). A large choice of executables suitable for various operations of the program are provided. Program summary Program title: SuSeFLAV Catalogue identifier: AEOD_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEOD_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License No. of lines in distributed program, including test data, etc.: 76552 No. of bytes in distributed program, including test data, etc.: 582787 Distribution format: tar.gz Programming language: Fortran 95. Computer: Personal Computer, Work-Station. Operating system: Linux, Unix. Classification: 11.6. Nature of problem: Determination of masses and mixing of supersymmetric particles within the context of MSSM with conserved R-parity with and without the presence of Type-I seesaw. Inter-generational mixing is considered while calculating the mass spectrum. Supersymmetry breaking parameters are taken as inputs at a high scale specified by the mechanism of supersymmetry breaking. RG equations including full inter-generational mixing are then used to evolve these parameters up to the electroweak breaking scale. The low energy supersymmetric spectrum is calculated at the scale where successful radiative electroweak symmetry breaking occurs. At weak scale standard model fermion masses, gauge couplings are determined including the supersymmetric radiative corrections. Once the spectrum is computed, the program proceeds to various lepton flavor violating observables (e.g., BR(mu -> e gamma), BR(tau -> mu gamma) etc.) at the weak scale. Solution method: Two loop RGEs with full 3 x 3 flavor mixing for all supersymmetry breaking parameters are used to compute the low energy supersymmetric mass spectrum. An adaptive step size Runge-Kutta method is used to solve the RGEs numerically between the high scale and the electroweak breaking scale. Iterative procedure is employed to get the consistent radiative electroweak symmetry breaking condition. The masses of the supersymmetric particles are computed at 1-loop order. The third generation SM particles and the gauge couplings are evaluated at the 1-loop order including supersymmetric corrections. A further iteration of the full program is employed such that the SM masses and couplings are consistent with the supersymmetric particle spectrum. Additional comments: Several executables are presented for the user. Running time: 0.2 s on a Intel(R) Core(TM) i5 CPU 650 with 3.20 GHz. (c) 2012 Elsevier B.V. All rights reserved.

Conformation of substituted polyacetylenes in solution : relationship between electronic properties and local order

A study of the chain conformation in solutions of polyphenylacetylene and poly(2-octyne) has been performed. The two polymers differ in many ways : polyphenylacetylene gives a red solution while poly(2-octyne) is transparent and, a marked difference on the chain rigidity is observed : the statistical length are 45 Å and 135 Å respectively. From the study of these two systems, one deduces that curvature fluctuations play a minor role on the π electrons localization, and that the torsion between monomer units is the pertinent parameter to understand the chain conformation and the π electrons localization.

Dissociation constants of weak organic acids in protic solvents obtained from their first hyperpolarizabilities in solution

The first hyperpolarizabilities (beta) of some weak aromatic organic acids have been measured in protic solvents by the hyper-Rayleigh scattering (HRS) technique at low concentrations. The measured hyperpolarizability (beta(m)) varies between the two extreme limits: the hyperpolarizability of the acid form (beta(HA)) at the lower side and that of the basic form (beta(A-)) at the higher side. The degree of dissociation (alpha) of the acid in a solvent is related to the measured hyperpolarizability, beta(m), by the following relationship: beta(m)(2)=(1-alpha)beta(HA)(2)+alpha beta(A-)(2). The calculated beta's including solvent effects in terms of an Onsager field do not reproduce the experimentally measured hyperpolarizabilities. Other solvent-induced effects like hydrogen bonding and van der Waals interactions seem to influence the first hyperpolarizability and, thus, indirectly the extent of dissociation of these weak acids in these protic solvents.

Spin Probe Signature of Freezing in Water Slow Cooling Experiment: Our Observation is Global or Local?

First systematic spin probe ESR study of water freezing has been conducted using TEMPOL and TEMPO as the probes. The spin probe signature of the water freezing has been described in terms of the collapse of narrow triplet spectrum into a single broad line. This spin probe signature of freezing has been observed at an anomalously low temperature when a milimoler solution of TEMPOL is slowly cooled from room temperature. A systematic observation has revealed a spin probe concentration dependence of these freezing and respective melting points. These results can be explained in terms of localization of spin probe and liquid water,most probably in the interstices of ice grains, in an ice matrix. The lowering of spin probe freezing point, along with the secondary evidences, like spin probe concentration dependence of peak-to-peak width in frozen limit signal, indicates a possible size dependence of these localizations/entrapments with spin probe concentration. A weak concentration dependence of spin probe assisted freezing and melting points, which has been observed for TEMPO in comparison to TEMPOL, indicates different natures of interactions with water of these two probes. This view is also supported by the relaxation behavior of the two probes.

Cavitation in materials with distributed weak zones: Implications on the origin of brittle fracture in metallic glasses

Several experimental studies have shown that fracture surfaces in brittle metallic glasses (MGs) generally exhibit nanoscale corrugations which may be attributed to the nucleation and coalescence of nanovoids during crack propagation. Recent atomistic simulations suggest that this phenomenon is due to large spatial fluctuations in material properties in a brittle MG, which leads to void nucleation in regions of low atomic density and then catastrophic fracture through void coalescence. To explain this behavior, we propose a model of a heterogeneous solid containing a distribution of weak zones to represent a brittle MG. Plane strain continuum finite element analysis of cavitation in such an elastic-plastic solid is performed with the weak zones idealized as periodically distributed regions having lower yield strength than the background material. It is found that the presence of weak zones can significantly reduce the critical hydrostatic stress for the onset of cavitation which is controlled uniquely by the local yield properties of these zones. Also, the presence of weak zones diminishes the sensitivity of the cavitation stress to the volume fraction of a preexisting void. These results provide plausible explanations for the observations reported in recent atomistic simulations of brittle MGs. An analytical solution for a composite, incompressible elastic-plastic solid with a weak inner core is used to investigate the effect of volume fraction and yield strength of the core on the nature of cavitation bifurcation. It is shown that snap-cavitation may occur, giving rise to sudden formation of voids with finite size, which does not happen in a homogeneous plastic solid. (c) 2012 Elsevier Ltd. All rights reserved.

Microscopic description of the evolution of the local structure and an evaluation of the chemical pressure concept in a solid solution

Extended x-ray absorption fine-structure studies have been performed at the Zn K and Cd K edges for a series of solid solutions of wurtzite Zn1-xCdxS samples with x = 0.0, 0.1, 0.25, 0.5, 0.75, and 1.0, where the lattice parameter as a function of x evolves according to the well-known Vegard's law. In conjunction with extensive, large-scale first-principles electronic structure calculations with full geometry optimizations, these results establish that the percentage variation in the nearest-neighbor bond distances are lower by nearly an order of magnitude compared to what would be expected on the basis of lattice parameter variation, seriously undermining the chemical pressure concept. With experimental results that allow us to probe up to the third coordination shell distances, we provide a direct description of how the local structure, apparently inconsistent with the global structure, evolves very rapidly with interatomic distances to become consistent with it. We show that the basic features of this structural evolution with the composition can be visualized with nearly invariant Zn-S-4 and Cd-S-4 tetrahedral units retaining their structural integrity, while the tilts between these tetrahedral building blocks change with composition to conform to the changing lattice parameters according to the Vegard's law within a relatively short length scale. These results underline the limits of applicability of the chemical pressure concept that has been a favored tool of experimentalists to control physical properties of a large variety of condensed matter systems.

Expansivity of liquid (cyclohexane + acetic anhydride) near its critical solution point

A simple volume dilatometer is described for the precise measurements of volume changes as a function of temperature in liquid mixtures. The expansivity of (cyclohexane + acetic anhydride) in the critical region was measured. The critical solution temperature Tc was approached to within 9 mK. For T > (Tc + 0.3 K), the results results follow both a logarithmic and a power-law behaviour with an exponent ≈ 1/8. But for T < (Tc + 0.3 K), the results seem to be affected possibly by gravity or temperature gradients. In this region, the expected expansivity anomaly is rounded off to a cusp. The expansivity shows a reduced anomaly for off-critical compositions. A discussion of the local extremum and a correlation between negative expansivity and the resistivity anomaly are also given.

X-ray and IR characterisation of Nasicon solid solution, Na1+xZr2SixP3−xO12

X-ray and IR studies on Nasicon solid solutions, Na1+xZr2SixP3−xO12, are carried out as a function of composition x. X-ray diffraction studies show that the unit cell volume increases as x increases and exhibits a maximum at x ≈ 2. On further increase in x the unit cell volume decreases. The infrared absorption peak positions and the splitting of these absorption peaks suggest a distortion of the PO4 and SiO4 tetrahedra. But the distortion is not large enough to change the local symmetry around the phosphorus or silicon ion from Td to C3v.

Solid state and solution conformation of Boc-L-Met-Aib-L-Phe-OMe. Beta-turn conformation of a sequence related to an active chemotactic peptide analog

The conformation of the peptide Boc-L-Met-Aib-L-Phe-OMe has been studied in the solid state and solution by X-ray diffraction and 1H n.m.r., respectively. The peptide differs only in the N-terminal protecting group from the biologically active chemotactic peptide analog formyl-L-Met-Aib-L-Phe-OMe. The molecules adopt a type-II beta-turn in the solid state with Met and Aib as the corner residues (phi Met = -51.8 degrees, psi Met = 139.5 degrees, phi Aib = 58.1 degrees, psi Aib = 37.0 degrees). A single, weak 4----1 intramolecular hydrogen bond is observed between the Boc CO and Phe NH groups (N---O 3.25 A, N-H---O 128.4 degrees). 1H n.m.r. studies, using solvent and temperature dependencies of NH chemical shifts and paramagnetic radical induced line broadening of NH resonances, suggest that the Phe NH is solvent shielded in CDCl3 and (CD3)2SO. Nuclear Overhauser effects observed between Met C alpha H and Aib NH protons provide evidence of the occurrence of Met-Aib type-II beta-turns in these solvents.

Weak discontinuities in relativisitic MHD

By using singular surface theory and ray theory the speeds of propagation of fast and slow waves, propagating into a medium in arbitrary motion, have been obtained in relativistic magnetohydrodynamics. The differential equation governing the growth of these waves along the rays has been derived and the solution has been presented in integral form.

R-parameter: A local truncation error based adaptive framework for finite volume compressible flow solvers

A residual-based strategy to estimate the local truncation error in a finite volume framework for steady compressible flows is proposed. This estimator, referred to as the -parameter, is derived from the imbalance arising from the use of an exact operator on the numerical solution for conservation laws. The behaviour of the residual estimator for linear and non-linear hyperbolic problems is systematically analysed. The relationship of the residual to the global error is also studied. The -parameter is used to derive a target length scale and consequently devise a suitable criterion for refinement/derefinement. This strategy, devoid of any user-defined parameters, is validated using two standard test cases involving smooth flows. A hybrid adaptive strategy based on both the error indicators and the -parameter, for flows involving shocks is also developed. Numerical studies on several compressible flow cases show that the adaptive algorithm performs excellently well in both two and three dimensions.

On the propagation of a weak shock front: Theory and application

In this paper the kinematics of a weak shock front governed by a hyperbolic system of conservation laws is studied. This is used to develop a method for solving problems, involving the propagation of nonlinear unimodal waves. It consists of first solving the nonlinear wave problem by moving along the bicharacteristics of the system and then fitting the shock into this solution field, so that it satisfies the necessary jump conditions. The kinematics of the shock leads in a natural way to the definition of ldquoshock-raysrdquo, which play the same role as the ldquoraysrdquo in a continuous flow. A special case of a circular cylinder introduced suddenly in a constant streaming flow is studied in detail. The shock fitted in the upstream region propagates with a velocity which is the mean of the velocities of the linear and the nonlinear wave fronts. In the downstream the solution is given by an expansion wave.

Intersoliton forces in weak-coupling quantum field theories

We offer a procedure for evaluating the forces exerted by solitons of weak-coupling field theories on one another. We illustrate the procedure for the kink and the antikink of the two-dimensional φ4 theory. To do this, we construct analytically a static solution of the theory which can be interpreted as a kink and an antikink held a distance R apart. This leads to a definition of the potential energy U(R) for the pair, which is seen to have all the expected features. A corresponding evaluation is also done for U(R) between a soliton and an antisoliton of the sine-Gordon theory. When this U(R) is inserted into a nonrelativistic two-body problem for the pair, it yields a set of bound states and phase shifts. These are found to agree with exact results known for the sine-Gordon field theory in those regions where U(R) is expected to be significant, i.e., when R is large compared to the soliton size. We take this agreement as support that our procedure for defining U(R) yields the correct description of the dynamics of well-separated soliton pairs. An important feature of U(R) is that it seems to give strong intersoliton forces when the coupling constant is small, as distinct from the forces between the ordinary quanta of the theory. We suggest that this is a general feature of a class of theories, and emphasize the possible relevance of this feature to real strongly interacting hadrons.