Symmetry properties of one- and two- electron molecular integrals

The maximum numbers of distinct one- and two-electron integrals that arise in calculating the electronic energy of a molecule are discussed. It is shown that these may be calculated easily using the character table of the symmetry group of the set of basis functions used to express the wave function. Complications arising from complex group representations and from a conflict of symmetry between the basis set and the nuclear configuration are considered and illustrated by examples.

Fast assembly of fock matrices utilising symmetry properties of the basis set

A method of assembling the elements of the Fock matrix is described which is a modification of that due to Dacre. Lists of symmetry equivalent one-electron integrals are used as pointers to abbreviate the process of collecting two-electron integrals into the Fock matrix.

Symmetry-breaking in chiral polymerisation

We propose a model for chiral polymerisation and investigate its symmetric and asymmetric solutions. The model has a source species which decays into left- and right-handed types of monomer, each of which can polymerise to form homochiral chains; these chains are susceptible to `poisoning' by the opposite handed monomer. Homochiral polymers are assumed to influence the proportion of each type of monomer formed from the precursor. We show that for certain parameter values a positive feedback mechanism makes the symmetric steady-state solution unstable. The kinetics of polymer formation are then analysed in the case where the system starts from zero concentrations of monomers and chains. We show that following a long induction time, extremely large concentrations of polymers are formed for a short time, during this time an asymmetry introduced into the system by a random external perturbation may be massively amplified. The system then approaches one of the steady-state solutions described above.

Adaptivity and a posteriori error control for bifurcation problems II: Incompressible fluid flow in open systems with Z_2 symmetry

In this article we consider the a posteriori error estimation and adaptive mesh refinement of discontinuous Galerkin finite element approximations of the bifurcation problem associated with the steady incompressible Navier-Stokes equations. Particular attention is given to the reliable error estimation of the critical Reynolds number at which a steady pitchfork or Hopf bifurcation occurs when the underlying physical system possesses reflectional or Z_2 symmetry. Here, computable a posteriori error bounds are derived based on employing the generalization of the standard Dual-Weighted-Residual approach, originally developed for the estimation of target functionals of the solution, to bifurcation problems. Numerical experiments highlighting the practical performance of the proposed a posteriori error indicator on adaptively refined computational meshes are presented.