Factors contributing to phantom bonds in inorganic molecules: Interpretations and predictions

The validity of various qualitative proposals for interpreting and predicting the existence of short contacts between formally non-bonded atoms, as in cyclodisiloxane and related inorganic ring systems, is critically evaluated. The models range from simple considerations of geometric constraints, lone pair repulsions and pi-complex formation to proposals such as the unsupported pi-bond model and the sigma-bridged-pi bond concept. It is pointed out that a unified description based on a combination of closed and open 3-centre 2-electron bonds is possible. The role of hybridisation is emphasized in the short phantom bond computed in an earlier model system. These insights are used to predict structures with exceptionally short Si..Si and B..B phantom bonds. The proposals are confirmed by ab initio calculations.

Parameter estimation in water-distribution systems by least squares

The weighted-least-squares method using sensitivity-analysis technique is proposed for the estimation of parameters in water-distribution systems. The parameters considered are the Hazen-Williams coefficients for the pipes. The objective function used is the sum of the weighted squares of the differences between the computed and the observed values of the variables. The weighted-least-squares method can elegantly handle multiple loading conditions with mixed types of measurements such as heads and consumptions, different sets and number of measurements for each loading condition, and modifications in the network configuration due to inclusion or exclusion of some pipes affected by valve operations in each loading condition. Uncertainty in parameter estimates can also be obtained. The method is applied for the estimation of parameters in a metropolitan urban water-distribution system in India.

Fracture Analysis of Concrete Structural Components Accounting for Tension Softening Effect

This paper presents methodologies for fracture analysis of concrete structural components with and without considering tension softening effect. Stress intensity factor (SIF) is computed by using analytical approach and finite element analysis. In the analytical approach, SW accounting for tension softening effect has been obtained as the difference of SIP obtained using linear elastic fracture mechanics (LEFM) principles and SIP due to closing pressure. Superposition principle has been used by accounting for non-linearity in incremental form. SW due to crack closing force applied on the effective crack face inside the process zone has been computed using Green's function approach. In finite element analysis, the domain integral method has been used for computation of SIR The domain integral method is used to calculate the strain energy release rate and SIF when a crack grows. Numerical studies have been conducted on notched 3-point bending concrete specimen with and without considering the cohesive stresses. It is observed from the studies that SW obtained from the finite element analysis with and without considering the cohesive stresses is in good agreement with the corresponding analytical value. The effect of cohesive stress on SW decreases with increase of crack length. Further, studies have been conducted on geometrically similar structures and observed that (i) the effect of cohesive stress on SW is significant with increase of load for a particular crack length and (iii) SW values decreases with increase of tensile strength for a particular crack length and load.

Additive Spanners and (alpha, beta)-Spanners

An (alpha, beta)-spanner of an unweighted graph G is a subgraph H that distorts distances in G up to a multiplicative factor of a and an additive term beta. It is well known that any graph contains a (multiplicative) (2k - 1, 0)-spanner of size O(n(1+1/k)) and an (additive) (1, 2)-spanner of size O(n(3/2)). However no other additive spanners are known to exist. In this article we develop a couple of new techniques for constructing (alpha, beta)-spanners. Our first result is an additive (1, 6)-spanner of size O(n(4/3)). The construction algorithm can be understood as an economical agent that assigns costs and values to paths in the graph, purchasing affordable paths and ignoring expensive ones, which are intuitively well approximated by paths already purchased. We show that this path buying algorithm can be parameterized in different ways to yield other sparseness-distortion tradeoffs. Our second result addresses the problem of which (alpha, beta)-spanners can be computed efficiently, ideally in linear time. We show that, for any k, a (k, k - 1)-spanner with size O(kn(1+1/k)) can be found in linear time, and, further, that in a distributed network the algorithm terminates in a constant number of rounds. Previous spanner constructions with similar performance had roughly twice the multiplicative distortion.

Stretch effects extracted from inwardly and outwardly propagating spherical premixed flames

Unsteady propagation of spherical flames, both inward and outward, are studied numerically extensively for single-step reaction and for different Lewis numbers of fuel/oxidizer. The dependence of flame speed ratio (s) and flame temperature ratio are obtained for a range of Lewis numbers and stretch (kappa) values. These results of s versus kappa show that the asymptotic theory by Frankel and Sivashinsky is reasonable for outward propagation. Other theories are unsatisfactory both quantitatively and qualitatively. The stretch effects are much higher for negative stretch than for positive stretch, as also seen in the theory of Frankel and Sivashinsky. The linearity of the flame speed ratio vs stretch relationship is restricted to nondimensional stretch of +/-0.1. It is shown further that the results from cylindrical flames are identical to the spherical flame on flame speed ratio versus nondimensional stretch plot thus confirming the generality of the concept of stretch. The comparison of the variation of (ds/dkappa)kappa=0 with beta(Lc - 1) show an offset between the computed and the asymptotic results of Matalon and Matkowsky. The departure of negative stretch results from this variation is significant. Several earlier experimental results are analysed and set out in the form of s versus kappa plot. Comparison of the results with experiments seem reasonable for negative stretch. The results for positive stretch are satisfactory qualitatively for a few cases. For rich propane-air, there are qualitative differences pointing to the need for full chemistry calculations in the extraction of stretch effects.

A thermodynamic criterion for selection of gas compositions for diamond deposition

Isoactivity lines for carbon with respect to diamond as the standard state have been calculated in the ternary system C-H-O at 1223 K to identify the diamond deposition domain. The gas composition is calculated by suppressing the formation of all condensed forms of carbon using the SOLGASMIX free-energy minimization program. Thirty six gas species were included in the calculation. From the gas composition, isoactivity lines are computed using recent data on the Gibbs energy of diamond. Except for activities less than 0.1, the isoactivity lines are almost linear on the C-H-O ternary diagram. Gas compositions which generate activity of diamond ranging from 1 to 100 at 1223 K fall inside a narrow wedge originating from the point representing CO. This wedge is very similar to the revised lens-shaped diamond growth domain identified by Bachman et al., using inputs from experiment. The small difference between the calculated and observed domains may be attributed to variation in the supersaturation required for diamond deposition with gas composition. The diamond solubility in the gas phase along the isoactivity line for a(di)=100 and P=6.7 kPa exhibits a minimum at 1280 K, which is close to the optimum temperature found experimentally. At higher supersaturations, non-diamond forms of carbon, including amorphous varieties, are expected. The results suggest that thermodynamic calculations can be useful for locating diamond growth domains in more complex CVD systems containing halogens, for which very little experimental data is available.

Electron donor-acceptor complexes of I-2 with diethyl ether and diethyl sulphide. An ab initio MO study

The charge-transfer complexes of I-2 with the n-donors diethyl ether and diethyl sulfide were studied at the Hartree-Fock and MP2 levels. The structures were fully optimized using the 3-21G((*)) basis set as well as with effective core potentials. The calculations consistently yield a C-2v structure for the ether-I-2 complex, but an unsymmetrical form for the sulfide-I-2 complex. A natural bond orbital analysis and the BSSE-corrected complexation energies reveal stronger interactions in the sulfide complex. The computed orbital energies of the monomers and complexes reproduce the trends in experimentally observed vertical ionization potentials.

Electron donor-acceptor complex of ICl with diethyl ether - He I photoelectron spectroscopy and ab initio molecular orbital study

The He I photoelectron spectrum of the diethyl ether-ICl complex has been obtained. The oxygen orbitals are shifted to higher binding energies and that of ICl to lower binding energies owing to complex formation. Ab initio molecular orbital (MO) calculations of the complex molecule showed that the bonding is between the sigma-type lone pair of oxygen and the I atom and that the complex has C-2v symmetry. The binding energy of the complex is computed to be 8.06 kcal mol(-1) at the MP2/3-21G* level. The orbital energies obtained from the photoelectron spectra of the complex are compared and assigned with orbital energies obtained by MO calculations. Natural bond orbital analysis (NBO) shows that charge transfer is from the sigma-type oxygen lone pair to the iodine atom and the magnitude of charge transfer is 0.0744 e.

Automatic identification of modal damping from Floquet analysis

A practical method is proposed to identify the mode associated with the frequency part of the eigenvalue of the Floquet transition matrix (FTM). From the FTM eigenvector, which contains the states and their derivatives, the ratio of the derivative and the state corresponding to the largest component is computed. The method exploits the fact that the imaginary part of this (complex) ratio closely approximates the frequency of the mode. It also lends itself well to automation and has been tested over a large number of FTMs of order as high as 250.

Basis set dependence of the molecular electrostatic potential topography. A case study of substituted benzenes

The basis set dependence of the topographical structure of the molecular electrostatic potential (MESP), as well as the effect of substituents on the MESP distribution, has been investigated with substituted benzenes as test cases. The molecules are studied at HF-SCF 3�21G and 6�31G** levels, with a further MESP topographical investigation at the 3�21G, double-zeta, 6�31G*, 6�31G**, double-zeta polarized and triple-zeta polarized levels. The MESP critical points for a 3�21G optimized/6�31G** basis are similar to the corresponding 6�31G** optimized/6�31G** ones. More generally, the qualitative features of the MESP topography computed at the polarized level are independent of the level at which optimization is carried out. For a proper representation of oxygen lone pairs, however, optimization using a polarized basis set is required. The nature of the substituent drastically changes the MESP distribution over the phenyl ring. The values and positions of MESP minima indicate the most active site for electrophilic attack. This point is strengthened by a study of disubstituted benzenes.

Calculation of elastic pinning strength of semicoherent Y2BaCuO5 precipitates in YBa2Cu3O7

The pinning energy due to the elastic interaction of a semicoherent Y2BaCuO5 precipitate with the YBa2Cu3O7 matrix is computed. This is achieved by setting up dislocation arrays at the interface. The elastic stresses generated by such arrays are integrated over a fluxoid volume to obtain the energy. It is seen that this elastic interaction energy makes an additive contribution to the total J(c) value.

Hydrography and circulation in the western Bay of Bengal during the northeast monsoon

The Bay of Bengal, a semienclosed tropical basin that comes under the influence of monsoonal wind and freshwater influx, is distinguished by a strongly stratified surface layer and a seasonally reversing circulation. We discuss characteristics of these features in the western Bay during the northeast monsoon, when the East India Coastal Current (EICC) flows southward, using hydrographic data collected during December 1991. Vertical profiles show uniform temperature and salinity in a homogeneous surface layer, on average, 25 m deep but shallower northward and coastward. The halocline, immediately below, is approximately 50 m thick; salinity changes by approximately 3 parts per thousand. About two thirds of the profiles show temperature inversions in this layer. Salinity below the halocline hardly changes, and stratification is predominantly due to temperature variation, The halocline is noticeably better developed and the surface homogeneous layer is thinner in a low-salinity plume that hugs the coastline along the entire east coast of India, The plume is, on average, 50 km wide, with isohalines sloping down toward the coast. Most prominent in the geostrophic velocity field is the equatorward EICC. Its transport north of about 13 degrees N, computed with 1000 dbar as the level of reference, varies between 2.6 and 7.1 x 10(6) m(3) s(-1); just south of this latitude, a northwestward flow from offshore recurves and merges with the coastal current. At the southern end of the region surveyed, the transport is 7.7 x 10(6) m(3) s(-1). Recent model studies lead us to conclude that the EICC during the northeast monsoon is driven by winds along the east coast of India and Ekman pumping in the interior bay. In the south, Ekman pumping over the southwestern bay is responsible for the northwestward flow that merges with the EICC.

Image analysis using multifractals

The efficacy of the multifractal spectrum as a tool for characterizing images has been studied. This spectrum has been computed for digitized images of the nucleus of human cervical cancer cells and it was observed that the entire spectrum is almost fully reproduced for a normal cell while only the right half (q<0) of the spectrum is reproduced for a cancerous cell. Cells in stages in between the two extremes show a shortening of the left half of the spectrum proportional to their condition. The extent of this shortening has been found to be sufficient to permit a classification between three classes of cells at varying distances from a basal cancerous layer-the superficial cells, the intermediate cells and the parabasal cells. This technique may be used for automatic screening of the population while also indicating the stage of malignancy

Parallel computing concepts and methods for floquet analysis of helicopter trim and stability

Floquet analysis is widely used for small-order systems (say, order M < 100) to find trim results of control inputs and periodic responses, and stability results of damping levels and frequencies, Presently, however, it is practical neither for design applications nor for comprehensive analysis models that lead to large systems (M > 100); the run time on a sequential computer is simply prohibitive, Accordingly, a massively parallel Floquet analysis is developed with emphasis on large systems, and it is implemented on two SIMD or single-instruction, multiple-data computers with 4096 and 8192 processors, The focus of this development is a parallel shooting method with damped Newton iteration to generate trim results; the Floquet transition matrix (FTM) comes out as a byproduct, The eigenvalues and eigenvectors of the FTM are computed by a parallel QR method, and thereby stability results are generated, For illustration, flap and flap-lag stability of isolated rotors are treated by the parallel analysis and by a corresponding sequential analysis with the conventional shooting and QR methods; linear quasisteady airfoil aerodynamics and a finite-state three-dimensional wake model are used, Computational reliability is quantified by the condition numbers of the Jacobian matrices in Newton iteration, the condition numbers of the eigenvalues and the residual errors of the eigenpairs, and reliability figures are comparable in both the parallel and sequential analyses, Compared to the sequential analysis, the parallel analysis reduces the run time of large systems dramatically, and the reduction increases with increasing system order; this finding offers considerable promise for design and comprehensive-analysis applications.

Phase equilibria in the system CaO-CoO-SiO2 and Gibbs energies of formation of the quaternary oxides CaCoSi2O6, Ca2CoSi2O7, and CaCoSiO4

Phase relations in the pseudoternary system CaO-CoO-SiO2 have been established at 1323 K. Three quaternary oxides were found to be stable: CaCoSi2O6 with clinopyroxene (Cpx), Ca2CoSi2O7 with melilite (Mel), and CaCoSiO4 with olivine (Ol) structures. The Gibbs energies of formation of the quaternary oxides from their component binary oxides were measured using solid-state galvanic cells incorporating yttria-stabilized zirconia as the solid electrolyte in the temperature range of 1000-1324 K. The results can be summarized as follows: CoO (rs) + CaO (rs) + 2SiO(2) (Qtz) --> CaCoSi2O6 (Cpx), Delta G(f)(0) = -117920 + 11.26T (+/-150) J/mol CoO (rs) + 2CaO (rs) + 2SiO(2) (Qtz) --> Ca2CoSi2O7 (Mel), Delta G(f)(0) = -192690 + 2.38T (+/-130) J/mol CoO (rs) + CaO (rs) + SiO2 (Qtz) --> CaCoSiO2 (Ol), Delta G(f)(0) = -100325 + 2.55T (+/-100) J/mol where rs = rock salt (NaCl) structure and Qtz = quartz. The uncertainty limits correspond to twice the standard error estimate. The experimentally observed miscibility gaps along the joins CaO-CoO and CaCoSiO4-Co2SiO4 were used to calculate the excess free energies of mixing for the solid solutions CaxCo1-xO and (CayCo1-y)CoSiO4:Delta G(E) = X(1 - X)[31975X + 26736 (1 - X)] J/mol and Delta G(E) = 23100 (+/-250) Y(1 - Y) J/mol. A T-X phase diagram for the binary CaO-CoO was computed from the thermodynamic information; the diagram agrees with information available in the literature. The computed miscibility gap along the CaCoSiO4-Co2SiO4 join is associated with a critical temperature of 1389 (+/-15) K. Stability fields for the various solid solutions and the quaternary compounds are depicted on chemical-potential diagrams for SiO2, CaO, and CoO at 1323 K.