Escape times for rigid Brownian rotators in a bistable potential from the time evolution of the Green function and the characteristic time of the probability evolution

The greatest relaxation time for an assembly of three- dimensional rigid rotators in an axially symmetric bistable potential is obtained exactly in terms of continued fractions as a sum of the zero frequency decay functions (averages of the Legendre polynomials) of the system. This is accomplished by studying the entire time evolution of the Green function (transition probability) by expanding the time dependent distribution as a Fourier series and proceeding to the zero frequency limit of the Laplace transform of that distribution. The procedure is entirely analogous to the calculation of the characteristic time of the probability evolution (the integral of the configuration space probability density function with respect to the position co-ordinate) for a particle undergoing translational diffusion in a potential; a concept originally used by Malakhov and Pankratov (Physica A 229 (1996) 109). This procedure allowed them to obtain exact solutions of the Kramers one-dimensional translational escape rate problem for piecewise parabolic potentials. The solution was accomplished by posing the problem in terms of the appropriate Sturm-Liouville equation which could be solved in terms of the parabolic cylinder functions. The method (as applied to rotational problems and posed in terms of recurrence relations for the decay functions, i.e., the Brinkman approach c.f. Blomberg, Physica A 86 (1977) 49, as opposed to the Sturm-Liouville one) demonstrates clearly that the greatest relaxation time unlike the integral relaxation time which is governed by a single decay function (albeit coupled to all the others in non-linear fashion via the underlying recurrence relation) is governed by a sum of decay functions. The method is easily generalized to multidimensional state spaces by matrix continued fraction methods allowing one to treat non-axially symmetric potentials, where the distribution function is governed by two state variables. (C) 2001 Elsevier Science B.V. All rights reserved.

Maximum likelihood estimation of higher-order integer-valued autoregressive processes

In this article, we extend the earlier work of Freeland and McCabe [Journal of time Series Analysis (2004) Vol. 25, pp. 701–722] and develop a general framework for maximum likelihood (ML) analysis of higher-order integer-valued autoregressive processes. Our exposition includes the case where the innovation sequence has a Poisson distribution and the thinning is binomial. A recursive representation of the transition probability of the model is proposed. Based on this transition probability, we derive expressions for the score function and the Fisher information matrix, which form the basis for ML estimation and inference. Similar to the results in Freeland and McCabe (2004), we show that the score function and the Fisher information matrix can be neatly represented as conditional expectations. Using the INAR(2) speci?cation with binomial thinning and Poisson innovations, we examine both the asymptotic e?ciency and ?nite sample properties of the ML estimator in relation to the widely used conditional least
squares (CLS) and Yule–Walker (YW) estimators. We conclude that, if the Poisson assumption can be justi?ed, there are substantial gains to be had from using ML especially when the thinning parameters are large.

Computerised Diagnostic Decision Support System for the Classification of Pre-Invasive Cervical Squamous Lesions.

Previous studies have revealed considerable interobserver and intraobserver variation in the histological classification of preinvasive cervical squamous lesions. The aim of the present study was to develop a decision support system (DSS) for the histological interpretation of these lesions. Knowledge and uncertainty were represented in the form of a Bayesian belief network that permitted the storage of diagnostic knowledge and, for a given case, the collection of evidence in a cumulative manner that provided a final probability for the possible diagnostic outcomes. The network comprised 8 diagnostic histological features (evidence nodes) that were each independently linked to the diagnosis (decision node) by a conditional probability matrix. Diagnostic outcomes comprised normal; koilocytosis; and cervical intraepithelial neoplasia (CIN) 1, CIN II, and CIN M. For each evidence feature, a set of images was recorded that represented the full spectrum of change for that feature. The system was designed to be interactive in that the histopathologist was prompted to enter evidence into the network via a specifically designed graphical user interface (i-Path Diagnostics, Belfast, Northern Ireland). Membership functions were used to derive the relative likelihoods for the alternative feature outcomes, the likelihood vector was entered into the network, and the updated diagnostic belief was computed for the diagnostic outcomes and displayed. A cumulative probability graph was generated throughout the diagnostic process and presented on screen. The network was tested on 50 cervical colposcopic biopsy specimens, comprising 10 cases each of normal, koilocytosis, CIN 1, CIN H, and CIN III. These had been preselected by a consultant gynecological pathologist. Using conventional morphological assessment, the cases were classified on 2 separate occasions by 2 consultant and 2 junior pathologists. The cases were also then classified using the DSS on 2 occasions by the 4 pathologists and by 2 medical students with no experience in cervical histology. Interobserver and intraobserver agreement using morphology and using the DSS was calculated with K statistics. Intraobserver reproducibility using conventional unaided diagnosis was reasonably good (kappa range, 0.688 to 0.861), but interobserver agreement was poor (kappa range, 0.347 to 0.747). Using the DSS improved overall reproducibility between individuals. Using the DSS, however, did not enhance the diagnostic performance of junior pathologists when comparing their DSS-based diagnosis against an experienced consultant. However, the generation of a cumulative probability graph also allowed a comparison of individual performance, how individual features were assessed in the same case, and how this contributed to diagnostic disagreement between individuals. Diagnostic features such as nuclear pleomorphism were shown to be particularly problematic and poorly reproducible. DSSs such as this therefore not only have a role to play in enhancing decision making but also in the study of diagnostic protocol, education, self-assessment, and quality control. (C) 2003 Elsevier Inc. All rights reserved.

Stochastic Sznajd Model in open community

We extend the Sznajd Model for opinion formation by introducing persuasion probabilities for opinions. Moreover, we couple the system to an environment which mimics the application of the opinion. This results in a feedback, representing single-state opinion transitions in opposite to the two-state opinion transitions for persuading other people. We call this model opinion formation in an open community (OFOC). It can be seen as "stochastic extension of the Sznajd model for an open community, because it allows for "special choice of parameters to recover the original Sznajd model. We demonstrate the effect of feedback in the OFOC model by applying it to a scenario in which, e.g., opinion B is worse then opinion A but easier explained to other people. Casually formulated we analyzed the question, how much better one has to be, in order to persuade other people, provided the opinion is worse. Our results reveal a linear relation between the transition probability for opinion B and the influence of the environment on B.

PROSTATIC INTRAEPITHELIAL NEOPLASIA (PIN) - PERFORMANCE OF BAYESIAN BELIEF NETWORK FOR DIAGNOSIS AND GRADING

Prostatic intraepithelial neoplasia (PIN) diagnosis and grading are affected by uncertainties which arise from the fact that almost all knowledge of PIN histopathology is expressed in concepts, descriptive linguistic terms, and words. A Bayesian belief network (BBN) was therefore used to reduce the problem of uncertainty in diagnostic clue assessment, while still considering the dependences between elements in the reasoning sequence. A shallow network was used with an open-tree topology, with eight first-level descendant nodes for the diagnostic clues (evidence nodes), each independently linked by a conditional probability matrix to a root node containing the diagnostic alternatives (decision node). One of the evidence nodes was based on the tissue architecture and the others were based on cell features. The system was designed to be interactive, in that the histopathologist entered evidence into the network in the form of likelihood ratios for outcomes at each evidence node. The efficiency of the network was tested on a series of 110 prostate specimens, subdivided as follows: 22 cases of non-neoplastic prostate or benign prostatic tissue (NP), 22 PINs of low grade (PINlow), 22 PINs of high grade (PINhigh), 22 prostatic adenocarcinomas with cribriform pattern (PACcri), and 22 prostatic adenocarcinomas with large acinar pattern (PAClgac). The results obtained in the benign and malignant categories showed that the belief for the diagnostic alternatives is very high, the values being in general more than 0.8 and often close to 1.0. When considering the PIN lesions, the network classified and graded most of the cases with high certainty. However, there were some cases which showed values less than 0.8 (13 cases out of 44), thus indicating that there are situations in which the feature changes are intermediate between contiguous categories or grades. Discrepancy between morphological grading and the BBN results was observed in four out of 44 PIN cases: one PINlow was classified as PINhigh and three PINhigh were classified as PINlow. In conclusion, the network can grade PlN lesions and differentiate them from other prostate lesions with certainty. In particular, it offers a descriptive classifier which is readily implemented and which allows the use of linguistic, fuzzy variables.

EXPERT-SYSTEM SUPPORT USING BAYESIAN BELIEF NETWORKS IN THE DIAGNOSIS OF FINE-NEEDLE ASPIRATION BIOPSY SPECIMENS OF THE BREAST

Aim-To develop an expert system model for the diagnosis of fine needle aspiration cytology (FNAC) of the breast.

Methods-Knowledge and uncertainty were represented in the form of a Bayesian belief network which permitted the combination of diagnostic evidence in a cumulative manner and provided a final probability for the possible diagnostic outcomes. The network comprised 10 cytological features (evidence nodes), each independently linked to the diagnosis (decision node) by a conditional probability matrix. The system was designed to be interactive in that the cytopathologist entered evidence into the network in the form of likelihood ratios for the outcomes at each evidence node.

Results-The efficiency of the network was tested on a series of 40 breast FNAC specimens. The highest diagnostic probability provided by the network agreed with the cytopathologists' diagnosis in 100% of cases for the assessment of discrete, benign, and malignant aspirates. A typical probably benign cases were given probabilities in favour of a benign diagnosis. Suspicious cases tended to have similar probabilities for both diagnostic outcomes and so, correctly, could not be assigned as benign or malignant. A closer examination of cumulative belief graphs for the diagnostic sequence of each case provided insight into the diagnostic process, and quantitative data which improved the identification of suspicious cases.

Conclusion-The further development of such a system will have three important roles in breast cytodiagnosis: (1) to aid the cytologist in making a more consistent and objective diagnosis; (2) to provide a teaching tool on breast cytological diagnosis for the non-expert; and (3) it is the first stage in the development of a system capable of automated diagnosis through the use of expert system machine vision.

UV VIS AND RESONANCE RAMAN SPECTROELECTROCHEMICAL PROPERTIES OF TRANSITION-METAL CENTERS IMMOBILIZED WITHIN A POLY(AMINO ACID) MATRIX - ILLUSTRATED WITH AN IRON PORPHYRIN

An iron prophyrin complex has been immobilized on the surfaces of platinum, silver, and indium doped-tin oxide coated glass by using the poly(gamma-ethyl L-glutamate)-N-(3-aminopropyl)imidazole derivative 1 as a linking agent, thus allowing-the surface-enhanced resonance Raman and UV-VIS absorption spectra and electrochemical properties of the porphyrin to be studied in solvents in which it is not normally soluble.

R-matrix Floquet theory of multiphoton processes:XIII. Resonances in H multiphoton detachment

The results of calculations investigating the effects of autodetaching resonances on the multiphoton detachment spectra of H are presented. The R-matrix Floquet method is used, in which the coupling of the ion with the laser field is described non-perturbatively. The laser field is fixed at an intensity of 10 W cm, while frequency ranges are chosen such that the lowest autodetaching states of the ion are excited through a two- or three-photon transition from the ground state. Detachment rates are compared, where possible, to previous results obtained using perturbation theory. An illustration of how non-lowest-order processes, involving autodetaching states, can lead to light-induced continuum structures is also presented. Finally, it is demonstrated that by using a frequency connecting the 1s and 2s states, the probability of exciting the residual hydrogen atom is significantly enhanced.

The near-edge structure in energy-loss spectroscopy: many-electron and magnetic effects in transition metal nitrides and carbides

We investigate the ability of the local density approximation (LDA) in density functional theory to predict the near-edge structure in electron energy-loss spectroscopy in the dipole approximation. We include screening of the core hole within the LDA using Slater's transition state theory. We find that anion K-edge threshold energies are systematically overestimated by 4.22 +/- 0.44 eV in twelve transition metal carbides and nitrides in the rock-salt (B1) structure. When we apply this 'universal' many-electron correction to energy-loss spectra calculated within the transition state approximation to LDA, we find quantitative agreement with experiment to within one or two eV for TiC, TiN and VN. We compare our calculations to a simpler approach using a projected Mulliken density which honours the dipole selection rule, in place of the dipole matrix element itself. We find remarkably close agreement between these two approaches. Finally, we show an anomaly in the near-edge structure in CrN to be due to magnetic structure. In particular, we find that the N K edge in fact probes the magnetic moments and alignments of ther sublattice.

An exactly solvable four transition point model

Heavy particle collisions, in particular low-energy ion-atom collisions, are amenable to semiclassical JWKB phase integral analysis in the complex plane of the internuclear separation. Analytic continuation in this plane requires due attention to the Stokes phenomenon which parametrizes the physical mechanisms of curve crossing, non-crossing, the hybrid Nikitin model, rotational coupling and predissociation. Complex transition points represent adiabatic degeneracies. In the case of two or more such points, the Stokes constants may only be completely determined by resort to the so-called comparison- equation method involving, in particular, parabolic cylinder functions or Whittaker functions and their strong-coupling asymptotics. In particular, the Nikitin model is a two transition-point one-double-pole problem in each half-plane corresponding to either ingoing or outgoing waves. When the four transition points are closely clustered, new techniques are required to determine Stokes constants. However, such investigations remain incomplete, A model problem is therefore solved exactly for scattering along a one-dimensional z-axis. The energy eigenvalue is b(2)-a(2) and the potential comprises -z(2)/2 (parabolic) and -a(2) + b(2)/2z(2) (centrifugal/centripetal) components. The square of the wavenumber has in the complex z-plane, four zeros each a transition point at z = +/-a +/- ib and has a double pole at z = 0. In cases (a) and (b), a and b are real and unitarity obtains. In case (a) the reflection and transition coefficients are parametrized by exponentials when a(2) + b(2) > 1/2. In case (b) they are parametrized by trigonometrics when a(2) + b(2) <1/2 and total reflection is achievable. In case (c) a and b are complex and in general unitarity is not achieved due to loss of flux to a continuum (O'Rourke and Crothers, 1992 Proc. R. Sec. 438 1). Nevertheless, case (c) coefficients reduce to (a) or (b) under appropriate limiting conditions. Setting z = ht, with h a real constant, an attempt is made to model a two-state collision problem modelled by a pair of coupled first-order impact parameter equations and an appropriate (T) over tilde-tau relation, where (T) over tilde is the Stueckelberg variable and tau is the reduced or scaled time. The attempt fails because (T) over tilde is an odd function of tau, which is unphysical in a real collision problem. However, it is pointed out that by applying the Kummer exponential model to each half-plane (O'Rourke and Crothers 1994 J. Phys. B: At. Mel. Opt. Phys. 27 2497) the current model is in effect extended to a collision problem with four transition points and a double pole in each half-plane. Moreover, the attempt in itself is not a complete failure since it is shown that the result is a perfect diabatic inelastic collision for a traceless Hamiltonian matrix, or at least when both diagonal elements are odd and the off-diagonal elements equal and even.

Breit-Pauli R-matrix calculation of fine-structure effective collision strengths for the electron impact excitation of Mg V

Context. Electron-impact excitation collision strengths are required for the analysis and interpretation of stellar observations.
Aims. This calculation aims to provide effective collision strengths for the Mg V ion for a larger number of transitions and for a greater temperature range than previously available, using collision strength data that include contributions from resonances.
Methods. A 19-state Breit-Pauli R-matrix calculation was performed. The target states are represented by configuration interaction wavefunctions and consist of the 19 lowest LS states, having configurations 2s22p4, 2s2p5, 2p6, 2s22p33s, and 2s22p33p. These target states give rise to 37 fine-structure levels and 666 possible transitions. The effective collision strengths were calculated by averaging the electron collision strengths over a Maxwellian distribution of electron velocities.
Results. The non-zero effective collision strengths for transitions between the fine-structure levels are given for electron temperatures in the range = 3.0 - 7.0. Data for transitions among the 5 fine-structure levels arising from the 2s22p4 ground state configurations, seen in the UV range, are discussed in the paper, along with transitions in the EUV range – transitions from the ground state 3P levels to 2s2p5?3P levels. The 2s22p4?1D–2s2p5?1P transition is also noted. Data for the remaining transitions are available at the CDS.

Ultra-fast spin dynamics in the 2s2p^2 configuration of C^+ studied by time-dependent R-matrix theory

We employ time-dependent R-matrix theory to study ultra-fast dynamics in the doublet 2s2p(2) configuration of C+ for a total magnetic quantum number M = 1. In contrast to the dynamics observed for M = 0, ultra-fast dynamics for M = 1 is governed by spin dynamics in which the 2s electron acts as a flag rather than a spectator electron. Under the assumption that m(S) = 1/2, m(2s) = 1/2 allows spin dynamics involving the two 2p electrons, whereas m(2s) = -1/2 prevents spin dynamics of the two 2p electrons. For a pump-probe pulse scheme with (h) over bar omega(pump) = 10.9 eV and (h) over bar omega(probe) = 16.3 eV and both pulses six cycles long, little sign of spin dynamics is observed in the total ionization probability. Signs of spin dynamics can be observed, however, in the ejected-electron momentum distributions. We demonstrate that the ejected-electron momentum distributions can be used for unaligned targets to separate the contributions of initial M = 0 and M = 1 levels. This would, in principle, allow unaligned target ions to be used to obtain information on the different dynamics in the 2s2p(2) configuration for the M = 0 and M = 1 levels from a single experime

Numerical study of the effects of reinforcement/matrix interphase on stress-strain behavior of YAl2 particle reinforced MgLiAl composites

For the potential influence produced by the reinforcement/matrix interphase in particle reinforced metal matrix composites (PMMCs), a unit cell model with transition interphase was proposed. Uniaxial tensile loading was simulated and the stress/strain behavior was predicted. The results show that a transition interphase with both appropriate strength and thickness could affect the failure mode, reduce the stress concentration, and enhance the maximum strain value of the composite.