The Mathematical Theory of Ito Diffusions on Hypersurfaces, with Applications to NMR Relaxation Problems

A mathematical framework for translational Brownian motion on hypersurfaces is presented, using an imbedding of the surface and Ito diffusions in the ambient space. This includes a survey of Ito calculus and differential geometry. Computational methods for time correlation functions relevant to spin relaxation studies on curved interfaces are given, and explicit calculations of time correlation functions and order parameters for a Rippled surface are presented.     

Molecular modeling of cytochrome P450 3A4

The three-dimensional structure of human cytochrome P450 3A4 was modeled based on crystallographic coordinates of four bacterial P450s: P450 BM-3, P450cam, P450terp, and P450eryF. The P450 3A4 sequence was aligned to those of the known proteins using a structure-based alignment of P450 BM-3, P450cam, P450terp, and P450eryF. The coordinates of the model were then calculated using a consensus strategy, and the final structure was optimized in the presence of water. The P450 3A4 model resembles P450 BM-3 the most, but the B helix is similar to that of P450eryF, which leads to an enlarged active site when compared with P450 BM-3, P450cam, and P450terp. The 3A4 residues equivalent to known substrate contact residues of the bacterial proteins and key residues of rat P450 2B1 are located in the active site or the substrate access channel. Docking of progesterone into the P450 3A4 model demonstrated that the substrate bound in a 6-orientation can interact with a number of active site residues, such as 114, 119, 301, 304, 305, 309, 370, 373, and 479, through hydrophobic interactions. The active site of the enzyme can also accommodate erythromycin, which, in addition to the residues listed for progesterone, also contacts residues 101, 104, 105, 214, 215, 217, 218, 374, and 478. The majority of 3A4 residues which interact with progesterone and/or erythromycin possess their equivalents in key residues of P450 2B enzymes, except for residues 297, 480 and 482, which do not contact either substrate in P450 3A4. The results from docking of progesterone and erythromycin into the enzyme model make it possible to pinpoint residues which may be important for 3A4 function and to target them for site-directed mutagenesis.     

Analytic study of power spectra of the tent maps near band-splitting transitions

Successive band-splitting transitions occur in the one-dimensional map x_i+1=g(x_i),i=0, 1, 2,... withg(x)=x, (0 x 1/2) –x +, (1/2 <x 1) as the parameter is changed from 2 to 1. The transition point fromN (=2ⁿ) bands to 2Nbands is given by=(2)^1/N (n=0, 1,2,...). The time-correlation function _i=x_ix₀/(x₀)²,x_i x_i–x_i is studied in terms of the eigenvalues and eigenfunctions of the Frobenius-Perron operator of the map. It is shown that, near the transition point=2, _i–[(10–42)/17] _i,0-[(102-8)/51]_i,1 + [(7 + 42)/17](–1)ⁱe^–yi, where2(–2) is the damping constant and vanishes at=2, representing the critical slowing-down. This critical phenomenon is in strong contrast to the topologically invariant quantities, such as the Lyapunov exponent, which do not exhibit any anomaly at=2. The asymptotic expression for _i has been obtained by deriving an analytic form of _i for a sequence of which accumulates to 2 from the above. Near the transition point=(2)^1/N, the damping constant of _i fori N is given by _N=2(^N-2)/N. Numerical calculation is also carried out for arbitrary a and is shown to be consistent with the analytic results.     

Dipolar relaxations in an epoxy-amine system

Dipolar relaxations in a reactive epoxy-amine system based on diglycidyl ether of bisphenol A with 4,4^′-diaminodiphenylsulfone were studied with the Havriliak-Negami function. The system was cured isothermally at 140 °C using simultaneous kinetic and microdielectric studies.The relaxation time was calculated from the frequency f_max of the peak of the loss factor ε^′′ versus frequency f. A linear relation exists between the logarithm of the relaxation time τ and the glass transition temperature T_g. Then, logτ follows the Di Benedetto equation revisited by Pascault and Williams allowing the prediction of the relaxation time τ during cure.The unrelaxed permittivity at high frequency ε_u, the relaxed or static permittivity ε_s and the skewness parameter β were found independent on the conversion or the curing time. The distribution parameter α decreases as curing time t increases. The difficulty of orientation of dipoles in the electric field due to vitrification is responsible of these behaviours.     

High-density properties of hard spheres within a modified Percus-Yevick theory: The role of thermodynamic consistency

Using an integral-equation approach based upon an approximation for the tail function, the equilibrium properties of a system of hard spheres are studied with special concern for the behavior in the region of close packing. The closure adopted is such that full, internal consistency is ensured in the thermodynamics of the model with respect to both the two zero-separation theorems as well as to the more standard virial and fluctuation routes to the equation of state. The scheme also makes use of the continuity properties of the tail function and of the cavity distribution function at contact. These properties are explictly tested in the low-density limit up to the fourth derivative. The theory generates an equilibrium branch bounded on the high-density side by a point corresponding to a packing fraction0.78, a value which closely matches Rogers' least upper bound for the densest packing of spheres. The pair structure of the fluid in the state of random close packing is also compared to the type of local order predicted by the theory at similar densities.     

Global optimization of univariate Lipschitz functions: II. New algorithms and computational comparison

We consider the following global optimization problems for a Lipschitz functionf implicitly defined on an interval [a, b]. Problem P: find a globally-optimal value off and a corresponding point; Problem Q: find a set of disjoint subintervals of [a, b] containing only points with a globally-optimal value and the union of which contains all globally optimal points. A two-phase algorithm is proposed for Problem P. In phase I, this algorithm obtains rapidly a solution which is often globally-optimal. Moreover, a sufficient condition onf for this to be the case is given. In phase II, the algorithm proves the-optimality of the solution obtained in phase I or finds a sequence of points of increasing value containing one with a globally-optimal value. The new algorithm is empirically compared (on twenty problems from the literature) with a best possible algorithm (for which the optimal value is assumed to be known), with a passive algorithm and with the algorithms of Evtushenko, Galperin, Shen and Zhu, Piyavskii, Timonov and Schoen. For small, the new algorithm requires only a few percent more function evaluations than the best possible one. An extended version of Piyavskii's algorithm is proposed for problem Q. A sufficient condition onf is given for the globally optimal points to be in one-to-one correspondance with the obtained intervals. This result is achieved for all twenty test problems.The research of the authors has been supported by AFOSR grants 0271 and 0066 to Rutgers University. Research of the second author has been also supported by NSERC grant GP0036426, FCAR grant 89EQ4144 and partially by AFOSR grant 0066. We thank Nicole Paradis for her help in drawing the figures.     

Convergence rates of a global optimization algorithm

Second order conditions for the (pseudo-) convexity of a function restricted to an affine subspace are obtained by extending those already known for functions on ⁿ . These results are then used to analyse the (pseudo-) convexity of potential functions of the type introduced by Karmarkar.This research was completed while the first author was on sabbatical leave at the Département d'Informatiques et de Recherche Opérationelle, Université de Montréal, and supported by NSERC (grant Q015807). This research was also supported by NSERC (grants A8312 and A5408) and la Coopération franco-québécoise (project 20-02-13).