Studying reactive processes with classical dynamics: rebinding dynamics in MbNO

A new surface-crossing algorithm suitable for describing bond-breaking and bond-forming processes in molecular dynamics simulations is presented. The method is formulated for two intersecting potential energy manifolds which dissociate to different adiabatic states. During simulations, crossings are detected by monitoring an energy criterion. If fulfilled, the two manifolds are mixed over a finite number of time steps, after which the system is propagated on the second adiabat and the crossing is carried out with probability one. The algorithm is extensively tested (almost 0.5 mu s of total simulation time) for the rebinding of NO to myoglobin. The unbound surface ((FeNO)-N-...) is represented using a standard force field, whereas the bound surface (Fe-NO) is described by an ab initio potential energy surface. The rebinding is found to be nonexponential in time, in agreement with experimental studies, and can be described using two time constants. Depending on the asymptotic energy separation between the manifolds, the short rebinding timescale is between 1 and 9 ps, whereas the longer timescale is about an order of magnitude larger. NO molecules which do not rebind within 1 ns are typically found in the Xenon-4 pocket, indicating the high affinity of NO to this region in the protein.

Bernoulli free-boundary problems in strip-like domains and a property of permanent waves on water of finite depth

We study weak solutions for a class of free-boundary problems which includes as a special case the classical problem of travelling gravity waves on water of finite depth. We show that such problems are equivalent to problems in fixed domains and study the regularity of their solutions. We also prove that in very general situations the free boundary is necessarily the graph of a function.

Attic pottery in early Classical Thessaly. A case study from Achaia Phthiotis

This paper attributes a previously unnoticed Attic black-figured lekythos kept in the Archaeological Museum of Volos to the Pholos Group (ca. 470 BC) and discusses its findspot in the peripheral Thessalian district of Achaia Phthiotis. Beyond an art-historical appreciation of the hastily-drawn chariot scene on this lekythos and a discussion of stylistic parallels, which include a lekythos in Prague, it is argued that such lekythoi were socially important for their shape and small size that made them easily transportable. The assumed scarcity of Attic pottery in Thessaly can be questioned given that a considerable amount of Attic pottery from Thessalian locations is mentioned only in passing or remains unpublished in museum storage. Small late black-figured lekythoi predominate amongst Attic pottery shapes in Thessaly. The popularity of such lekythoi can become indicative of human mobility across the landscape and the consumption of imported (grave) goods by social groups other than elites.

An improved finite element space for discontinuous pressures

We consider incompressible Stokes flow with an internal interface at which the pressure is discontinuous, as happens for example in problems involving surface tension. We assume that the mesh does not follow the interface, which makes classical interpolation spaces to yield suboptimal convergence rates (typically, the interpolation error in the L(2)(Omega)-norm is of order h(1/2)). We propose a modification of the P(1)-conforming space that accommodates discontinuities at the interface without introducing additional degrees of freedom or modifying the sparsity pattern of the linear system. The unknowns are the pressure values at the vertices of the mesh and the basis functions are computed locally at each element, so that the implementation of the proposed space into existing codes is straightforward. With this modification, numerical tests show that the interpolation order improves to O(h(3/2)). The new pressure space is implemented for the stable P(1)(+)/P(1) mini-element discretization, and for the stabilized equal-order P(1)/P(1) discretization. Assessment is carried out for Poiseuille flow with a forcing surface and for a static bubble. In all cases the proposed pressure space leads to improved convergence orders and to more accurate results than the standard P(1) space. In addition, two Navier-Stokes simulations with moving interfaces (Rayleigh-Taylor instability and merging bubbles) are reported to show that the proposed space is robust enough to carry out realistic simulations. (c) 2009 Elsevier B.V. All rights reserved.

Hyperbolic unit groups and quaternion algebras

We classify the quadratic extensions K = Q[root d] and the finite groups G for which the group ring o(K)[G] of G over the ring o(K) of integers of K has the property that the group U(1)(o(K)[G]) of units of augmentation 1 is hyperbolic. We also construct units in the Z-order H(o(K)) of the quaternion algebra H(K) = (-1, -1/K), when it is a division algebra.

Bicyclic units, Bass cyclic units and free groups

Let G be a finite group and ZG its integral group ring. We show that if alpha is a nontrivial bicyclic unit of ZG, then there are bicyclic units beta and gamma of different types, such that and are non-abelian free groups. In the case when G is non-abelian of order coprime to 6 we prove the existence of a bicyclic unit u and a Bass cyclic unit v in ZG, such that < u(m), v > is a free non-abelian group for all sufficiently large positive integers m.

REALIZING PROFINITE REDUCED SPECIAL GROUPS

Special groups are an axiomatization of the algebraic theory of quadratic forms over fields. It is known that any finite reduced special group is the special group of some field. We show that any special group that is the projective limit of a projective system of finite reduced special groups is also the special group of some field.

The lower central and derived series of the braid groups of the projective plane

In this paper, we determine the lower central and derived series for the braid groups of the projective plane. We are motivated in part by the study of Fadell-Neuwirth short exact sequences, but the problem is interesting in its own right. The n-string braid groups B(n)(RP(2)) of the projective plane RP(2) were originally studied by Van Buskirk during the 1960s. and are of particular interest due to the fact that they have torsion. The group B(1)(RP(2)) (resp. B(2)(RP(2))) is isomorphic to the cyclic group Z(2) of order 2 (resp. the generalised quaternion group of order 16) and hence their lower central and derived series are known. If n > 2, we first prove that the lower central series of B(n)(RP(2)) is constant from the commutator subgroup onwards. We observe that Gamma(2)(B(3)(RP(2))) is isomorphic to (F(3) X Q(8)) X Z(3), where F(k) denotes the free group of rank k, and Q(8) denotes the quaternion group of order 8, and that Gamma(2)(B(4)(RP(2))) is an extension of an index 2 subgroup K of P(4)(RP(2)) by Z(2) circle plus Z(2). As for the derived series of B(n)(RP(2)), we show that for all n >= 5, it is constant from the derived subgroup onwards. The group B(n)(RP(2)) being finite and soluble for n <= 2, the critical cases are n = 3, 4. We are able to determine completely the derived series of B(3)(RP(2)). The subgroups (B(3)(RP(2)))((1)), (B(3)(RP(2)))((2)) and (B(3)(RP(2)))((3)) are isomorphic respectively to (F(3) x Q(8)) x Z(3), F(3) X Q(8) and F(9) X Z(2), and we compute the derived series quotients of these groups. From (B(3)(RP(2)))((4)) onwards, the derived series of B(3)(RP(2)), as well as its successive derived series quotients, coincide with those of F(9). We analyse the derived series of B(4)(RP(2)) and its quotients up to (B(4)(RP(2)))((4)), and we show that (B(4)(RP(2)))((4)) is a semi-direct product of F(129) by F(17). Finally, we give a presentation of Gamma(2)(B(n)(RP(2))). (C) 2011 Elsevier Inc. All rights reserved.

Some classes of semisimple group (and loop) algebras over finite fields

We determine the structure of the semisimple group algebra of certain groups over the rationals and over those finite fields where the Wedderburn decompositions have the least number of simple components We apply our work to obtain similar information about the loop algebras of mdecomposable RA loops and to produce negative answers to the isomorphism problem over various fields (C) 2010 Elsevier Inc All rights reserved

CLASSIFICATION OF SUBALGEBRAS OF THE CAYLEY ALGEBRA OVER A FINITE FIELD

We classify all unital subalgebras of the Cayley algebra O(q) over the finite field F(q), q = p(n). We obtain the number of subalgebras of each type and prove that all isomorphic subalgebras are conjugate with respect to the automorphism group of O(q). We also determine the structure of the Moufang loops associated with each subalgebra of O(q).

On Groups Whose Maximal Cyclic Subgroups Are Maximal

We consider the problem of classifying those groups whose maximal cyclic subgroups are maximal. We give a complete classification of those groups with this property and which are either soluble or residually finite.

Classification of the virtually cyclic subgroups of the pure braid groups of the projective plane

We classify the ( finite and infinite) virtually cyclic subgroups of the pure braid groups P(n)(RP(2)) of the projective plane. The maximal finite subgroups of P(n)(RP(2)) are isomorphic to the quaternion group of order 8 if n = 3, and to Z(4) if n >= 4. Further, for all n >= 3, the following groups are, up to isomorphism, the infinite virtually cyclic subgroups of P(n)(RP(2)): Z, Z(2) x Z and the amalgamated product Z(4)*(Z2)Z(4).

Dynamical erosion of asteroid groups in the region of the Phocaea family

In a previous paper, the current state of knowledge of the region containing the Phocaea dynamical family was revised. Here, the dynamical evolution and possible origin of the Phocaea dynamical family and asteroid groups in the region are investigated. First, I study the case of asteroids at high eccentricity (e > 0.31). I find that these objects are unstable because of encounters with Mars on time-scales of up to 270 Myr. The minimum time needed by members of the Phocaea classical family to reach the orbital locations of these objects, 370 Myr, can be used to set a lower limit on the age of the Phocaea family.Next, attention is focused on the chaotic layer previously identified near the nu(6) secular resonance border. Using analytical and numerical tools, I find that the presence of the nu(6) secular resonance forces asteroids with vertical bar g-g(6)vertical bar < 2.55 arcsec yr(-1) to reach eccentricities high enough to allow them to experience deep, close encounters with Mars. Results of the analytical model of Yoshikawa and of my numerical simulations fully explain the low-inclination chaotic region found by Carruba.Finally, I investigate the long-term stability of the minor families and clumps identified in the previous paper, with particular emphasis on a clump only identifiable in the domain of proper frequencies (n, g, g - s) around (6246) Komurotoru. I find that while the clumps identified in the space of proper elements quickly disperse when the Yarkovsky effect is considered, the family around (19536) is still observable for time-scales of more than 50 Myr. The (6246) clump, characterized by its interaction with the nu(5) + nu(16) and 2 nu(6) - nu(16) secular resonances, is robust on time-scales of 50 Myr. I confirm that this group may be the first clump ever detected in the frequency domain that can be associated with a real collisional event.

Dynamical erosion of asteroid groups in the region of the Pallas family

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Finite Element Modeling for Development and Optimization of a Bone Plate for Mandibular Fracture in Dogs

This study aimed to develop a plate to treat fractures of the mandibular body in dogs and to validate the project using finite elements and biomechanical essays. Mandible prototypes were produced with 10 oblique ventrorostral fractures (favorable) and 10 oblique ventrocaudal fractures (unfavorable). Three groups were established for each fracture type. Osteosynthesis with a pure titanium plate of double-arch geometry and blocked monocortical screws offree angulanon were used. The mechanical resistance of the prototype with unfavorable fracture was lower than that of the fcworable fracture. In both fractures, the deflection increased and the relative stiffness decreased proportionally to the diminishing screw number The finite element analysis validated this plate study, since the maximum tension concentration observed on the plate was lower than the resistance limit tension admitted by the titanium. In conclusion, the double-arch geometry plate fixed with blocked monocortical screws has sufficient resistance to stabilize oblique,fractures, without compromising mandibular dental or neurovascular structures. J Vet Dent 24 (7); 212 - 221, 2010