### Submissions

##### Asymptotic solution of natural convection in a uniformly Joule-heating shallow cavity

xiaohui zhang, Soochow University, , China

Abstract: The steady laminar two-dimensional Joule heating natural convection is investigated using asymptotical analysis, the fluid is in a rectangular cavity, the direct current contributes heat for heating the process medium by a pair of plate electrodes, the top wall is cooled with atmosphere and all the other walls are kept thermally insulated. The asymptotic solution is obtained in the core region in the limit as the aspect ratio, which is defined as the ratio of the vertical dimension of cavity to the horizontal dimension of cavity, goes to zero. The numerical experiments are also carried out to compare with the asymptotic solution of the steady two-dimensional Joule heating convection. The asymptotic results indicate that the expressions of velocity and temperature fields in the core region are valid in the limit of the small aspect ratio.
##### Ab-initio simulation of the ionization and fragmentation of acetylene by strong femtosecond laser pulses

Kalman Varga, Vanderbilt, Physics, United States

Abstract: The electron and nuclear dynamics of acetylene when interacting with strong short laser pulses has been simulated in the framework of real–space Time Dependent Density Functional Theory (TDDFT) and molecular dynamics. The stretching and dissociation of individual bonds are reported, and are shown to depend on the laser field intensity and orientation relative to the laser polarization. The ionization dynamics, including ionization from individual Kohn–Sham orbitals, is also reported. The orbital ionization dynamics are shown to vary with an increase in the intensity of the laser field.
##### Ab-initio study of the structural and electronic properties of cadmium telluride CdTe and its alloy cadmium zinc telluride Cd1-xZnxTe

Louiza Messaadia , Physics Laboratory at Guelma, University 8 Mai 1945 Guelma, Department of Materials Science, Algeria

Abstract: The scope of this investigation is to make a clear contrast between the structural and electronic properties of cadmium telluride CdTe and its alloy cadmium zinc telluride Cd1-xZnxTe using first principles calculations based on density functional theory within the local density approximation (LDA). A supercell with 2×2×2 unit cells has been considered in the calculations. Structural parameters and electronic densities of states are shown to discuss the effects of local environment induced by Zn impurity on the structure of zinc- blend CdTe. Our results are in good agreement with some theoretical studies.
##### Ab-initio study of the structural and electronic properties of cadmium telluride CdTe and its alloy cadmium zinc telluride Cd1-xZnxTe

Louiza Messaadia , Physics Laboratory at Guelma, University 8 Mai 1945 Guelma, Department of Materials Science, Algeria

Abstract: The scope of this investigation is to make a clear contrast between the structural and electronic properties of cadmium telluride CdTe and its alloy cadmium zinc telluride Cd1-xZnxTe using first principles calculations based on density functional theory within the local density approximation (LDA). A supercell with 2×2×2 unit cells has been considered in the calculations. Structural parameters and electronic densities of states are shown to discuss the effects of local environment induced by Zn impurity on the structure of zinc- blend CdTe. Our results are in good agreement with some theoretical studies.
##### Canonical Quantization of the Reissner-Nordstr\"om Black Hole via Conditional Symmetries

Evangelos Melas, Technological Educational Institution of Patras, Department of Management, Greece

Abstract: We use the conditional symmetry approach to study the $r$-evolution of a minisuperspace spherically symmetric model both at the classical and quantum level. After integration of the coordinates $t$, $\theta$ and $\phi$ in the gravitational plus electromagnetic action the configuration space dependent dynamical variables turn out to correspond to the $r$-dependent metric functions and the electrostatic field. In the context of the formalism for constrained systems (Dirac - Bergmann, ADM) with respect to the radial coordinate $r$, we set up a point-like reparametrization invariant Lagrangian. It is seen that, in the constant potential parametrization of the lapse, the corresponding minisuperspace is a Lorentzian three-dimensional flat manifold which obviously admits six Killing vector fields plus a homothetic one. The weakly vanishing $r$-Hamiltonian guarantees that the phase space quantities associated to the six Killing fields are linear holonomic integrals of motion. The homothetic field provides one more rheonomic integral of motion. These seven integrals are shown to comprise the entire classical solution space, i.e. the space-time of a Reissner-Nordstr\"om black hole, the $r$-reparametrization invariance since one dependent variable remains unfixed, and the two quadratic relations satisfied by the integration constants. We then quantize the model using the quantum analogues of the classical conditional symmetries, and show that the existence of such symmetries yields solutions to the Wheeler-DeWitt equation which, as a semiclassical analysis shows, exhibit a good correlation with the classical regime. Finally, we use the resulting wave functions to investigate the possibility of removing the classical singularities.
##### A mathematical characterization of the gel point in sol-gel transition

Ayse Humeyra Bilge, Kadir Has University, Fculty of Engineerin and Natural Sciences, Turkey

Abstract: We model the sol-gel transition in terms of Susceptible-Infected-Removed (SIR) and Susceptible-Exposed-Infected-Removed (SEIR) models and compare with experimental results. We show, numerically, that the gel point" described as the onset of the gelation phenomena and measured experimentally, corresponds to an accumulation point of the extreme values of the derivatives of the gelation curve. We define the critical point of a sigmoidal curve" as the limit of the points where the derivatives reach their extreme values, provided that this limit exists.
##### Multiscale continuation algorithms for binary Rydberg-dressed Bose-Einstein condensates

Cheng-Sheng Chien, Chien Hsin University of Science and Technology, Department of Computer Science and Information Engineering, Taiwan

Abstract: We present two multiscale continuation algorithms for binary Rydberg-dressed Bose-Einstein condensates which are governed by a system of the Gross-Pitaevskii equations (GPEs). First we describe a three-parameter continuation algorithm to trace the ground state solution curve of the GPEs, where we use the chemical potentials λ1 and λ2 together with an artificial parameter ν0 belong to [0, ν*] as the three continuation parameters. Here ν0 is used to control the value of the parameter μ12 for the integration term, say μ12 = ν0 × μ0 for some μ0. Besides, the parameters μ11 = μ22 are fixed. Next, we describe a four-parameter continuation algorithm to trace the ground state solutions of the GPEs. At the beginning we use the two-parameter continuation algorithm described by Wang and Chien [CiCP, 2013] to trace the ground state solution curves of the GPEs with the coefficients of the integration terms μ11 = μ22 = μ12 = 0. When the constraint conditions ||ψ1|| = ||ψ2|| = 1 are satisfied, we use the chemical potentials λ1 and λ2 together with two additional parameters ν0 belong to [0, ν*] and ν1 belong to [0, ν**] as the four continuation parameters and proceed again to trace the ground state solution curves, where μ11 = μ22 = ν1 × μ1 for some μ1. Since the parameters μ11 = μ22 and μ12 have different scales, the parameter ν1 reaches the desired value earlier than the parameter ν0. Then we reduce the four-parameter continuation algorithm to the three-parameter one, and proceed to trace the ground state solution curves until the desired ground state solutions are obtained. The proposed multiscale/multi-parameter continuation algorithms have the advantage that we can obtain the contours of the wave functions ψ1 and ψ2 for various values of μ11 = μ22 belong to [0, μ*] and μ12 belong to [0, μ**], where we only need to trace the ground state solution curve once. Our numerical experiments show that the proposed algorithms outperform the classical continuation algorithm, and is very competitive compared to other numerical methods for treating similar problems.
##### BICEP2, inflation and quantum fluctuations

Prof. Dr. Emilio Elizalde, ICE-CSIC and IEEC Barcelona, , Spain

Sergei Odintsov, ICREA and ICE (CSIC-IEEC), Theoretical Physics and Cosmology , Spain

Abstract: After a short review of a few basic notions around the Big Bang theory, a summary will be given of the impacting results of BICEP2, together with the subsequent reactions of cosmologists and theoretical physicists. While some basic models of inflation can accommodate the new results, they seem to exclude, on the contrary, other popular and very important families of inflationary models. And it is quite difficult to explain, simultaneously, the previous results of PLANCK together with the new data. The fundamental question has also arisen: have we finally found indeed the footprints of the very elusive quantum gravity?
##### Unification of inflation with dark energy era in modified gravity

Sergey Odintsov, ICREA and ICE(CSIC-IEEC), Cosmology, Spain

Abstract: We review the unification of inflation with dark energy within modified F(R) gravity. Several models are presented, their properties are described. Some applications- like massive neutron stars from F(R) gravity are also mentioned. Finally, the unification of early-time with late-time acceleration is also considered in modified Gauss-Bonnet gravity.
##### Supersymmetric extension of a coupled Korteweg-de Vries system

Adrián Sotomayor, Antofagasta University, Mathematics, Chile

Alvaro Restuccia, Antofagasta University, Physics, Chile

Abstract: We obtain a supersymmetric extension of a coupled KdV system. This system arises from a complexiﬁcation of real KdV equation. The system has an inﬁnite sequence of conserved quantities involving the even and odd ﬁelds describing the new system. This inﬁnite se- quence can be obtained from a supersymmetric Gardner equation via a Gardner integrable deformation. We also ﬁnd the hamiltonian and the Poisson structure of the system.
##### One dark matter mystery: halos in the cosmic web

Abstract: The current cold dark matter cosmological model explains the large scale cosmic web structure but is challenged by the observation of a relatively smooth distribution of matter in galactic clusters. We consider various aspects of modeling the dark matter around galaxies as distributed in smooth halos and, especially, the structure of the dark matter halos seen in cosmological N-body simulations. We conclude that the problems of the cold dark matter cosmology on small scales are more serious than normally admitted.
##### Stabilizability of stochastic nonlinear hybrid systems

Ewelina Seroka, Cardinal Stefan Wyszynski University in Warsaw, Faculty of Mathematics and Natural Sciences, Poland

Abstract: The problem of the asymptotic stabilizability in probability of a class of stochastic nonlinear control hybrid systems (with a linear dependence of the control) with any, state dependent and markovian switching rule is considered in the paper. It is assumed that the trivial solution of unforced hybrid system is stable in probability (wherein some of subsystems of unforced hybrid systems can be unstable). By applying the stabilizing control the trivial solution of hybrid system becomes asymptotically stable in probability. To solve the issue, the Lyapunov technique including a common, single, multiple Lyapunov function, the hybrid control theory and some results of P. Florchinger for stochastic non-hybrid systems, are used. Moreover, in the case of hybrid systems with markovian switching rule, results of R.Z. Khasminskii, C. Zhu, G. Yin are applied. Sufficient conditions for the asymptotic stabilizability in probability for a considered class of hybrid systems are formulated. Also the stabilizing control in a feedback form is considered. Furthermore, in the case of hybrid systems with the state dependent switching rule, a method for a construction of stabilizing switching rules is proposed. Obtained results are illustrated by examples and numerical simulations.
##### Microscopic black hole stabilization via the uncertainty principle

Constantinos G. Vayenas, University of Patras, , Greece

Abstract: We study the properties of bound states consisting of relativistic rotating particles [1] and formed via stabilization of microscopic black holes due to the uncertainty principle [1-3]. We show that the combination of Schwarzschild or Reissner-Nordström geodesics with the uncertainty principle, expressed via the Compton wavelength of the confined mass, leads to Planckian or sub-Planckian states with the mass and other properties of hadrons, formed via the confinement of much lighter relativistic particles. The concepts of generalized Uncertainty Principle (GUP) [4] and Generalized Event Horizon (GEH) [3] in relation to these states are also discussed. 1. C.G. Vayenas, S. Souentie, A. Fokas, Physica A, 405 (2014) 360-379. 2. C.G. Vayenas, S. Aretakis, A. Fokas and D. Grigoriou, in preparation (2014). 3. B.J. Carr, arXiv:1402.1427v1 [gr-qc] (2014). 4. S. Das, E.C. Vagenas, Can. J. Phys. 87 (2009) 233-240.
##### Computer Simulation of Fire Dynamics in Industrial Hall

Jan Glasa, Institute of Informatics, Slovak Academy of Sciences, , Slovakia (Slovak Republic)

Abstract: In this paper, computer simulation of smoke spread dynamics in industrial hall is investigated. A set of simulations of fire in three industrial halls with the same geometry varying in the height of ceiling is realized using the CFD-based fire field model, FDS, version 6. The obtained simulation results are described focusing on the impact of the ceiling height and fire barriers on the fire course and smoke spread dynamics.
##### Feedback Stabilization of an Oscillating Vertical Cylinder by POD Reduced-Order Model

Gilles Tissot, PPRIME Institute, Fluides, Thermique, Combustion, France

Laurent CORDIER, PPRIME Institute, Fluides, Thermique, Combustion, France

Bernd NOACK, PPRIME Institute, Fluides, Thermique, Combustion, France

Abstract: The objective of this paper is to demonstrate the use of Reduced-Order Models (ROM) based on Proper Orthogonal Decomposition (POD) to stabilize by vertical oscillations the flow over a circular cylinder for a Reynolds number equal to 60. The 2D Navier-Stokes equations are first solved by a finite element method with COMSOL Multiphysics in which the moving mesh is introduced via ALE. Since in Fluid-Structure Interaction, the POD algorithm cannot be applied directly, we then implement the fictitious domain method of Glowinski et al. (1999) where the solid domain is treated as a fluid undergoing an additional constraint. The POD-ROM is then classically obtained by projecting the Navier-Stokes equations on a small number of POD modes. At this level, the cylinder movement is enforced in the POD-ROM through the introduction of Lagrange multipliers. A Linear Quadratic Regulator framework is used to determine the optimal control law, in our case the vertical velocity of the cylinder, such that the flow is stabilized. After linearization of the POD-ROM around the steady flow state, the optimal linear feedback gain is obtained as solution of a Generalized Algebraic Riccati Equation. Finally, when the optimal feedback control is applied, it is shown that the flow converge rapidly to the steady state. In addition, a vanishing control is obtained proving the efficiency of the control approach.
##### Inference of Transcriptional Network for Pluripotency in mouse Embryonic Stem Cells

Sachiyo Aburatani, National Institute of Advanced Science and Technology, Computational Biology Research Center, Japan

Abstract: In embryonic stem cells, some transcription factors (TFs) are known to maintain the pluripotent process. To gain insights into the regulatory system to control pluripotency, I inferred regulatory relationships between TFs, which expressed in ES cells. In this study, I applied a method based on structural equation modelling (SEM), combined with factor analysis, to 649 expression profiles of 19 TF genes measured in mouse ES (mES) cells. By the factor analysis, 19 TF genes were regulated by several unmeasured factors. Since the known cell reprogramming TF genes (Pou5f1, Sox2, Nanog and Klf4) were regulated by the different factors, the each estimated factor is considered to be an input signal transduction to control pluripotency in mES cells. In the inferred network model, TF proteins were also arranged as unmeasured factors which control the other TFs. The interpretation of the inferred network model allowed us to reveal the regulatory mechanism for controlling pluripotency in ES cells.
##### MPR Supported Multiscale Initialization of Transfer Functions for Interactive 3D Visualization

Merve Ozdemir, Dokuz Eylül University, Electrical and Electronics Engineeing, Turkey

Alper Selver, Dokuz Eylül University, Electrical and Electronics Engineering, Turkey

Oguz Dicle, Dokuz Eylül University, Radiology, Turkey

Abstract: Visualization aims to produce clear and informative pictures of the important structures in a dataset. Depending on the application, this requires interactive determination of visual parameters such as opacity and color. In volume rendering technique, combinations of these visual parameters can be determined during the rendering pipeline. During the generation of volume rendered images, Transfer Function (TF) specification is the step where these adjustments can be done. Therefore, it is crucial and important to design accurate TFs to produce meaningful and intelligible 3-D images. However, TF design is a very difficult task because of the availability of various possibilities in extensive search spaces of TFs. Since this flexibility of search space cannot be kept in strict bounds, specification of an appropriate TF is a challenging problem where effective initial TF designs should be generated prior to the optimization that is controlled by the user. Moreover, advanced user interaction interfaces and data exploration tools should be provided for fulfilling user expectations. To overcome the difficulty of initial TF generation generally a number of predefined TF presets are used as starting point (so called initial TF design). The main idea behind this approach is that certain types of volume data are standardized in the range of data values and special sub-ranges are assigned to the same type of structure (Thus, predefined TFs are adjusted due to these ranges). However, volumetric data usually have varying characteristics even in different samples of the same application. For instance, in medical imaging, depending on different modality settings, injection of a contrast media or environmental circumstances, the sub-ranges of the tissues may vary significantly. Similarly in natural stone analysis, the ingredients of brecciated rocks differ significantly from each other which prevent the use of pre-defined ranges for analysis. For these reasons, a limited number of TF presets cannot be enough to cover all possible cases and to provide useful initial TFs. In order to create a useful initial TF that provides a good basis prior to optimization, an automatic sub-range detection method that finds the intensity range for each structure of interest is needed. Moreover, it is necessary to integrate the developed method into the TF design procedure without losing user control and interaction over the search space. By addressing this problem, a semi-automatic method for initial generation of TFs is introduced in this study. The proposed approach is based on modeling the search space in a hierarchical manner using gaussian functions. The search space is constructed with generalized Volume Histogram Stack (VHS) which can be determined by the user by interacting with Multi Planar Reconstruction (MPR) images of the volume data. VHS is recently introduced in as a new domain which is created by aligning the histograms of the image slices of a CT/MR series. Histograms were generated from orthogonal directions of slice planes, namely, axial, coronal and saggital. Thus, VHS can represent the intensity values of the tissues as well as their spatial information and local distributions (via lobes in VHS) which are not available in conventional volume histograms. The tissues which are at different slices but with similar gray level distributions can clearly be distinguished by using this spatial information. Then, a tissue (a structure of interest) can effectively be visualized by determining its corrsponding lobe(s) in VHS, which represents that structure of interest, and by assigning a color-opacity value to that lobe. In this study, VHS data is further generalized in such a way that it is possible to calculate VHS on an arbitrary aligned axis (i.e. a slice plane with an arbitrary normal vector, not aligned with x, y or z) which can be obtained with the help of not only MPR but also of Curved MPR or Oblique sectioning techniques. A generalized VHS can be generated not only for slice based medical image series but all volumetric data as long as spacing values in 3D is known. With the help of this expansion, the VHS becomes an effective new domain as a search space for TF specification on any kind of 3D data. Although extending the implementation of VHS to work along an arbitrary axis or with an arbitrary surface would provide many advantages, it should still be supported by a lobe detection method. Because, when two or more structure of interest have overlapping intensity values and if their spatial location also overlaps at some degree, they construct overlapping regions which produce hardly recognizable minor lobes (or side parts of main lobes) due to the domination by the major ones. The proposed multiscale and hierarchical modeling strategy applied on generalized VHS allows recognizing suppressed lobes corresponding to suppressed structures and also for representing the overlapping regions which are parts of the lobes but can not be represented by the clusters associated to the lobes due to the overlapping. The developed strategy allows the integration of spatial knowledge, local distribution of the structures and their intensity information into the TF while preserving the user control. The proposed method is applied to 10 MRI datasets for abdominal tissue/organ visualization. The results show that, the proposed approach effectively increases performance on classification of abdominal organs.
##### Some of the cosmological funds gravitational theories F(T,\Theta)

Faezeh, Kiani, , Iran (Islamic Republic of)

Abstract: we study a new type of modified teleparallel gravity of the form F(T,\Theta) in which T, the torsion scalar, is coupled with $\Theta$, the trace of the stress-energy tensor. We also study the minimal and non-minimal coupling by the dynamical system approach. In this Model, the crossing the phantom divide, the existence of de Sitter solution and its stability are investigated in the non-minimal version.
##### Automatic Method to Classify Images Based on Multiscale Fractal Descriptors and Paraconsistent Logic

Eduardo Pavarino, São Paulo State University (UNESP), Department of Computer Science and Statistics (DCCE), Brazil

Leandro Alves Neves, São Paulo State University (UNESP), Department of Computer Science and Statistics - DCCE, Brazil

Marcelo Zanchetta do Nascimento, Federal University of Uberlândia, FACOM, Brazil

Moacir Fernandes de Godoy, FAMERP - São José do Rio Preto, Transdisciplinary Center for Study of Chaos and Complexity (NUTECC), Brazil

Pedro Francisco de Arruda, FAMERP - São José do Rio Preto, , Brazil

Leandro Neves, São Paulo State University, DCCE, Brazil

Dalísio de Santi Neto, Hospital de Base de São José do Rio Preto, Department of Pathology, Brazil

Abstract: In this study is presented an automatic method to classify images from fractal descriptors as decision rules, such as multiscale fractal dimension and lacunarity. The proposed methodology was divided in three steps: quantification of the regions of interest with fractal dimension and lacunarity, techniques under a multiscale approach; definition of reference patterns, which are the limits of each studied group; and, classification of each group, from the combination of the reference patterns with signals maximization (an approach commonly considered in paraconsistent logic). As a first application, the proposed method was used to classify histological prostatic images with 40x of magnification, aiming the diagnostic of prostate cancer. The investigated groups were normal, hyperplasia and cancer. The accuracy levels were important, overcoming those obtained with Support Vector Machine (SVM), Radial Basis Function Network (RBFNetwork) and Best-first Decicion Tree (BFTRee) classifiers: well know and widely applied techniques to recognize patterns. The proposed technique was better than SVM and RBFNetwork in all tests with both stroma and lumen and only with stroma, respectively. Regarding the BFTree classifier, the developed method was better in, approximately, 67% of the comparisons. The groups with better results were normal versus hyperplasia and normal versus cancer, being the region of stroma the more significant to distinguish the prostate cancer. Thus, the proposed approach allows recognize and classify patterns in studied context, offering the advantage of giving comprehensive models to the specialists.
##### Approximate best proximity points of cyclic self-mapsand cyclic asymptotic regularity

Manuel de la Sen, University of the Basque Country, IIDP, Spain

Abstract: This manuscript relies on 2-cyclic self-mappings by giving a simple formal development to join some useful concepts, and related results, like those of approximate best proximity points of cyclic self-mappings, the approximate best proximity (respectively, partial best proximity) point property or the cyclic asymptotic regularity. A formal mathematical development is given which brings together the concepts of approximate best proximity points of cyclic self-mappings, approximate best proximity (respectively, partial best proximity) point property and cyclic asymptotic regularity of cyclic self-mappings. There are also further induced results available for the approximate fixed points of the composite self-mapping obtained from the composition of the original one with itself which follow from the approximate best proximity point properties and those related ones of cyclic asymptotic regularity .
##### Clusters and molecules in extreme light

Eric Suraud, Universite Paul Sabatier, Lab. Phys. Theorique, France

Abstract: The progress in laser technology over the last decades has opened up new avenues for the explo- ration of properties of clusters and molecules. A laser pulse is characterized by its frequency but also by the laser intensity as well as the laser time profile. While for years the variations of these parame- ters were heavily constrained by technology, the last two decades and even more so the last years have seen tremendous increases in the range of attainable parameters. This is true for intensity, which since the 1990’s can reach huge values which can lead to very large energy deposits and possibly violent disintegration of the irradiated species. But this is also true for the tuning of the time profile which can now be tailored up to time scales of the order of magnitude of electronic motion and even below. This allows the follow up of the detail of electronic dynamics at its own ”natural” time. The latest breaktroughs were attained in terms of laser frequency with the ongoing possibility of reaching very large frequencies up the X domain. This opens up new possibilities of imaging which are progressively being explored. We shall discuss some of these directions of investigation, taking examples in cluster and molecular physics. We shall especially discuss the case of high intensity and short time pulses for which a sizable amount of results have already been attained. We shall also discuss in detail the case of very short times (attoseconds) which are becoming more and more studied. [1] Mechanisms of cluster ionization in strong laser pulses, U. Saalmann, C. Siedschlag,, J. M. Rost, J. Phys. B 39 (2006) R39. [2] Non linear electron dynamics in metal clusters, F. Calvayrac, P. G. Reinhard, E. Suraud, C. Ullrich, Phys. Reports 337(2000)493-578 [3] Laser-driven nonlinear cluster dynamics, Th. Fennel, K.-H. Meiwes-Broer, J. Tiggesbumker, P.-G. Reinhard, P. M. Dinh, E. Suraud, Rev. Mod. Phys. in press, 2010, http://arxiv.org/abs/0904.2706 [4] Probing Time-Dependent Molecular Dipoles on the Attosecond Time Scale, C. Neidel et al Phys. Rev. Lett 111 (2013) 033001.
##### Multi-parameter continuation and collocation methods for rotating multi-component Bose-Einstein condensates

Yun-Shih Wang, Chien Hsin University of Science and Technology, , Taiwan

Abstract: We describe multi-parameter continuation methods combined with spectral collocation methods (SCM) for computing numerical solutions of rotating two-component Bose-Einstein condensates (BECs), which are governed by the Gross-Pitaevskii equations (GPEs). Various types of orthogonal polynomials are used as the basis functions for the trial function space. A novel multi-parameter/multiscale continuation algorithm is proposed for computing the solutions of the governing GPEs, where the chemical potential of each component and angular velocity are treated as the continuation parameters simultaneously. The proposed algorithm can effectively compute numerical solutions with abundant physical phenomena. Numerical results on rotating two-component BECs are reported.
##### Two reliable approaches involving Haar wavelet method and Optimal Homotopy Asymptotic method for the solution of fractional Fisher type equation

SANTANU SAHA ROY, NATIONAL INSTITUTE OF TECHNOLOGY ROURKELA, MATHEMATICS, India

Abstract: In this article, two reliable techniques, Haar wavelet method and optimal homotopy asymptotic method (OHAM) are presented. Haar wavelet method is an efficient numerical method for the numerical solution of fractional order partial differential equation like Fisher type. The approximate solutions of the fractional Fisher type equation are compared with the optimal homotopy asymptotic method as well as with the exact solutions. Comparisons between the obtained solutions with the exact solutions exhibit that both the featured methods are effective and efficient in solving nonlinear problems. However, the results indicate that optimal homotopy asymptotic method provides more accurate value than Haar wavelet method.
##### Atomistic and Continuum Simulations of Nanoscale Pattern formation on Alloy Surfaces by Ion Beam Sputtering

Bharathi Srinivasan, Institute of High Performance Computing, Engineering Mechanics, Singapore

Ramanarayan Hariharaputran, Insititute of High Performance Computing, Materials Science & Engineering, Singapore

Yong Wei Zhang, Institute of High Performance Computing, EM, Singapore

Abstract: Nanoscale surface pattern creation is a topic of recent research interest due to its relevance in optoelectronic devices. Sputtering by low energy ion beam is one of the techniques used to achieve self-assembled morphological patterns on surfaces. In alloy surfaces, both morphological and compositional patterns are created due to interplay of various factors such as diffusion, sputter removal and elastic misfit. Temperature, flux of the incident ion beam and the elastic interactions arising out of misfit stresses dominate the self-assembly of these nanostructures. Here we present the results from both kinetic Monte Carlo simulations and a nonlinear continuum model on the influence of these parameters on the dynamics of pattern formation in alloy surfaces by ion beam sputtering.
##### Some properties of light mesons in a strong magnetic field in context of the SU(3) lattice gauge theory

Elena Luschevskaya, Institute of Theoretical and Experimental Physics, Laboratory of theoretical and nuclear physics, Russian Federation

Abstract: We calculated the correlators of vector, axial and pseudoscalar currents in external strong abelian magnetic field in $SU(3)$ gluodynamics. The masses of neutral and charged $\rho$ and $A$ mesons with various spin projections to the external magnetic field $B$ have been calculated. We found that the masses of neutral mesons with zero spin $s=0$ decrease, while the masses of the $\rho$ and $A$ mesons with spin $s=\pm 1$ increase in increasing magnetic field. The masses of charged vector and axial mesons with nonzero spins diminish, but the masses of the same mesons with zero spin grow with the field value.
##### Trans-dimensional Monte Carlo sampling applied to the magnetotelluric inverse problem

Eric Mandolesi, Dublin Institute for Advanced Studies, Geophysics, Ireland

Nicola Piana Agostinetti, Dublin Institute for Advanced Studies, Geophysics, Ireland

Abstract: The data required to build geological models of the subsurface are often unavailable from direct measurements or well logs. In order to image the subsurface geological structures several geophysical methods have been developed. The magnetotelluric (MT) method uses natural, time-varying electromagnetic (EM) fields as its source to measure the EM impedance of the subsurface. The interpretation of these data is routinely undertaken by solving inverse problems to produce 1D, 2D or 3D electrical conductivity models of the subsurface. In classical MT inverse problems the investigated models are parametrized using a fixed number of unknowns (i.e. fixed number of layers in a 1D model, or a fixed number of cells in a 2D model), and the non-uniqueness of the solution is handled by a regularization term added to the objective function. This study presents a different approach to the 1D MT inverse problem, by using a trans-dimensional Monte Carlo sampling algorithm, where trans-dimensionality implies that the number of unknown parameters is a parameter itself. This construction has been shown to have a built-in “Occam’s razor”, so that the regularization term is not required to produce a simple model. The influences of subjective choices in the interpretation process can therefore be sensibly reduced. The inverse problem is solved within a Bayesian framework, where posterior probability distribution of the investigated parameters are sought, rather than a single best-fit model, and uncertainties on the model parameters, and their correlation, can be easily measured.
##### Classification of Histological Images Based on the Stationary Wavelet Transform

Sidon, Federal University of ABC , CMCC Centro de Matemática, Computação e Cognição, Brazil

Yan Duarte, Federal University of ABC , CMCC Centro de Matemática, Computação e Cognição, Brazil

Marcelo Zanchetta do Nascimento, Federal University of Uberlândia, FACOM, Brazil

Leandro Neves, São Paulo State University, DCCE, Brazil

Valério Batista, Federal University of ABC, CMCC, Brazil

Abstract: Non-Hodgkin lymphomas are of many distinct types, and different classification systems make it difficult to diagnose them correctly. Many of these systems classify lymphomas only based on how they look under a microscope. In 2008 the World Health Organisation (WHO) introduced the most recent system, which also considers the chromosome features of the lymphoma cells and the presence of certain proteins on their surface. The WHO system is the one that we apply in this work. Herewith we present an automatic method to classify histological images of three types of non-Hodgkin lymphoma. Our method is based on the Stationary Wavelet Transform (SWT), and it consists of three steps: 1) extracting sub-bands from the histological image through SWT, 2) applying Analysis of Variance (ANOVA) to clean noise and select the most relevant information, 3) classifying it by the Support Vector Machine (SVM) algorithm. The kernel types Linear, RBF and Polynomial were evaluated with our method applied to 339 images of lymphoma from the National Institute on Aging. We concluded that the following combination led to the most relevant results: detail sub-band, ANOVA and SVM with Linear and RBF kernels.
##### Family of higher order exponential variational integrators

Odysseas Kosmas, University of Erlangen-Nuremberg, Mechanical Engineering, Germany

Sigrid Leyendecker, University of Erlangen-Nuremberg, Mechanical Engineering, Germany

Abstract: In the present work we derive a family of higher order exponential variational integrators for the numerical integration of systems with oscillatory solutions. To increase the order of variational integrators, first the discrete Lagrangian in any time interval is defined as a weighted sum of the evaluation of the continuous Lagrangian at intermediate time nodes while expressions for configurations and velocities are obtained using interpolating functions that can depend on free parameters. Secondly, in order to chose those parameters appropriately, exponential integration techniques are embedded. Focusing on highly oscillatory problems, we then split their potential energy into a fast and a slow component and use different quadrature rules for the different potentials. Finally, we study the behavior of this family of integrators in numerical tests.
##### Analysis on the use of Multi-Sequence MRI Series for Segmentation of Abdominal Organs

Esref Selvi, Dokuz Eylül University, Institute of Natural and Applied Sciences, Turkey

Emre Kavur, Dokuz Eylül University, Institute of Natural and Applied Sciences, Turkey

Alper Selver, Dokuz Eylül University, Electrical and Electronics Engineering, Turkey

Oguz Dicle, Dokuz Eylül University, Radiology, Turkey

Abstract: Segmentation of abdominal organs from MRI data sets is a challenging task due to various limitations and artefacts. During the routine clinical practice, radiologists use multiple MR sequences in order to analyze different anatomical properties. These sequences have different characteristics in terms of acquisition parameters (such as contrast mechanisms and pulse sequence designs) and image properties (such as pixel spacing, slice thicknesses and dynamic range). For a complete understanding of the data, computational techniques should combine the information coming from these various MRI sequences. These sequences are not acquired in parallel but in a sequential manner (one after another). Therefore, patient movements and respiratory motions change the position and shape of the abdominal organs. In this study, the amount of these effects is measured using three different symmetric surface distance metrics performed to three dimensional data acquired from various MRI sequences. The results are compared to intra and inter observer differences and discussions on using multiple MRI sequences for segmentation and the necessities for registration are presented.
##### Spherical Target Recognition using Time-Domain Multiscale Approximation of Scattered Signals

Mehmet Taygur, Dokuz Eylül University, Electrical and Electronics Engineeing, Turkey

Alper Selver, Dokuz Eylül University, Electrical and Electronics Engineering, Turkey

Yesim Zoral, Dokuz Eylül University, Electrical and Electronics Engineeing, Turkey

Abstract: Classification of similar shaped objects from scattered electromagnetic waves is a difficult problem to solve, as it heavily depends on the aspect angle. Eliminating the effects of the aspect angle is possible by extracting distinguishable features from the scattered signals. These features should be robust to noise effects especially at SNR levels, where noise effects become dominant on the scattered signal. In this paper, we propose a target classification method, which uses a structural feature set extracted from scattered signal. Prior to feature extraction, a multi-scale approximation is performed using hierarchical radial basis function network topology to suppress the effects of noise on scattered signal. After principle component analysis, k-fold cross validation based experiments is performed. Results show that spherical targets are recognized successfully up to -10dB SNR.
##### A Mathematical Model for Analyzing Surface Plasmon Polariton Propagation through a Schottky-Junction-Based Plasmonic Amplifier

Abdolber Mallah Livani, Amirkabir University, Electrical Engineering, Iran (Islamic Republic of)

Hassan Kaatuzian, Amirkabir University of Technology, Photonics Research Laboratory (PRL), Electrical Engineering Dept., Iran, Islamic Republic Of

Abstract: We have devised a simple model for describing an amplifier which operates on Surface Plasmon Polaritons (SPPs). A semiconductor is considered instead of dielectric which its interface with metal can support TM polarized SPP propagation. Metal-Semiconductor interface in particular conditions can be regarded as a Schottky junction that has the capability of being pumped electrically. So compensation of propagation losses is possible and beyond that, amplification occurs. This configuration has advantages such as simple fabrication process and compact size. However, designing such an amplifier with too many effects that arise in a Schottky junction may be an extremely difficult process. So a simplified model which regards essential effects and ignores non important ones is useful. In this work, gold has brought in contact with n+-doped In0.53Ga0.47As as semiconductor to form a Schottky junction. Concentration of doping in semiconductor is 9.3×1017cm-3. To estimate optical gain, at first, we solve drift-diffusion currents, continuity and Poisson’s equations to calculate carrier densities then using wave analysis, power density will be obtained. In 0.8V forward bias, optical gain of In0.53Ga0.47As is estimated around 2000cm-1 which is required for compensation of ohmic losses in metal for SPP. Current density in this forward bias is 104A/cm2 which is a typical current density.
##### Multicriteria optimization for wind farm layout design and assessment: Effective implementation of the Weighted Sum and Lexicographical methods

Daniela Borissova, Information and Communication Technologies – Bulgarian Academy of Sciences, Information Processes and Decision Support Systems, Bulgaria

Ivan Mustakerov, Information and Communication Technologies – Bulgarian Academy of Sciences, Information Processes and Decision Support Systems, Bulgaria

Abstract: The design of a reliable and cost-effective industrial wind power farm is a prerequisite for the effective use of wind power as an alternative resource of renewable energy. The optimal layout could be defined by proper mathematical modeling involving explicitly multiple criteria in decision-making environment. The multicriteria approach is well suited for testing various design scenarios. Having in mind that designers usually prefer more flexible and transparent approaches, the multicriteria optimization could be used as a wind farm design simulation tool. In contrast to single criterion optimization, multiple criteria decision making is concerned with solving decision problems to choose the “best” alternative from a set of available alternatives, where the “best” can be interpreted as “the most preferred alternative” of decision maker. The multicriteria methods reflect more precisely the decision maker preferences in respect of wind farm layout design. The paper presents a multicriteria optimization approach for wind farm layout design by formulation of mixed integer problem for simultaneously determination of number, type and turbines placement taking into account the specifics of given wind area and decision maker preferences. A methodology for assessment of wind farm layout design by multicriteria optimization is proposed. The solutions results assist the project manager to assess the wind farm project on the early design stage. The weighted sum and lexicographical method are used to assess the effectiveness of wind farm layout design. Numerical case study examples of wind farm layout design demonstrate that both methods adequately reflect different criteria importance for different wind conditions of a particular site. During the simulation process, different points of view are analyzed to get Pareto-optimal alternatives. The wind farm energy output and costs corresponding to different Pareto-optimal alternatives under different wind directions are illustrated.
##### An algorithm for optimal design of wind farm layout over complex terrain

Ivan Mustakerov, Information and Communication Technologies – Bulgarian Academy of Sciences, Information Processes and Decision Support Systems, Bulgaria

Daniela Borissova, Information and Communication Technologies – Bulgarian Academy of Sciences, Information Processes and Decision Support Systems, Bulgaria

Abstract: The paper concerns problem for optimal design of wind farm layout over complex terrain. When trying to make wind energy project economically viable, wind turbines type and placement are important design parameters that constitute a significant part of the overall project costs. Generally, the energy production is the major factor that has to be considered for economic success of a wind farm project, but in reality it has to be evaluated toward associated costs. It is essential to evaluate the relation of costs and energy production of wind farm on the early design stages. A combinatorial optimization model for determination of the wind turbines type and wind farm layout is described. This model is used in an algorithm for evaluating of different layouts over given terrain with preferable and not preferable zones for wind turbines placement. The basic steps of the algorithm are: definition of first k optimal types and corresponding numbers of wind turbines; determination of wind farm layouts for each type of turbines for given wind site specifics; calculation of the expected farm output potential and related costs for each layout; choice of the turbine’s type toward the best ratio of wind farm output potential and related costs. The algorithm is tested numerically for real turbines type’s data for given site dimensions and terrain specifics. The numerical illustration showed the applicability of the proposed algorithm. Using of this approach contributes the preliminary evaluation of wind farm effectiveness on the design stage of wind farm project development. It can be implemented as reasonable decision making tool for effective wind farm design. Some illustrative layouts for different wind turbines types defined by means of the described algorithm for uniform and predominant wind directions are shown.
##### Cooperative particle motion and dissipative solitary waves in complex (dusty) plasmas

Sergey Zhdanov, MPE, , Germany

Abstract: As an important element of self-organization, cooperative particle motion is present in many physical, astrophysical and biological systems. A particularly interesting and challenging topic is to study dynamic cooperativity at local and intermediate scales. As a rule, cooperative dynamics, bringing to life ’abnormal’ effects like enhanced diffusion, self-dragging, or self-propelling of particles, hold aspects of ’strange’ kinetics. Such kind of cooperative behavior was evidenced for string-like formations of colloidal rods, dynamics of mono- and di-vacancies in 2d colloidal crystals and in complex plasmas. Strongly coupled complex (dusty) plasmas give us a unique opportunity to go beyond the limits of continuous media and study various generic processes occurring in liquids or solids in real time and at the kinetic level. There is a certain advantage to experiment with complex plasmas merely because these systems are easy to manipulate in a controllable way. Externally manipulated ’dust molecules’, self-assembled strings in driven 3d particle clusters and dissipative solitary waves are noticeable examples. The results of recent experiments showing microparticle cooperative movements occurring under natural conditions are reviewed and interpreted.
##### GPU Computing in Bayesian Inference of Realized Stochastic Volatility Model

Tetsuya Takaishi, Hiroshima University of Economics, , Japan

Abstract: Recently the realized stochastic volatility model has been proposed to infer volatility of financial time series. We perform the Bayesian inference of the realized stochastic volatility model by the Hybrid Monte Carlo algorithm. The Hybrid Monte Carlo algorithm can be parallelized and thus performed on the GPU. The GPU code is developed with CUDA Fortran. We compare the computational time in performing the Hybrid Monte Carlo algorithm on GPU (GTX 760) and CPU (Intel i7-4770 3.4GHz) and find that the GPU can be up to 17 times faster than the CPU. We also code the program with OpenACC and find the similar speedup with OpenACC.
##### Magnetic phase transition in the spin-fluctuation theory

Nikolai Melnikov, Lomonosov Moscow State University, Computational Mathematics and Cybernetics, Russian Federation

Georgiy Paradezhenko, Lomonosov Moscow State University, Computational Mathematics and Cybernetics, Russian Federation

Abstract: The effect of spin fluctuations on the magnetic phase transition is studied by the functional integral method. The pair interaction between magnetic moments at finite temperatures is replaced by the interaction with the fluctuating exchange field. Magnetic characteristics are obtained by the integrataion over the fluctuating field configurations, which requires an appropriate approximation. A characteristic feature of the Gaussian approximation is the first-order phase transition. We present a renormalization method that takes into account the fourth-order terms of the free energy in the fluctuating field. By the example of the Ising model, we demonstrate that the renormalized Gaussian approximation yields the second-order phase transition, which is observed in experiment.
##### Interaction of strong laser pulses and nanostructures

Kalman Varga, Vanderbilt, Physics, United States

Abstract: Attosecond science has brought “an era of control of the quantum world” enabling scientists to observe the electron dynamics in molecules and solids directly on their natural length (Angstrom) and time (sub-femtosecond) scales. We will present our study of electron and nuclear dynamics induced by strong laser pulses in the framework of the time-dependent density functional theory in real-time and real-space. Several prototypical examples will be used to highlight the correlated electron and nuclear dynamics in strong fields, including Coulomb explosion of clusters, laser-enhanced field emission from nanostructures, and laser-assisted desorption of hydrogen from surfaces of silicon clusters and graphene flakes.
##### On the computer simulations of electron states of ionic nanosystems with extended defects

Yury K. Timoshenko, Voronezh State University, Department of Mathematical and Applied Analysis, Faculty of Applied Mathematics, Informatics, and Mechanics, Russian Federation

Abstract: Yury K. Timoshenko The electron states of nanosystems of ionic compounds AgCl and KCl with extended defects (charged defects, edge dislocations) were under consideration. The semi-empirical tight-binding approximation and different calculation methods were used. The obtained results and efficiency of calculation schemes are discussed. This work is an improvement of our early investigations [1-3]. References 1. Timoshenko Yu K and Shunina V A Surface Sci. 603 2564 (2009) 2. Timoshenko Yu K and Shununa V A Proc. of SPIE 7396 73960T (2009) 3. Timoshenko Yu K Journal of Physics: Conference Series 490 (2014) 012173 doi:10.1088/1742-6596/490/1/012173
##### Computational Simulation of Semiconductor Laser with Optical Injection

gulnaz, Kazakh National Technical university, Department of Applied Physics, Kazakhstan

Abstract: In this paper we study the dynamics of a semiconductor laser with optical injection. The time behaviour of solutions of a system of three coupled nonlinear rate equations, describing the electric field amplitude and the carrier concentration and the phase difference within the resonator, is discussed both qualitatively and numerically. We then concentrate on the periodic orbits that emanate from Hopf bifurcations. Depending on the injection strength and the phase difference two types of oscillations can be found, such as relaxation and periodic oscillations.
##### Numerical solution of a single charged particle under the Spatial Cyclotron Autoresonance Conditions

JESUS GONZALEZ, UNIVERSIDAD DEL MAGDALENA, DEPARTAMENTO DE FISICA, Colombia

Eduardo Orozco, Universidad Industrial de Santander, Physics, Colombia

VALERIY DUGAR ZHABON, Universidad Industrial de Santander, Escuela de Física, Colombia

Abstract: Relativistic dynamics of a single charged particle, accelerated in the cyclotron resonance conditions by TE11P waves in the stationary and inhomogeneous electric and magnetic fields is studied. The magnetic field is increased in the direction of the charged particle movement, so that the electron are in the spatial cyclotron autoresonance conditions. The electrostatic field has a suitable profile, which opposes the diamagnetic force associated with the mirror effect. By using a numerical solution of the Newton-Lorentz equation of motion, with the Leap-frog Boris scheme, a simulation of the system was made. The electric and magnetic fields at particle positions are calculated by using a bilinear interpolation method. It is shown that the resonant interaction of the electron with the electromagnetic wave can be prolonged in the static electric field; as a result the beam gets a larger energy. The numerical experiments are carried out with the TE113, and TE114 microwave fields of 0.1 GHz and 2.45 GHz, respectively, of 6 kV/cm amplitude.
##### Poisson structure on the phase space associated to the hamiltonian dynamics of coupled Korteweg-de Vries type equations

Alvaro Restuccia, Antofagasta University, Physics, Chile

Adrián Sotomayor, Antofagasta University, Mathematics, Chile

Abstract: We present the hamiltonian structures for a wide class of coupled Korteweg-de Vries systems, including the Gear and Grimshaw system that models the strong interaction of internal waves in a stratified liquid and the system of Lou , Tong , Hu and Tang that describes a two layer fluid model. Among the hamiltonian structures of these systems we found new Poisson algebras which define consistent algebras of observables.
##### Modeling of strain in magnetic media by method movable cellular automata

Victor V. Usachev, Far Eastern Federal University, The School of Natural Sciences, Russian Federation

Petr Andriushchenko, Far Eastern Federal University, , Russian Federation

Leonid Afremov, Far Eastern Federal University, Department of theoretical and experimental physics, Russian Federation

Abstract: This paper represents ideas of the application of the method Movable Cellular Automata for modeling the distribution of magnetic particles in an elastic non magnetic matrix. The basic principles of mathematical model and algorithms is shown below. The model is represented as an ensemble of N discrete elements (automata), which identified the characteristic physical properties, rules of interaction and the possible states. Nearby automata capable of forming bonds with each other. Each machine can have 4 neighbors in a square package or 6 neighbors in the dense packing in the 2D model. Bound pair belong the same piece of material, and pair of unrelated automata belong to different fragments. The criterion for the formation of link and the link gap is the distance between the automata. The part of automata uniformly magnetized and interact with each other like magnetic dipoles. The special conditions can be given for certain groups of automata, for example, the lower layer can be fixed, and the upper layer can be influenced by an external force. Critical integration time step dt affects the accuracy of the calculations. Speed of sound in the medium venv possible to use for evaluate dt: dt = 0.25(d/ venv), venv = √ (K/ρ + 4G/3ρ), ρ - density, K - bulk compression modulus, G - shear modulus. In the case of different types of automata the smallest of their respective times is accepted. Offset the automaton from its position at the previous step is calculated according to the Werle’s scheme. At each step of integration the forces, acting on the automaton by the system is calculated, after a new speed and new coordinates of automata is calculated. The resultant force is the sum of the elastic and viscous forces, as well as force of the magnetic interaction automata.
##### Time Dependent Quantum Transport in Nano Structures

Hang XIE, The University of Hong Kong,, Department of Chemistry, Hong Kong

Abstract: In the modern electronics industry, the devices’ scale is in the range of nanometers. In this range the classical description for the electron transport is not suitable and the quantum mechanics will be implemented. We have developed a type of time-dependent-quantum-transport method based on the non-Equilibrium Green’s function theory. This method introduces some auxiliary density matrices to replace the complicated combinations of self-energies and the Green’s functions. So a series of coupling equation of motions are derived instead of solving the differential-integral equations in the traditional quantum dissipation systems, which is very time-consuming. In the practical calculations we use a Lorentzian-Pade expansion scheme to transform the energy integral of spectrum functions into finite summations. An automatic Lorentzian fitting algorithm for the complicated lead spectrum matrices is also developed. For the large systems, we implement a parallel algorithm to deal with the huge memory of the auxiliary density matrices. With this method, we have calculated some dynamic electron response in nano structures. We observe different ‘overshooting’ phenomena in graphene nanoribbons which result in the profiles of the density of states near the Fermi level. We also calculated the transient current in some carbon nanotubes. We find that under a rapid-rising stepwise bias, there exist a lot of high frequency current excitation in these carbon nanotubes.
##### Computational analysis of difenoconazole interaction with soil chitinase

Diana-Larisa VLADOIU, West University of Timisoara, Biology-Chemistry, Romania

Alecu Aurel Ciorsac, Politehnica University of Timișoara, Physical Education and Sport, Romania

Marioara Nicoleta Filimon, West University of Timisoara, Biology-Chemistry, Romania

Vasile Ostafe, West University of Timișoara, Biology-Chemistry, Romania

ISVORAN ADRIANA, West University of Timisoara, Biology-Chemistry, Romania

Abstract: Difenoconazole (DFC) is a fungicide managing to control a comprehensive range of fungi causing diseases of field crops. It usually comes into contact with soil, where it may undergo a variety of transformations and/or causes the loss of soil fertility or environmental damages. Within this study we use a computational approach to investigate the potential binding of DFC to soil chitinases. Computational characterization of the substrate binding sites of Serratia marcescens and Bacillus cereus chitinases using Fpocket tool reflects the role of hydrophobic residues for the substrate binding and the high local hydrophobic density of both sites. Molecular docking study reveals that difenoconazole is able to bind to Serratia marcescens and Bacillus cereus chitinases active sites, the binding energies being comparable. These results are in good agreement with the logP value (4.86) of DFC and published data reflecting that difenoconosale and its stereoisomers are degraded by soil microflora.
##### Ring Dark Solitons in Toroidal Bose-Einstein Condensates

Lauri A. Toikka, University of Turku, Department of Physics and Astronomy, Finland

Kalle-Antti Suominen, University of Turku, Department of Physics and Astronomy, Finland

Abstract: Ring dark solitons (RDSs) are examples of nonlinear quantum states that can be supported by a Bose-Einstein condensate (BEC). Unlike the well-known planar dark solitons, exact analytic expressions for RDSs are not known. We address this problem by presenting exact localised soliton-like solutions to the radial Gross-Pitaevskii equation. To date, RDSs have not been experimentally observed in cold atomic gases, either. To this end, we propose two protocols for their creation in experiments. In addition, we present results regarding the dynamics and stability of RDSs. Under certain trapping potentials, we show that the decay of RDSs into circular arrays of vortex-antivortex pairs can be reversible, but eventually the decay leads to a state with some properties of quantum turbulence.
##### Ferromagnetic/superconducting interface in a hybrid nanoscopic disc

JESUS GONZALEZ, UNIVERSIDAD DEL MAGDALENA, DEPARTAMENTO DE FISICA, Colombia

JOSE BARBA, Universidad Nacional de Colombia, Departamento de Física, Colombia

Abstract: The oscillations exhibited by the magnetization of superconducting nano and mesoscopics structures as a function of the applied magnetic field has attracted attention in the last few years. Nano/Mesoscopic hybrid systems in contact with ferromagnets exhibit interesting transport properties related to the influence of the exchange field on the density of states of clean ferromagnetic structures in contact with superconductor. In this work we study the vortex configurations in a superconducting disk with one central defect inside. The sample is surrounded by a ferromagnetic medium. We calculate the spatial distribution of the Cooper pairs and the phase of the superconducting order parameter and obtain magnetization, free energy and vorticity curves as a function of the magnetic applied field. We found that the first vortex penetration field decrease when a ferromagnetic/superconducting interface is used.
##### The Solutions of the equations of pulse propagation in optical fibers

Yuliya Zhiganshevna Pchelkina, Samara State Aerospace University (SSAU), The Department of applied mathematics, Russian Federation

Ivan Alimenkov, Samara State Aerospace University (SSAU), Applied mathematics, Russian Federation

Abstract: The solution of the equation of propagation of optical pulses in optical fibers is founded in quadratures. The expanded equation of propagation of optical pulses in silica optical fibers is displayed. The localized solution is founded. The decision of the extended equation of propagation of optical pulses in optical fibers for arbitrary functions of the response of a nonlinear medium for external harmonic perturbation is founded in quadratures .
##### Mathematical Modelling on Chemically Reacting Boundary Type Flow under Reduced Gravity Conditions

Wang Hui Ying, University of Poitiers, Pprime-CNRS, France

Abstract: A fire on board spacecraft is, though unlikely owe to strict safety measures, one of the scenarios with the highest damage potential for hardware and crew. In spacecraft, where buoyancy is negligible, the flow is limited to that induced by the ventilation system with a characteristic velocity of 0.2 m/s. Therefore, the flow is expected to be laminar. The most likely scenario of a flame incidentally initiated and spread on board spacecraft is a wall of combustible material adjacent to a low Reynolds number gas flow (forced by the venting system) and exposed to a critical ignition source. When analysing a combustion process, transport equations can be scaled and characteristic parameters can be extracted. In normal-gravity, a buoyantly induced velocity generates an upper limit for all time scales and dominates all forms of transport. All forms of transport slower than natural convection can thus be neglected. In micro-gravity this upper limit disappears and different phenomena, negligible in normal gravity, become relevant and sometimes of great importance. This is the case for all combustion processes related to fire and in particular the estimation of material flammability. Ignition includes solid and gas phase processes. Those related to the solid phase are relevant to the material and can be extrapolated to micro-gravity. The processes related to the gas phase are affected by buoyancy and thus require further attention. Under conditions of reduced gravity where the effects of buoyancy are small, due to the absence of natural convection, time scales associated with combustion processes are much longer and radiation can be the predominant mode of heat transfer even for small diffusion flames. Soot production is enhanced with increasing residence times, further emphasizing the role of radiation. Thus the ability to control soot formation and radiative emissions in nonpremixed combustion is of extreme importance to spacecraft fire safety. The research presented tends to provide fundamental numerical knowledge on co-current, boundary type flame spread. The objective of this work is to use the detail from the numerical simulations to gain further insight into the effects of the fuel injection rate on the soot formation and associated heat flux on a boundary layer diffusion flame. To understand and analyze the flame spread under conditions of reduced gravity, the most important issue is the identification of the flame position/ geometry and associated heat flux. A numerical study is performed to give a quantitative description of soot formation and associated radiation flux with different fuel injection rate. The strong coupling between soot and radiation are studied in three-dimensional flames by solving the Navier-Stokes fluid dynamic equations of elliptic, reacting flow. Binary diffusion coefficients, thermal conductivity and viscosity coefficients are all temperature and composition dependent. Investigations on the effects of fuel (ethylene) injection rate on the geometric characteristics, localisation and stability of the boundary diffusion flame are performed. The computed soot volume fraction is compared with experimental data from a laminar ethylene diffusion flame established over a horizontal flat plate in micro-gravity environment, and a relatively good agreement is attained. An increase of fuel injection velocity results in a significant increase in integrated soot formation. A stable symmetrical flame and counter-rotating vortex in the cross-section is generated even under microgravity conditions with an increase of fuel injection velocity. The extent of the flame in cross-stream plane surpasses significantly the pyrolysis zone. Development of the counter-rotating vortex leads to the formation of the kidney-shape structure on the cross-stream plane for the temperature. This allows soot particle for long residence times before crossing the flame sheet at fuel-rich conditions, and enhances soot formation.
##### A different model for open vehicle routing problem

Semih ONUT, Yildiz Technical University, , Turkey

Mehmet Kamber, Yildiz Technical University, , Turkey

Abstract: Open vehicle routing problem (OVRP) is a generalization of vehicle routing problem (VRP) where vehicles do not require to return to the depot. A new variant of OVRP, heterogeneous fleet balanced open vehicle routing problem (HFBOVRP), is studied in this study. There are more than one types of vehicles which can be assigned to tours, so the problem is a heterogeneous fleet OVRP. First objective is to minimize the total transportation cost. The second objective is to balance the tour lengths. Some theoretical instances are generated and problems are solved by optimization methods.
##### The Analysis of a Reactive Hydromagnetic fluid flow in a channel through a Porous Medium with Convective Cooling.

Hassan Anthony R., Tai Solarin University of Education, Ijagun, Ogun State, Department of Mathematics, Nigeria

Abstract: This study investigates the analysis of a reactive hydromagnetic fluid flow between two parallel plates through a porous medium with convective boundary conditions. Neglecting the consumption of the material which is exothermic under Arrhenius kinetics; it is assumed that the flow system exchanges heat with the ambient following Newton’s cooling law. Analytical solutions of the nonlinear dimensionless equations governing the fluid flow are obtained using the traditional perturbation method and Adomian Decomposition Method (ADM) together with the diagonal Pade Approximation which is used to determine the thermal criticality values as well as bifurcation conditions. Effects of all – important flow properties on the fluid flow are also presented and discussed.

Dr. Haider F Abdul Amir, Universiti Malaysia Sabah, Physics with Electronic, Malaysia

Abstract: Space and ground level electronic equipment with semiconductor devices are subjected to the deleterious effects by radiation. This paper is attempted to present the transient and post-irradiation response of optoelectronic devices to gamma (γ) rays utilizing cobalt-60. In situ measurements were made on the devices under test (DUTs) up to a total dose of 60 krad followed by a post-irradiation not in-flux test for eight hours. Current transfer ratio (CTR) with is the vital merit of the optoelectronic system is found to decrease remarkably with the absorbed dose. This degradation is induced by the interaction of the energetic photons from gamma rays via two main mechanisms. The dominant effect is the mechanism by ionization while the secondary is by displacement. This radiation effect is found to arouse either a permanent or temporarily damage in the DUTs depending on their current drives and also the Total Ionizing Dose (TID) absorbed. The TID effects by gamma rays are cumulative and gradually take place throughout the lifecycle of the devices exposed to radiation. The full damage cascade phenomenon in the DUTs is calculated via the simulation.
##### MR Image Intensity Inhomogeneity Correction

Punga (Visan) Mirela, Dunarea de Jos University, Galati, Faculty of Sciences and Environment, Department of Chemistry, Physics and Environment, Dunarea de Jos University of Galati, Romania

Simona Moldovanu, Dunarea de Jos University, Galati, Faculty of Sciences and Environment, Department of Chemistry, Physics and Environment, Dunarea de Jos , Romania

Luminita Moraru, Dunarea de Jos University of Galati, Faculty of Sciences and Environment, Department of Chemistry, Physics and Environment, Romania

Abstract: MR technology is one of the best and most reliable way to study the brain. Its main drawback is the so-called intensity inhomogeneity or bias field which impairs the visual inspection and the medical proceedings for diagnosis and strongly affects the quantitative image analysis. In order to accurately and effectively restore the original signal, we address here to filtering, bias correction and quantitative analysis of correction. In this report, we used two denoising algorithms; (i) Basis rotation ﬁelds of experts (BRFoE) and (ii) Anisotropic Diffusion (when Gaussian noise, the Perona-Malik and Tukey's biweight functions and the standard deviation of the noise of the input image are considered). Then, the bias field correction is done by using the level set methods. Various metrics, such as Mean Square Error, Peak Signal to Noise Ratio, Normalized Cross-Correlation, Average Difference, Structural Content and Normalized Absolute Error are computed in order to assess the performance of the proposed method.
##### Analytical Models for the Interactions between a Gas Mixture and Quartz Micro Balance Sensor Array

Abbas Abdelaziz, Centre de recherche en Soudage et Contrôle (CSC), , Algeria

Abstract: The aim of this work is to develop analytical models for the thermodynamic equilibrium at the interfaces (gas mixture / Quarz Micro Balance sensor arrays . Differential equations, which describe the change in the partial sensitivities of the sensor array elements depending on the gas mixture components concentrations, and the sensor array parameters, have been developed. Moreover, the responses of the sensor array as a function of the concentrations of the gas mixture components have been modeled.
##### CFD-Exergy Analysis of the Flow in a Supersonic Steam Ejector

Mohamed BOULENOUAR, USTO-MB, Génie Maritime, Algeria

Ahmed OUADHA, USTO-MB, Génie Maritime, Algeria

Abstract: Owing to the growing price of energy and environmental concerns, it is not surprising that the refrigeration industry is continuously looking for alternative solutions to conventional refrigeration systems. A promising alternative is ejector refrigeration systems using steam water as the working fluid. In an ejector refrigeration system, the mechanical compression is replaced by a thermal compression using an ejector. Ejector simplicity (i.e., it has no moving parts) compared to conventional systems and its low cost make it very attractive. Ejector refrigeration systems have the advantage of being simple, easy to install and require less maintenance. In addition, in an ejector system, the compression can be directly obtained without using mechanical power. They can be driven by thermal energy from solar energy or waste heat from industrial processes. However, because of their low coefficient of performance, ejector refrigeration systems are less dominant as compared to conventional refrigeration systems. Therefore, in order to promote the use of these systems, it is essential to improve their performance. This necessarily requires the understanding and improvement of the entrainment process in the ejector. Understanding the entrainment process and the mechanisms responsible for the losses in an ejector is essential in order to construct more efficient ejector refrigeration systems. Methods used to improve its efficiency present complicated issues for both experimental and numerical investigations, especially those aimed at optimisation of the ejector geometry. Experimental techniques usually do not cover the whole operating parameters range of an ejector. These techniques can be successfully complemented using appropriate numerical simulations. The advances in computational fluids dynamics methods permitted the development and implementation of mathematical models for the simulation of complex phenomena encountered in supersonic steam ejectors used in refrigeration systems. The flow field and the associated heat transfer can be described with less restriction using numerical methods offering possibilities for performance improvements of either components or systems. Numerical simulations of the flow can significantly reduce the time and cost in designing of ejectors before the first prototype is produced. In order to use efficiently and effectively the energy resources, recent analyses join quality to quantity of the energy used to realize a given purpose. In other words, associate the second law of thermodynamics which deals with the quality to the first law of thermodynamics dealing with the quantity of energy. The exergy, defined as the maximum useful work that can be extracted from a system in any process which brings the system into equilibrium with its environment, can be used to assess the quality of energy resource. The exergy concept which is based on the second law of thermodynamics is a powerful tool that identifies the location and the magnitude of irreversibilities in a process. The identification of irreversibilities will help to find and design more efficient processes. Exergy analysis presents also the advantage of evaluating meaningful efficiencies that measure how far the efficiency of a system or process deviates from ideality. Further, this technique permits to identify the causes, locations and magnitudes of exergy losses. Generally, exergy analysis uses thermodynamics-based models. These models are based only on mean values of inlet and outlet parameters. Although, these models are simple and cheap in terms of computation, they fall in predicting the fluid motion features. Fluid motion details can be easily obtained using CFD models. In this case, an exergy analysis can be carried using CFD-based models. These models provide all fluid motion details by solving a set of differential equations representing conservation of mass, momentum and energy and turbulence models. Only few studies on exergy analysis have been found in the open literature to the best of the authors’ knowledge, although there are many studies related to ejectors. The current study aims to carry out a CFD-exergy based analysis to assess the main areas of loss in a supersonic steam ejector encountered in ejector refrigeration systems. The governing equations for a compressible flow are solved using finite volume approach based on SST k- model to handle turbulence effects. Flow rates and the computed mean temperatures and pressures have been used to calculate the rate of exergy at the ejector inlets and outlet as well as the exergy losses within the ejector. Furthermore, the exergy efficiency of the ejector has been calculated.
##### Applying Enhancement Filters in the Pre-processing of Images of Lymphoma

Sergio Silva, Federal University of Uberlândia, Faculty of Mechanical Engineering, Brazil

Marcelo Zanchetta do Nascimento, Federal University of Uberlândia, FACOM, Brazil

Valério Batista, Federal University of ABC, CMCC, Brazil

Leandro Neves, São Paulo State University, DCCE, Brazil

Abstract: Lymphoma is a type of cancer that affects the immune system, and is classified as Hodgkin or non-Hodgkin. It is one of the ten types of cancer that are the most common on earth. Among all malignant neoplasms diagnosed in the world, lymphoma ranges from three to four percent of them. Our work presents a study of some filters devoted to enhancing images of lymphoma at the pre-processing step. Here the enhancement is useful for removing noise from the digital images. We have analysed the noise caused by different sources, like room vibration, scraps and defocusing, and in the following classes of lymphoma: follicular, mantle cell and B-cell chronic lymphocytic leukemia. The filters Gaussian, Median and Mean-Shift were applied to different colour models (RGB, LAB and HSV). Afterwards, we performed a quantitative analysis of the images by means of the Structural Similarity Index and the Mean Squared Error. This was done in order to evaluate the similarity between the images. In all cases we have obtained a certainty of at least 75%, which rises to 99% if one considers only HSV. Namely, we have concluded that HSV is an important choice of colour model at pre-processing histological images of lymphoma, because in this case the resulting image will get the best enhancement.
##### Ab initio study of the electronic band structure of the valence band of II-VI C(2×2) reconstructed surfaces

Alberto Rubio, Universidad Autónoma Metropolitana - Azcapotzalco, Ciencias Básicas, Mexico

Daniel Olguín, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Departamento de Física, Mexico

Abstract: The structural and electronic properties of the CdTe(001), CdSe(001), and ZnSe(001) C(2×2) reconstructed surfaces have been studied in the framework of the first principles calculations. To simulate the surface, a seventeen atomic layers slab was used, where also we add four atomic layers of vacuum. Using the experimentally determined lattice parameter as input we have relaxed the internal atomic positions of the outer atomic layers. We show that our model reproduces properly the surface structural parameters, as well as known electronic properties found in bulk for these compounds. Then, we discuss our results of the projected bulk bands and the surface and resonance states found for these surfaces.
##### Basic idea of Corbino-type single-electron transistor

Akira Suzuki, Tokyo University of Science, Physics, Japan

Hisao Taira, Tokyo University of Science, Department of Physics, Japan

Abstract: We have formulated the transmission probability of an electron in a Corbino quantum disk by taking into account charging effect. The confining potential of the Corbino disk has a singularity at the centre of the disk. In order to avoid this singularity problem, we have to reformulate the Schroedinger equation in the Riemannian manifold. The Schroedinger equation describing the electron in Corbino disk must be expressed by introducing a momentum operator reformed by the metric tensor since the original momentum operator is not Hermitian. In order to obtain a Hermitian momentum operator, it requires the deformation of the Hilbert space by introducing a new wave function. This deformation leads to the extra effective potential term in the Schroedinger equation, which depends on the metric, i.e. the geometry of the disk. It should be noted that the charging energy due to the confinement of electrons in the Corbino disk strongly depends on the geometry of the disk. Thus, we discuss the quantum tunneling of a particle confined in the Corbino disk to investigate the effect of both the effective and charging energy of the effective confining potential by using the Wentzel-Kramers-Brillouin (WKB) approximation. We will present that the transmission probability strongly depends on the charging energy. This suggests that the formulated transmission probability is applicable to the analysis of the single-electron transistor.
##### A new type of Wavelets and its relation to Gabor analysis

Ehud Moshe Baruch, Tecnion, Mathematics, Israel

Abstract: We introduce a new type of wavelets. These wavelets interpolate between the classical wavelets frames and a Gabor analysis type frames. We introduce a family of new wavelets which are suitable for localized signals and another family which is suitable for localized frequency.
##### Quantum Magnets and Matrix Lorenz Systems

Stam Nicolis, CNRS-LMPT Tours, Physics Department, France

Abstract: The Landau–Lifshitz–Gilbert equations for the evolution of the magnetization, in presence of an external torque, can be cast in the form of the Lorenz equations and, thus, can describe chaotic fluctuations. To study quantum effects, we describe the magnetization by matrices, that take values in a Lie algebra. The finite dimensionality of the representation encodes the quantum fluctuations, while the non-linear nature of the equations can describe chaotic fluctuations. We identify a criterion, for the appearance of such non-linear terms. This depends on whether an invariant, symmetric tensor of the algebra can vanish or not. This proposal is studied in detail for the fundamental representation of u(2) = u(1) × su(2). We find a knotted structure for the attractor, a bimodal distribution for the largest Lyapunov exponent and that the dynamics takes place within the Cartan subalgebra, that does not contain only the identity matrix, thereby can describe the quantum fluctuations.
##### Multiscale cross-correlation analysis of Shanghai stock market

Dan Li, Wuhan University of Technology, , China

Qingju Fan, School of Science, Wuhan University of Technology, Department of Statistics, China

Abstract: In this paper, we investigate the correlation and cross-correlation behaviors in Shanghai stock market by combining the traditional detrended fluctuation analysis (DFA) and detrended cross-correlation analysis（DCCA） method with moving fitting windows, respectively. The new method can not only show detailed scale exponent properties of non-stationary time series in small and large scale simultaneously, but also provide a more faithful and more interpretable description of series under investigation. Using the moving fitting windows, we find that the correlations in Shanghai B-share is stronger than Shanghai A-share on the whole, and we also show the dynamic long–range cross-correlations behaviors between Shanghai A-share and B-share index series.
##### High order multi-step variational integrators

Dimitrios Vlachos, University of Peloponnese, Department of Informatics and Telecommunications, Greece

Abstract: In this work, we present a new approach to the construction of variational integrators. In the general case, the estimation of the action integral in a time interval [qk,qk+1] is used to construct a symplectic map (qk,qk+1) → (qk+1,qk+2). The basic idea here, is that only the partial derivatives of the estimation of the action integral of the Lagrangian are needed in the general theory. The analytic calculation of these derivatives, give raise to a new integral which depends not on the Lagrangian but on the Euler–Lagrange vector, which in the continuous and exact case vanishes. Since this new integral can only be computed through a numerical method based on some internal grid points, we can locally fit the exact curve by demanding the Euler–Lagrange vector to vanish at these grid points. Thus the integral vanishes, and the process dramatically simplifies the calculation of high order approximations. The new technique is tested for high order solutions in the two-body problem with high eccentricity (up to 0.99) and in the outer solar system.
##### On the Integrability and Linearization of the Second Order Differential Equations via First Integrals and Sundman Transformation

ÖZLEM, ORHAN, , Turkey

Teoman Özer, Istanbul Technical University, Civil Engineering , Turkey

Abstract: In this study, we consider the ordinary differential equations of the form x ̈+a2(t,x)x ̇2+a1(t,x)x ̇+ a0(t, x) = 0. The second order ordinary differential equations in this form have first integrals of the form A(t, x)x ̇ +B(t, x), λ-symmetries and integrating factors. The purpose of this study is to characterize the second order differential equations and integrate these equations by using first integrals and λ-symmetries. We use an algorithm to calculate these symmetries. Furthermore, we characterize equations that can be linearized by means of nonlocal transformation which is called Sundman transformation. We can apply these algorithm and transformation to nonlinear fin equation where has thermal conductivity and heat transfer coefficients, which are considered as functions of temperature. We obtain first integrals for different coefficients of fin equation. Finally, we compare the results obtained by different methods.
##### Function projective synchronization of two identical discrete-time hyperchaotic systems

Suping Qian, Changsha institute of technology, School of Mathematics and statics, China

Abstract: A function projective synchronization is defined in discrete-time dynamical systems, in which the drive and response state vectors evolve in a proportional scaling function matrix. Based on backstepping design with three controllers, a systematic, concrete and automatic scheme is developed to investigate the function projective synchronization of two identical discrete-time hyperchaotic systems. Numeric simulations are used to verify the effectiveness of our scheme.
##### Cluster analysis of word frequency dynamics

Yulia Maslennikova, Kazan federal university, Institute of Physics, Russian Federation

Inna, Kazan Federal University, , Russian Federation

Abstract: This paper describes the analysis and modelling of word usage frequency time series. During one of the previous studies, an assumption was put forward that all word usage frequencies have uniform dynamics that is near to a form of a Gaussian function. This assumption can be checked using frequency dictionaries of Google Books Ngram database. Google Books database includes 5.2 million books published between 1500 and 2008. The corpus contains over 500 billion of American English, British English, French, German, Spanish, Russian, Hebrew, and Chinese words. In this research, time series of word usage frequencies were clustered using a Kohonen's neural network. The similarity between input vectors was estimated using both algorithms: the Spearman rank correlation coefficients and the maximum likelihood method. As a result of the neural network training procedure, more than ten different forms of time series were found. They correspond to dynamics of word usage frequencies from the word birth to the word death. Сonsequently, different groups of word forms have different dynamics of word usage frequency variations.
##### Evaluation of outliers in acquired brain MR images

Simona Moldovanu, Dunarea de Jos University, Galati, Faculty of Sciences and Environment, Department of Chemistry, Physics and Environment, Dunarea de Jos , Romania

Luminita Moraru, Dunarea de Jos University of Galati, Faculty of Sciences and Environment, Department of Chemistry, Physics and Environment, Romania

Punga (Visan) Mirela, Dunarea de Jos University, Galati, Faculty of Sciences and Environment, Department of Chemistry, Physics and Environment, Dunarea de Jos University of Galati, Romania

Abstract: Pre-processing is an important stage in the analysis of magnetic resonance images (MRI), because the effect of specific image artefacts, such as intensity inhomogeneity, noise and low contrast can adversely affect quantitative image analysis. Image histogram is a useful tool to analyse MR images as it allows a close relationship with important image features as contrast and noise. The noise and variable contrast are elements that locally modify the quality of images. The key issue of this study derives from the fact that spatial histogram can contain outliers indicating corrupted image information through the disorder of the bins. These aberrant errors should be excluded from the studied data sets. Here, the outliers are evaluated by using rigorous methods based on the probability theory and Chauvenet (CC), Grubbs (GC) and Prices (PC) criteria. In order to check the quality of the MR images, the Minkowsky (MD), Euclidean (ED) and cosine (CD) distance functions were used. They act as similarity scores between histogram of the acquired MRI and processed image. This analysis is necessary because, sometimes, the distance function exceeds the co-domain due to the outliers. In this paper, 32 MRIs are tested and the outliers are removed so that the distance functions generate uncorrupted and real value.
##### Biochemical and histopathological response of liver rats induced by antilipolytic drugs.

Nofal, Menoufia university Egypt, Zoology, Egypt

Abstract: The present study provides evidence that 1,5-Bis (3,5-Dimethylpyrazol-1-yl)-3-oxapentane-diacetatocopper has an antidiabetic effect, as hypoglycemic agent and as antilipolytic agent, but with many abnormalities. Animals treated with 1,5-Bis (3,5-Dimethylpyrazol-1-yl)-3-oxapentane-diacetatocopper revealed many histopathological, biochemical, histochemical and immunohistochemical alterations on the liver. Sera of animals treated with 1,5-Bis (3,5-Dimethylpyrazol-1-yl)-3-oxapentane-diacetatocopper revealed a significant increase in ALT and AST. The present study revealed many histopathological alterations on the liver; inflammatory infiltration, marked vacuolated cytoplasm in cells, congestion of blood vessels, hemorrhage, pyknotic cells and binucleated cells, as well as, some of the degenerated cells showed karyorhexis, pyknosis and area of necrosis, as well as, significant decrease in proteinic contents on the liver. the present study revealed negative expression of Bcl-2 proteins in all hepatic cells of control rats and positive expression of Bcl-2 proteins in hepatic cells of treated rats with 1,5-Bis(3,5-Dimethylpyrazol-1-yl)-3-oxapentane-diacetatocopper.
##### RELAXED ELASTIC LINES OF SECOND KIND IN SEMI-DUAL SPACE

Gülşah Aydın, süleyman demirel university, , Turkey

A.Ceylan ÇÖKEN, University of Süleyman Demirel, Department of Mathematics, Turkey

Abstract: Theory of elasticity is a topic that keeps improving by using on many fields such as geometry, physics, chemistry and engineering. Energy density is given as some functions of curvature and torsion. If the curve of the α will be an external for the variation problem that minimizes the value of energy density; then this curve is called as relaxed elastic line. The relaxed elastic line on an oriented surface is considered as a model of DNA molecule. In this study, we worked on the second type relaxed elastic lines on the semi-dual spaces which has an important point on kinematic and Einstein’s relativity theory. We also obtained boundary conditions for this type of curves. Moreover, the minimization problem of the energy which occurs with an applied force on an elastic line was discussed. Then, we researched the formed potential energy due to the applied force. Also, during the calculation of the potential energy on the elastic line, the amount of the potential energy for unit length of the elastic line was used. Afterwards, by integrating that amount, total potential energy calculated. So, we study to make a contribute both Einstein’s relativity theory and kinematic.
##### Computer simulation of electron states of the isoelectron substitutional impurities in quantum wire of KCl with edge dislocation

Yury K. Timoshenko, Voronezh State University, Department of Mathematical and Applied Analysis, Faculty of Applied Mathematics, Informatics, and Mechanics, Russian Federation

Valentina Shunina, Voronezh State Technical University, Department of Information Technologies and Computer Safety, Russian Federation

Alexander I. Shashkin, Voronezh State University, Department of Mathematical and Applied Mathematics, Faculty of Applied Mathematics, Informatics, and Mechanics, Russian Federation

Abstract: Yury K. Timoshenko, Valentina A. Shunina, Alexander I. Shashkin The electron states of quantum wires KCl, KCl:Br, and KCl:J with an edge dislocations were investigated. The tight-binding semi-empirical band approximation and non-empirical cluster approach were used. Semi-empirical calculations were carried out in framework of model [1, 2]. Besides, the algorithms for electronic levels calculations of polar nanosystems with the partial self-consistency [3] were used. Using results of computer modelling we discuss the problem of localization of electronic states near the line of edge dislocation. References 1. Timoshenko Yu K and Shunina V A 2005 Phys. Stat. Sol. (c) 2 1788 2. Timoshenko Yu K and Shunina V A 2009 Surface Sci. 603 2564 3. Timoshenko Yu K Journal of Physics: Conference Series 490 (2014) 012173 doi:10.1088/1742-6596/490/1/012173
##### The generalized uncertainty principle as quantum gravitational friction

Pedro Bargueno, Universidad de los Andes, Física, Colombia

Abstract: The Generalized Uncertainty Principle gives place to deformed commutation relations which are linear or quadratic in particle momenta. In this talk we show that, in the linear case, which corresponds to double special relativity theories, this deformation is equivalent to a gravitationally-induced damping process in an Ohmic environment at zero temperature. Therefore, both minimum length and maximum momentum give place to quantum gravitational friction.
##### Wormholes Threaded by Chiral Fields

Theodora Ioannidou, Aristotle University of Thessaloniki, Department of Mathematics, Physics and Computational Sciences, Greece

Abstract: We consider Lorentzian wormholes with a phantom field and chiral matter fields. The chiral fields are described by the non-linear sigma model with or without a Skyrme term. When the gravitational coupling of the chiral fields is increased, the wormhole geometry changes. The single throat is replaced by a double throat with a belly inbetween. For a maximal value of the coupling, the radii of both throats reach zero. Then the interior part pinches off, leaving a closed universe and two (asymptotically) flat spaces. A stability analysis shows that all wormholes threaded by chiral fields inherit the instability of the Ellis wormhole.
##### The peculiarities of magnetization processes in layered high-temperature superconductors with ferromagnetic defects under applying of transport current and external magnetic field

Vladimir Kashurnikov, National Research Nuclear University MEPhI, Department of Physics of the Solid State and Nanoscaled Systems (№70), Russian Federation

Anastasiia Maksimova, National Research Nuclear University MEPhI, Department of Physics of the Solid State and Nanoscaled Systems (№70), Russian Federation

Igor Rudnev, National Research Nuclear University MEPhI, Department of Physics of the Solid State and Nanoscaled Systems (№70), Russian Federation

Abstract: The Monte-Carlo method was used for study of magnetization processes in 2D high-temperature superconductors with internal ferromagnetic defects under application of transport current and external dc magnetic field. The magnetic part of the pinning force was calculated as an energy of interaction of magnetic dipole of ferromagnetic defect with the field of Abrikosov vortex. The self-consistent interaction between magnetic moments of the particles and the vortex system was taken into account. The current-voltage characteristics (VACH) in presence of external dc magnetic field were calculated. A novel S-type of current-voltage characteristics of the superconductor/ferromagnet system in external magnetic field was demonstrated. It was shown that the S-type nonlinearity is due to the local reversal magnetization of magnetic particles by the field of vortices. The effect of temperature, value of external magnetic field and defect concentration on magnetization processes were analyzed. The H-T phase diagram which demonstrates the region of existence VACH nonlinearity was obtained. Vortex configurations arising during magnetization were obtained. The conditions for electromagnetic generation at the region of nonlinearity were found and the frequency of such a generation was estimated.
##### Sharpening technique to decrease comb passband deviation

Gordana Jovanovic Dolecek, Institute INAOE, Electronics, Mexico

Abstract: This paper deals with the wideband compensation of the comb decimation filter pass-band droop using the sharpened technique. The sharpening technique is introduced for simultaneous improvements of both the pass-band and stop-band of a linear-phase FIR (Finite Impulse Response) filters. The technique uses the amplitude change function (ACF) which is a polynomial relationship between the amplitudes of the sharpened and the prototype filters. One-stage and two two-stage sharpened-based compensated combs are proposed taking the cascaded narrowband and wideband compensators. The choice among structures offers a trade-off between the complexity and the pass-band improvement. The resulting pass-band characteristic is almost flat for the one-stage decimation comb structure. In the two-stage structures the sharpening is performed only at the second stage and as a result the pass-band deviation is slightly increased but it is less than 0.05db. The proposed structures are compared with the compensated combs introduced in the open literature.
##### COMPTON OPERATOR IN QUANTUM ELECTRODYNAMICS

Enrique Ortega Muñoz, Instituto Politécnico Nacional, Ingenieria Quimica, Mexico

Abstract: In the frame in the quantum electrodynamics exist four basic operators; the electron self-energy Σ, vacuum polarization Π, vertex correction Λ, and the Compton operator, denoted by C. The first three operators are very important by its relation with renormalized ans Ward identity. However, the Compton operator has equal impotance, but without divergence, and Little attention has been given it. We have calculated the Compton operator and obtained the closed expression for it in the frame of dimensionally continuos integration and hypergeometric functions.
##### Time adaptive phase fitted variational integrators using a space-time geodesic approach

Odysseas Kosmas, University of Erlangen-Nuremberg, Mechanical Engineering, Germany

Abstract: In this work, we investigate the use of space-time geodesic approach of classical mechanics in order to derive time adaptive phase-fitted variational integrators. The proposed technique is tested in systems with separable Lagrangians. Towards this end, at first we unfold the standard Euler-Lagrange system to its space-time manifold and, then, we recast it as a geodesic problem in which the potential energy part is absent. Preliminary simulated results on the Kepler problem (without optimizing the choice of step sizing) show that, one can use the space-time geodesic formulation to generate an adaptive scheme that still preserves some underlying geometric structure.
##### Study on EP resin modified NANO-SILICA and mechanical properties of composites

Ju Chunhua, Harbin Institute of Technology, Material Science, China

Abstract: Owing to their light, excellent mechanical properties and a series of other advantages, and obviously advantage of lightweight of polymer composite Liquid Oxygen (LOX) tank, composite materials have been critical and broadly studied in the development of the launch vehicles of next generation, especially in the propelling system. In this present research, prepared polymer matrix of composite was modified, and influence of composite materials interface and mechanical properties at room and low temperature was studied. We want to prepare composite with good properties and wish to the lay a good foundation of this area for the following researchers. Curing kinetics and mechanism of polymer matrix are exmined by means of TG-DSC and FTIR. Using the result of dynamic contact angle, micro-debonding and mechanical properties of composite materials and matrix, influence of modified with nano-silica resin to composite at room and low temperature, and low temperature to matrix, interface, mechanical properties of composite, are exmined.
##### Modeling deformation and fracture of Ni nanofilm

Mohammed, Physics department, Faculty of science, Menoufia university, Egypt, Physics , Egypt

Abstract: Morse potentials were employed to carry out three dimensional molecular dynamics simulations. A computer experiment is performed at a temperature corresponding to 300K and 1000K. MD simulation used to investigate the effect of cross-sectional area of Ni nanofilm on the nature of deformation and fracture. The engineering stress–time diagrams obtained by the MD simulations of the tensile specimens of these Ni nanofilms show a rapid increase in stress up to a maximum followed by a gradual drop to zero when the specimen fails by ductile fracture. The feature of deformation energy can be divided into four regions: quasi-elastic, plastic, flow and failure. The results showed that breaking position depended on the nanofilm length cross-sectional area.
##### Modeling and simulation of a batch reactor with automatic temperature control applied to fine chemicals.

Hachama Kamel, Laboratory valorization on natural substances, Khemis Miliana University, Sciences and Technology, , Algeria

BEZZINA MOHAMED, University of Khemis-Miliana, , Algeria

ZOUANTI MUSTAPHA, University of Khemis-Miliana, , Algeria

Abstract: The approach presented in this work aims to develop a methodology for modeling of a batch chemical reactor, which takes place in these, two consecutive reactions, irreversible, highly exothermic and in a homogeneous liquid phase. We are interested in this study of the thermal stability of the system and the evolution of concentration of the reactants and products as a function of time. This methodology is the first step to develop a comprehensive model that takes into account, the coupling between the equations of material balance, energy balance equations and the system of control and regulation. The results are used to study the operation of the reactor and predict its behavior in various situations, and allow us to optimize its operation to increase efficiency and ensure quality products, while avoiding thermal runaway.
##### Two-Step Syntheses of 3-Methyl-1,2,4-Benzotriazine

Hachama Kamel, Laboratory valorization on natural substances, Khemis Miliana University, Sciences and Technology, , Algeria

KHODJA MOHAMED, University Blida 1, , Algeria

ZOUIKRI MOHAMED, University Blida 1, , Algeria

Abstract: The reduction cyclization of nitrophenylhydrazones has been widely used to afford nitrogen-based heterocyclic molecules. To name but a few, there are 4H-pyrazolo-[1,5-a] benzimidazoles, ethyl-(1,2,4-benzotriazine-3-yl) acetate and its derivatives, benzo[1,2-b:5,4-b]bis(1H)-imidazo[1,2-b]pyrazoles. Several pathways for the synthesis of 1,2,4-Benzotriazines have been reported. In continuation of our work on heterocyclic synthesis, we developed short pathways to produce 2-Methyl-1,2,4-benzotriazine compounds and some of its derivatives, were prepared in yields to (50-70%), via a reductive cyclization by a PtO2-catalyzed hydrogenation, of the corresponding 2-nitrophenylhydrazones of the pyrivic acid. The later compounds were obtained in yields higher than 90% by reacting 2-nitrophenylhydrazines with sodium pyruvate salt.
##### High Performance Computing based on the CPU and GPU systems for the Quantum Trajectory Method

Joanna Wiśniewska, Military University of Technology, Faculty of Cybernetics, Poland

Marek Sawerwain, University of Zielona Góra, Institute of Control and Computation Engineering, Poland

Wiesław Leoński, University of Zielona Góra, Quantum Optics and Engineering Division, Institute of Physics, Poland

Abstract: Nowadays, the dynamic progress in computational techniques allows for a development of various methods which offer significant speed-up of computations, especially those related to the problems of quantum optics and quantum computing. In this work, we propose computational solutions which re-implement the Quantum Trajectory Method (QTM) algorithm in the modern parallel computation environment where multi-core CPUs and modern many-core GPUs can be used. As the effect, more effective computational routines - than those applied in other commonly used packages, such as Quantum Optics Toolbox for Matlab or QuTIP for Python - are obtained. The QTM algorithm is still widely used method for examination of open quantum systems, but its application requires a special attention at its reimplementation in modern hardware systems, especially those comprising GPUs. The main problem in application of QTM is to focus on a proper definition and use of computational grid (or computational index space), and it is especially important in implementation for CUDA technology discussed here. The proposed implementation of QTM is also a hybrid method, in which power of traditional CPU and modern GPU is combined to obtain highly effective and precise routines for simulation of the dynamics of open quantum systems. It should be pointed out that discussed here solution, based on GPUs, is about 20 or even 50 times faster than other implementations of QTM prepared for traditional serial CPU systems. The calculation speed-up in proposed implementation depends on numerical method chosen for solving ordinary differential equations (Runge-Kutta and Backward Differentiation methods were simultaneously applied in the proposed solution).
##### A new group theoretical approach to ultrahyperbolic general relativity

Evangelos Melas, Technological Educational Institution of Patras, Department of Management, Greece

Abstract: The Bondi$-$Metzner$-$Sachs (BMS) group $B$ is the common asymptotic group of all asymptotically flat (lorentzian) space$-$times, and is the best candidate for the universal symmetry group of General Relativity (G.R.). $B$ admits generalizations to real space$-$times of any signature, to complex space$-$times, and supersymmetric generalizations for any space$-$time dimension. Generalizations of $B$ are important both for quantum and classical G.R.. %and general relativity. In 1939 Wigner laid the foundations of special relativistic quantum mechanics and relativistic quantum field theory by constructing the Hilbert space strongly continuous unitary irreducible representations (IRs) of the (universal cover) of the Poincare group $P$. In a quantum setting the universal property of $B$ for G.R. make it reasonable to attempt to lay a similarly firm foundation for quantum gravity by following through the analogue of Wigner's programme with $B$ replacing $P$. With this motivation McCarthy constructed the IRs of $B$ some time ago. $B(2,2)$ is the generalization of $B$ appropriate to the to the ultrahyperbolic signature' ($+$,$+$,$-$,$-$) and asymptotic flatness in null directions. Following a historical introduction and some motivating remarks on the correspondence of the IRs of $B$ and its generalisations to the solutions of the Einstein equations in all signatures as well as in complex space$-$times we continue this programme by introducing a new group $\mathcal H \mathcal B$ in the group theoretical study of ultrahyperbolic G.R. which happens to be a proper subgroup of $B(2,2)$. We give the first general results on the representation theory of $\mathcal H \mathcal B$. In particular it is proved that {\it all} little groups of $\mathcal H \mathcal B$ are compact and that the Wigner$-$Mackey's inducing construction is exhaustive despite the fact that $\mathcal H \mathcal B$ is not locally compact in the employed Hilbert topology. The significance of these results for this approach to quantum gravity as well as their repercussions for other approaches to quantum gravity is analysed. Finally relation to other work is outlined.
##### Application of nonlinear methods to the study of ionospheric plasma

Alexander Chernyshov, Space Research Institute, , Russian Federation

Mikhail Mogilevsky, Space Research Institute, , Russian Federation

Boris Kozelov, Polar Geophysical Institute, , Russian Federation

Abstract: Most of the processes taking place in the auroral region of Earth's ionosphere are reflected in a variety of dynamic forms of the aurora borealis. In order to study these processes it is necessary to consider temporary and spatial variations of the characteristics of ionospheric plasma. Most traditional methods of classical physics are applicable mainly for stationary or quasi-stationary phenomena, but dynamic regimes, transients, fluctuations, self-similar scaling could be considered using the methods of nonlinear dynamics. Therefore, development of methods for studies of characteristics of nonlinear processes in open dissipative systems is now actual problem, in particular, to study the Earth's ionosphere. Special interest is the development of the methods for describing the spatial structure and the temporal dynamics of auroral ionosphere based on the ideas of percolation theory and fractal geometry. The fractal characteristics (the Hausdorff fractal dimension and the index of connectivity) of Hall and Pedersen conductivities are used to the description of fractal patterns in the ionosphere. To obtain the self-consistent estimates of the parameters the Hausdorff fractal dimension and the index of connectivity in the auroral zone, an additional relation describing universal behavior of the fractal geometry of percolation at the critical threshold is applied. Also, in this work Tsallis statistics is used to study the properties of the auroral ionosphere.
##### Geometry Design of an Electron Dump with Simple Geo

Mert Sekerci, Süleyman Demirel University, , Turkey

Fatma OZTURK, Istanbul University, Department of Physics, Turkey

Nurgul HAFIZOGLU, Istanbul University, Department of Physics, Turkey

Suat OZKORUCUKLU, Istanbul University, Department of Physics, Turkey

Abstract: The goal of this study is to design a geometry of an electron dump with Simple Geo code which is a freeware product and provides the ability of designing complex geometric systems easily. Also, Simple Geo can output the designed geometry in many different formats. Desired design of the electron dump is to stop the 40 – 42 MeV electron beams. To reach this aim, requested geometric design with the possible material was done with Simple Geo and a FLUKA output format file created to run the simulations in FLUKA code.
##### Residual symmetries and B\"acklund transformations

Lou SenYue, Ningbo University , Physics Department, China

Abstract: It is known that for a given truncated Painlev\'e expansion of an arbitrary nonlinear Painlev\'e integrable system, the residue with respect to the singularity manifold is a nonlocal symmetry. It is proved that the residual symmetries can be localized to Lie point symmetries for suitable prolonged systems. For the Korteweg-de Vries equation, the $n^{th}$ Binary Darboux transformation is re-obtained by the Lie point symmetry approach companied by the localization of the $n$-fold residual symmetries.
##### Born Oppenheimer Molecular Dynamics calculation of the nuO-H IR spectra for acetic acid cyclic dimers

Benmalti Mohamed el Amine, Université abdelhamid ibn badis-mostaganem, Chimie, Algeria

Abstract: Both ab initio molecular dynamics simulations based on the Born-Oppenheimer approach calculations and a quantum theoretical model are used in order to study the IR spectrum of the acetic acid dimer in the gas phase. The theoretical model is taking into account the strong anharmonic coupling, Davydov coupling, multiple Fermi resonances between the first harmonics of some bending modes and the rst excited state of the symmetric combination of the two O-H modes and the quantum direct and indirect relaxation. The IR spectra obtained from DFT-based molecular dynamics is compared with our theoretical lineshape and with experiment. Note that in a previous work we have shown that our approach reproduces satisfactorily the main futures of the IR experimental lineshapes of the acetic acid dimer[Mohamed el Amine Benmalti, , Paul Blaise, H. T. Flakus, Olivier Henri- Rousseau, Chem Phys, 320(2006) 267-274.].
##### KP Predication by Using mixture of MLP-experts

Mina Ghoorchian, HIRBODAN ARIANA Professional Engineers Training Institute, , Iran (Islamic Republic of)

Abstract: For space weather phenomenon that is searched always,in the paper ,solar activity effects on power system are considered.According research,solar plasma can have destructive effect on power system. such as saturation of the transformer core,extra harmonic that can [switch out of service such quebec blackout of march 13,1989 witch left the whole province without power for over 9 hours].also there are major problems in restoring power after grid collapse such as cold load pick-up,no spare transformer. The solar wind driven magnetosphere is a complex dynamical system with highly nonlinear and chaotic behavior. A large number of studies have been carried out to provide appropriate dynamical models of magnetosphere, and to predict various geomagnetic indices, e.g. Dst storm time index and AE auroral electrojet index. But the most popular indicator of geomagnetic disturbances, the Kp index, which is used mainly in warning and alert systems for power system, has not been considered as much. Kp is used the same as alarm in 5 level, Kp =9,8,7,6 and 5. Kp is predicted and compared with mixture of MLP - experts this method can work similar an alarm intelligent system for power system.
##### FUZZY INTERACTIVE GOAL PROGRAMMING METHOD FOR SUM OF LINEAR FRACTIONAL MULTI-OBJECTIVE PROGRAMMING PROBLEMS

Pitam Singh, Motilal Nehru National Institute of Technology Allahabad, Mathematics, India

Abstract: This paper is concerned with an efficient algorithm for solving the multi-objective program of sum of linear fractional functions whose constraints are linear. The algorithm proposed here for the solution of sum linear ratios multi-objective programming (SOLR-MOP)is integrate the merits of three approaches , (i) fuzzy approach, (ii) fuzzy goal programming, (iii) Interactive programming. Each membership function of fuzzy objective is approximated into linear function by using first order Taylor theorem about the vertex of the feasible region where the objective function has maximum value. Then the resulted approximated linearized membership function may be used for the formulation of interactive fuzzy goal programming.The efficiency of the method is measured by numerical a example.
##### Some More Solutions of Burgers’ Equation

Mukesh Kumar, MNNIT Allahabad, Mathematics, India

Abstract: In the present study, one-dimensional Burgers' equation has been studied. An attempt has been made to search Lie symmetries and construction of exact solutions of viscous Burgers' equation through Lie group properties. The symmetry generators are used for constructing Lie group of point symmetries with commuting infinitesimal operators which lead the governing partial differential equation to ordinary differential equation. To the best of author’s knowledge, exact solutions of Burgers' equation obtained in terms of Bessel functions are appearing very first time in the existing literature. Effect of various parameters through graph is shown and discussed physically.
##### Symbolic computations by the Wick's theorem.

Igor V Beloussov, Institute of Applied Physics, Academy of Sciences of Moldova, Theoretical Physics, Moldova, Republic of

Abstract: The new formulation of Wick's theorem that allows one to present the vacuum or thermal averages of the chronological product of an arbitrary number of field operators as a determinant (permanent) of the matrix is proposed. Each element of the matrix is the average of the chronological product of only two operators. This formulation is extremely convenient for practical calculations in quantum field theory and statistical physics by the methods of symbolic mathematics using computers.
##### Asymptotic solution of natural convection in a uniformly Joule-heating shallow cavity

xiaohui zhang, Soochow University, , China

Abstract: The steady laminar two-dimensional Joule heating natural convection is investigated using asymptotical analysis, the fluid is in a rectangular cavity, the direct current contributes heat for heating the process medium by a pair of plate electrodes, the top wall is cooled with atmosphere and all the other walls are kept thermally insulated. The asymptotic solution is obtained in the core region in the limit as the aspect ratio, which is defined as the ratio of the vertical dimension of cavity to the horizontal dimension of cavity, goes to zero. The numerical experiments are also carried out to compare with the asymptotic solution of the steady two-dimensional Joule heating convection. The asymptotic results indicate that the expressions of velocity and temperature fields in the core region are valid in the limit of the small aspect ratio.
##### Gravitational radiation as radiation same level of electromagnetic and his generation of pulsed high-current discharge

Fisenko, Rusthermosynthesis JSC, , Russian Federation

Abstract: The notion of gravitational radiation as a radiation of the same level as the electromagnetic radiation is based on theoretically proved and experimentally confirmed fact of existence of stationary states of an electron in its gravitational field characterized by the gravitational constant K = 10^42G (G is the Newtonian gravitational constant) and unrecoverable space-time curvature Λ. This paper gives an overview of the authors' works [1, 2, 3, 4], which set out the relevant results. Additionally, data is provided on the broadening of the spectra characteristic radiation. The data show that this broadening can be explained only by the presence of excited states of electrons in their gravitational field. What is more, the interpretation of the new line of X-ray emission spectrum according to the results of observation of MOS-camera of XMM-Newton observatory is of interest. The given work contributes into further elaboration of the findings considering their application to dense high-temperature plasma of multiple-charge ions. This is due to quantitative character of electron gravitational radiation spectrum such that amplification of gravitational radiation may take place only in multiple-charge ion high-temperature plasma.
##### The origin of self-focusing effect in terahertz quantum cascade lasers

Jian Wang, Beijing Jiaotong University, Department of Physics, China

Yanfang Li, Institute of Applied Physics and Computational Mathematics, , China

Ning Yang, Institute of Applied Physics and Computational Mathematics, , China

Abstract: The terahertz quantum cascade lasers (THz-QCLs) are the compact and coherent terahertz light source based on the inter-subband transition and resonant tunneling of carriers in semiconductor superlattice. In the recent studies on tapered THz-QCLs, it was found that the self-focusing effect in the active region of the devices may cause the abnormal increase of the far-field divergence of the laser beam. By simulating the propagation of optical model in QCL waveguide and considering both the nonlinearity effect and thermal accumulation in the active region, we demonstrate that the refractive index change caused by the third-order nonlinearity of multi-quantum-wells in active region is the key reason for the self-focusing in THz-QCLs. This result indicates that the nonlinear effect has great impact on the beam quality of QCLs which must be carefully considered in applications of THz-QCLs, such as the TH- imaging systems.
##### A simple model of hysteresis behavior using spreadsheet analysis

Andrea Ehrmann, Niederrhein University of Applied Sciences, Faculty of Textile and Clothing Technology, Germany

Tomasz Blachowicz, Silesian University of Technology, , Poland

Abstract: Hysteresis loops occur in many scientific and technical problems, especially as field dependent magnetization of ferromagnetic materials, but also as stress-strain-curves of materials measured by tensile tests including thermal effects, liquid-solid phase transitions, in cell biology or economics. While several mathematical models exist which aim to calculate hysteresis energies and other parameters, here we offer a simple model for a general hysteretic system, showing different hysteresis loops depending on the defined parameters. The calculation which is based on basic spreadsheet analysis plus an easy macro code can be used by students to understand how these systems work and how the parameters influence the reactions of the system on an external field. Importantly, in the step-by-step mode, each change of the system state, compared to the last step, becomes visible. The simple program can be developed further by several changes and additions, enabling the building of a tool which is capable of answering real physical questions in the broad field of magnetism as well as in other scientific areas, in which similar hysteresis loops occur.
##### Anatomy of Demagnetizing and Exchange Fields in Magnetic Nano-Dots Influenced by 3D Shape Modifications

Tomasz Blachowicz, Silesian University of Technology, , Poland

Andrea Ehrmann, Niederrhein University of Applied Sciences, Faculty of Textile and Clothing Technology, Germany

Abstract: Hysteresis loops of 3D ferromagnetic permalloy nano-half-balls (dots) with 100 nm base diameter have been examined by means of LLG micromagnetic simulations and finite element methods. Tests were carried out with two orthogonal directions of the externally applied field at 10 kA/(m.ns) field sweeping speed. The comparison of samples with different 3D modifications at the sub-10nm scale, accessible by nowadays lithographic techniques, enables conclusions about different mechanisms of competition between demagnetizing and exchange fields. Design paradigms provided here can be applied, e.g., in bit-patterned media used as novel magnetic storage systems.
##### Ideal Fluids, the Quark Gluon Plasma, and Hadronic Gases

Dr. Nasser S Demir, Kuwait University, Physics, Kuwait

Abstract: Ultrarelativistic heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) are believed to have created a state of matter known as the Quark Gluon Plasma (QGP), where quarks are thought to be deconfined in a "quark soup." The QGP thought to have been produced at RHIC and LHC has properties resembling those of an ideal fluid. As such, significant interest arose in quantifying the shear viscosity coefficient and the shear viscosity to entropy density ratio of the QGP. A conjecture on the lower bound of the shear viscosity to entropy density ratio of the "most perfect fluid" is that its minimum value should be $\frac{\hbar}{4 \pi k_B}$. However, in order to quantify the shear viscosity of the QGP, it is necessary to quantify the shear viscosity of the hadronic phase of a heavy ion reaction. As the QGP cools and decays, a hadronic gas is formed, which is assumed to be highly viscous. Several calculations of the shear viscosity exist for hadronic systems, but a systematic comparison between two disparate methods is necessary for validation and for testing the limits of any method involved. A systematic comparison between two methods (the Green-Kubo and Chapmann-Enskogg methods) is performed for two sets of hadronic systems. One system involves a system of chiral pions interacting with an energy dependent cross section and the other system involves a pion rho mixture. Differences are highlighted between the results.
##### Persistent Dirac fermion state on Si surface

Kehui Wu, Institute of Physics, Chinese Academy of Sciences, , China

Abstract: Recently, silicene, a single sheet of Si atoms arranged in honeycomb lattice with sp2 bonding, has been proposed and successfully fabricated on Ag(111) and other substrates in ultrahigh vacuum. The existence of Dirac fermion in silicene on Ag(111) surface with (√3×√3)R30° superstructure has been proven by the observation of linear energy-momentum dispersion and quasiparticle chirality by scanning tunneling microscopy (STM) and spectroscopy (STS). In addition to monolayer silicene sheet, we reported the existence of “multilayer silicene” film on Ag(111), comparable to multilayer graphene. In this talk we will present our new results on silicene, mainly concerning the detailed features of the electronic states, and electronic states on multilayer silicene film. For silicene film with increasing thickness above 30 monolayers (ML), the film always exhibits a (√3×√3)R30° honeycomb superstructure on the surface, which has never been observed on bulk Si(111) surface. More interestingly, we observed persistent Dirac Fermion state on the surface, independent on the film thickness. The origination of this Dirac fermion state is discussed
##### Study of basic structural, sorption and diffusion properties of novel microporous materials based on a combination of molecular simulation methods

Stylianos Karozis, National Center for Scientific Research Demokritos, Environmental Research Laboratory, Greece

Georgia Charalambopoulou, National Center for Scientific Research Demokritos, Environmental Research Laboratory, Greece

Nikolaos Papadimitriou, National Center for Scientific Research Demokritos, Environmental Research Laboratory, Greece

Theodore Steriotis, National Center for Scientific Research Demokritos, Institute of Nanoscience and Nanotechnology, Greece

Michael Kainourgiakis, National Center for Scientific Research Demokritos, Environmental Research Laboratory, Greece

Athanasios Stubos, National Center for Scientific Research Demokritos, Environmental Research Laboratory, Greece

Abstract: A new class of microporous materials with extremely interesting properties for a great number of applications is that of Metal Organic Frameworks (MOFs). MOFS are among others extensively investigated for gas and liquid separations and they are typically constructed by connecting metal clusters with organic linkers resulting in crystalline porous materials. In this work we used a combination of stochastic and deterministic molecular simulation techniques, in order to examine structural features as well as thermodynamic and dynamic properties in a MOF/gas system. Focus was placed in well-established MOF structures such as IRMOF-1 and UiO-66(Zr) and gases of technological relevance such as N2, CH4 and CO2. For each MOF/gas system, Monte Carlo simulations were performed in order to calculate e.g. the specific surface area, the percentage of loading, gas adsorption isotherms at various temperatures, etc. In a further step, molecular dynamics simulations were performed in order to calculate the respective self-diffusion coefficient of the gas molecules and examine the mechanism of selective separation due to the atomistic interactions between the the gas molecules and the MOF framework. The simulations were carried out using the parameters of UFF forcefield for a rigid MOF structure, while specialized parameters were taken into account for each gas. All forcefields used in this work, were parameterized according to the Lennard-Jones potential equation. The obtained results show that the combination of Monte Carlo and MD techniques can offer a valuable tool for the study of microporous materials such as MOFs. References 1. Yang, Q., Wiersum, A. D., Jobic, H., Guillerm, V., Serre, C., Llewellyn, P. L., Maurin, G. (2011). The Journal of Physical Chemistry C, 115(28), 13768–13774. 2. Allen, M. P., Tildesley, D. J., Banavar, J. R. (1989). Physics Today, 42(3), 105. 3. Frenkel, D., Smit, B., Ratner, M. a. (1997). Physics Today (Vol. 50, pp. 1–628). 4. Gotzias, A., Heiberg-Andersen, H., Kainourgiakis, M., Steriotis, T. (2010). Applied Surface Science, 256(17), 5226–5231.
##### Passivity-based control of discrete-time SIS epidemics in complex networks

Luis Alarcon-Ramos, Universidad Autonoma Metropolitana - Cuajimalpa, Applied Mathematics, Mexico

Roberto Bernal-Jaquez, Universidad Autonoma Metropolitana - Cuajimalpa, Matematicas Aplicadas y Sistemas, Mexico

Alexander Schaum, Christian-Albrechts-Universität zu Kiel, Institute of Electrical and Information Engineering, , Germany

Carlos Rodriguez-Lucatero, Universidad Autonoma Metropolitana - Cuajimalpa, Tecnologias de la Informacion, Mexico

Abstract: In this paper, we investigate and analyze the complex dynamics of an epidemic model, in the light of passivity properties for feedback control design. For this purpose, a framework of discrete time Markov process dynamical systems is employed, to propose a control mechanism that allows to bring the system to the extinction state. Considering individual controls at single nodes with relative degree one, an approach is presented to identify the nodes to be controlled in the network, based on the property of feedback equivalence to a passive system. Therefore, taking the set of controlled nodes as degree of freedom, sufficient conditions are derived for an asymptotically stable associated zero dynamics. This implies that it is not necessary to control all nodes, but only a minimal set of nodes for which the associated zero dynamics are asymptotically stable, and the complete network dynamics become passiviable via state-feedback control. Simulation results are presented for large free-scale and regular networks, that corroborate our theoretical findings.
##### Monte Carlo studies of clathrate hydrates for applications that involve energy gases

Nikolaos Papadimitriou, National Center for Scientific Research Demokritos, Environmental Research Laboratory, Greece

Ioannis Tsimpanogiannis, National Center for Scientific Research Demokritos, Environmental Research Laboratory, Greece

Ioannis Economou, Texas A&M University at Qatar, Chemical Engineering Program, Qatar

Athanasios Stubos, National Center for Scientific Research Demokritos, Environmental Research Laboratory, Greece

Abstract: Clathrate hydrates are ice-like crystalline materials whose crystal lattice contains cavities that can trap gas molecules. Recently, they have attracted strong research interest as materials with potential use in applications that involve the storage, transportation, or separation of energy gases. The present work was motivated by two applications in this field: storage of H2, and CH4/CO2 separation. The concept of our work is to use Grand Canonical Monte Carlo (GCMC) simulations to calculate the number of gas molecules that can be trapped in every cavity of the hydrate (cavity occupancy). In this approach, the gas content of the hydrate at equilibrium is calculated as the result of an adsorption process where the rigid hydrate structure is the adsorbent. The simulations are performed over a wide range of conditions (temperature, pressure, gas mixture composition) and all of the most common hydrate structures (sI, sII, sH) are studied. All simulations were carried out with the MCCCS Towhee code. In the case of H2, the knowledge of cavity occupancies allows the estimation of the amount of gas that can be stored in the hydrate at specified conditions and the evaluation of the hydrogen-storage capacity of each hydrate structure [1]. For reasons of validation, several water models (e.g. SPC/E, TIP4P/Ice, TIP5P) are used along with various sets of interaction parameters for the H2 molecule. The effect of the lattice constant on the H2 content of the hydrate is also studied. H2 is one of the few gases that present the effect of multiple occupancy (i.e. more than one molecules in the same cavity) [2] and this effect is extensively studied. In the case of CH4/CO2 mixtures, the cavity occupancies for each of the gases determine the selectivity of the hydrate for the separation of this mixture. CH4 fills both the small and large cavities of the sI hydrate while CO2 is assumed to fill only the large ones. This different behavior between the two gases is an extra challenge from a computational point of view. To investigate the CH4/CO2 system, first we study the pure gases [3] and then we proceed to binary mixtures of various compositions. References: 1. N. I. Papadimitriou, I. N. Tsimpanogiannis, A. K. Stubos, Colloids. Surf. A 357:67-73 (2010). 2. N. I. Papadimitriou, I. N. Tsimpanogiannis, A. Th. Papapioannou, A. K. Stubos, J. Phys. Chem. C 112:10294–10302 (2008). 3. M. Lasich, A. H. Mohammadi, K. Bolton, J. Vrabec, D. Ramjugernath, Fluid Phase Equilib. 369:47-54 (2014).
##### Influence of Mobility Models in Precision Spray Aided by Wireless Sensor Networks

Leandro Bertini Lara Gonçalves, Universidade Estadual Paulista "Júlio de Mesquita Filho", , Brazil

Fausto Costa, ICMC-USP, , Brazil

Leandro Alves Neves, São Paulo State University (UNESP), Department of Computer Science and Statistics - DCCE, Brazil

Geraldo Francisco Donega Zafalon, Sao Paulo State University - UNESP/SJRP, Computer Science and Statistics - DCCE, Brazil

Jó Ueyama, ICMC-USP, , Brazil

Carlos Montez, Universidade Federal de Santa Catarina - UFSC, , Brazil

Alex Sandro Pinto, Universidade Federal de Santa Catarina, Câmpus de Blumenau, Brazil

Abstract: Applications of agricultural inputs using Unmanned Aerial Vehicles (UAVs) aided by Wireless Sensor Networks (WSN) is an approach related to the spray problem in precision agriculture. The efficiency of this approach depends of the correct communication between UAV and WSN. This communication might be affected by the method the UAV uses the network. Thus, in this work we perform comparisons among classical mobility models as Random Walk, Random Waypoint, Random Direction, Manhattan Grid, and we analyze how these models act in the the communication performance between UAV and WSN using the IEEE 802.15.4 protocol for random and grid networks. The tests were performed using the OMNeT++ simulator with MiXiM package and BonnMotion for the generation of mobility scenarios. The mobility models were set to be the same conditions for comparisons. In this study were collected and analized in each test: the amount of received frames by UAV; the amount of sent frames; the amount of backoffs from WSN; and the efficiency of network communication. The obtained results show that the network has similar performance both using random and grid implementation. The efficiency decreases as the density of nodes increases, mainly due to the concurrency by the transmission medium and saturation of the received message threshold of the UAV for densities above 3.67 X 10-4 nodes/m2. In the results it possible to check that Random Direction and Manhattan Grid had both similar performance when analyzed the communication efficiency of the mobility models and better results when compared with other models.
##### Frequency comb generation for wave transmission through the nonlinear dimer

Sadreev Almas, L.V. Kirensky Institute of Physics, Laboratory of theory of nonlinear processes, Russian Federation

Konstantin Pichugin, L.V. Kirensky Institute of Physics, , Russian Federation

Abstract: We study dynamical response of the nonlinear dimer relative to monochromatic wave injected via the waveguide. We show existence of a domain in space of frequency and injected amplitude where the stationary solutions of the temporal equations do not exist. We present time dependent solutions which show that scattering waves carry multiple harmonics with frequencies spaced equidistantly.
##### From single- to multi-layer silicene

Patrick Vogt, Technische Universität Berlin, Institut für Festkörperphysik, Germany

Abstract: Silicene, a novel silicon allotrope, which does not exist in nature, was theoretically conjectured a few years ago [1] as a stand-alone material. However, for the synthesis of silicene a substrate material is essential that prevents the formation of bulk silicon and, at the same time, does not react with the Si-atoms. Recently, silicene mono-layers could be synthesized on different substrates such as Ag(111), Ir(111) or ZrB2 [2]. The successful preparation of single layer silicene rises the question if the multi-layer silicene structure can also be synthesized which might be understood as the Si analogue to graphite. In this talk the epitaxial formation of single layer silicene on a Ag(111) substrate will be briefly summarized. Based on these results the formation of multilayer silicene will be discussed that forms upon further deposition of Si onto the first silicene layer [3]. Atomic structure, stacking order, electronic properties and growth aspects are discussed to demonstrate the differences of theses silicene multi-layer structure to bulk Si(111). The formation of multi-layer silicene could open new possibilities for the technological applicability of silicene. 1) G. G. Guzmán-Verri, L. C. Lew Yan Voon, Phys. Rev. B 76, 075131 (2007); S. Lebègue and O. Eriksson, Phys. Rev. B, 79 115409 (2009); S. Cahangirov et al., Phys. Rev. Lett. 102, 236804 (2009). 2) P. Vogt et al., Phys. Rev. Lett. 108, 55501 (2012); C.-L. Lin et al., Appl. Phys. Exp. 5, 045802 (2012); B. Feng et al., Nano Lett. 11, 3507 (2012); A. Fleurence et al., Phys. Rev. Lett. 108, 245501 (2012); L. Meng et al., Nano Lett. 13, 685 (2013). 3) P. Vogt et al., Appl. Phys. Lett. 104, 021602(2014), A. Resta et al. , Sci. Rep. 3, 2399 (2013), De Padova et al., Appl. Phys. Lett. 102, 163106 (2013).
##### Universal quasi-equilibrium state of self-gravitating systems not depending on dimensions

Tohru Tashiro, Ochanomizu University, Department of Physics, Japan

Abstract: Universal quasi-equilibrium state not depending on dimensions has been observed in the universe. Globular clusters which consist of about hundreds of thousands of stars are well-known to have a spherically symmetric number density. From a survey of the Herschel space observatory, it has been unveiled that there are about 30 filamentary structures of molecular clouds in IC 5146, and they have a cylindrically symmetric density profile in real space [1]. In both cases, a major factor constructing these structures is the gravity, which means that we can treat the systems as self-gravitating systems: Owing to the symmetries, filamentary structures of molecular clouds and globular clusters can be treated as the two and three dimensional self-gravitating system, respectively. The density profiles of these self-gravitating systems are constant around the center of the system and decrease in a power law at large radius. This is the universal quasi-equilibrium state. Until now, no one can explain the physics behind the universality. In this presentation, we will show that the same density profile as the observations are obtained by N-body simulations of self-gravitating system. Next, we will improve the dimension of our model, which can describe the density profile of three-dimensional self-gravitating system [2], and then we will exhibit that the universal distributions can be derived uniformly by the improved model. [1] D. Arzoumanian, et al. A&A 529 (2011) L6. [2] T. Tashiro and T. Tatekawa, J. Phys. Soc. Jpn. 79 (2010) 063001. Numerical Simulations of Physical and Engineering Processes, Jan Awrejcewicz (Ed.) (INTECH, Croatia) (2011) 301-318.
##### Shubnikov state in a Technetium film

JESUS GONZALEZ, UNIVERSIDAD DEL MAGDALENA, DEPARTAMENTO DE FISICA, Colombia

Jose Jose Barba Ortega, Universidad Nacional de Colombia, Physics, Colombia

Abstract: We numerically investigate the flux quantum configurations and some thermodynamical properties of a Technetium film by using the link variables technique for one shape of circular geometry. The technetium exhibit superconductivity properties indicated by to an extrapolated critical magnetic field value of Hc (T=0)=1410 Oe , and a Ginzburg-Landau parameter of k=0.92, being the magnetic behavior of this material characteristic of a type-II superconductor with a weak-coupling superconductor of the BCS type. The studied sample is a circular sector with angular width surrounded by a dielectric material and submitted to external magnetic field applied perpendicular to its plane. We evaluate the magnetic moment density, Shubnikov state and thermodynamical fields as a function of the external magnetic field.
##### New trace theorems for neutronic function spaces

Boulanouar, LMCM-RSA, , France

Abstract: Let $X\subset\R^n$ ($n\ge1$) be an open subset with $C^1$ boundary $\partial X$ and let $d\mu$ be a Radon measure on $\R^n$ ($n\ge1$) whose support is $V$. Let $\Gamma_\mp=\left\{(x,v)\in\partial X\times V \; : \; \mp v\cdot n(x)>0\right\}$ where, $n(x)$ is the outer unit normal at $x\in\partial X$. For each $(x,v)\in X\times V$, we set $t(x,v)=\inf\{t,\;\; x-tv\not \in X\}$ and $\theta(x,v)=t(x,v)+t(x,-v).$ Similarly, if $(x,v)\in \Gamma_+$ we set $\tau(x,v)=\inf\{t,\;\; x-tv\not \in X\}.$ Let $$\widetilde{W}^p(X\times V) =\aco{ \varphi\in L^p(X\times V),\;\; v\cdot\nabla_x\varphi\in L^p(X\times V) }$$ whose norm is $\norm{\varphi}_{\widetilde{W}^p(X\times V)}= \left[\norm{\varphi}_p^p+\norm{v\cdot\nabla_x\varphi}_p^p \right]^{\frac{1}{p}}$ and let $L^p(\Gamma_\pm, d\widetilde{\xi}_\pm)$ be the weighted trace spaces whose norm is \begin{equation*}\label{} \norm{\varphi}_{L^p(\Gamma_\pm, d\widetilde{\xi}_\pm)}= \left[\int_{\Gamma_{\pm}}\abs{\varphi(x,v)}^pd\widetilde{\xi}_\pm\right]^{\frac{1}{p}} \quad \text{where,}\quad d\widetilde{\xi}_\pm=\min\{\tau(x,\pm v), 1\} \abs{v\cdot n(x)}d\gamma d\mu(v). \end{equation*} Finally, let $\gamma_\mp\, :\; \varphi \longrightarrow \varphi\mid_{\Gamma_\mp}$ be the trace mappings whose continuity is given by \begin{Lem}[\cite{Cessenat1}\cite{Cessenat2}] \label{CES} Both trace mappings $\gamma_\mp : \widetilde{W}^p(X\times V) \longrightarrow L^p(\Gamma\!_\mp, d\widetilde{\xi}_\mp)$ are continuous. \end{Lem} \noindent It is well known that Lemma~\ref{CES} is not suitable in practice. Indeed, let $\varphi\in \widetilde{W}^1(X\times V)$ such that $\gamma_+\varphi$ or $\gamma_-\varphi$ exists. Green's formula yields that $$\int_{X\times V} v\cdot\nabla_x\abs{\varphi}dxd\mu(v) =\int_{\Gamma_+}\abs{\gamma_+\varphi(x,v)}d\xi -\int_{\Gamma_-}\abs{\gamma_-\varphi(x,v)}d\xi.$$ One obviously sees that $L^1(\Gamma_\pm, d\xi)$ are natural trace spaces and therefore Lemma~\ref{CES} can not be used because of $d\xi\not=d\widetilde{\xi}$. For instance, let $X\times V=(0,1)\times(0,\infty)$ for which $\Gamma_\pm=(0,\infty)$ and let $\varphi(x,v)=\frac{1}{v^2+1}$ be in $\widetilde{W}^p(X\times V)$. We obviously have $\gamma_\pm\varphi=\varphi\in L^p(\Gamma_\pm, d\widetilde{\xi}_\pm)$ and $\gamma_\pm\varphi\not\in L^p(\Gamma_\pm)$. \par\noindent It is therefore natural to raise the following question : \emph{Does there exist a suitable space $W^p(X\times V)$ $(p \ge 1)$ on which the trace mappings $\gamma_\pm$ are continuous into their natural trace spaces $L^p(\Gamma_\pm)$ ?} A positive answer is given by \begin{The}[\cite{Boulanouar}]\label{THE} Let $W^p(X\times V)$ $(p \ge 1)$ be the following Banach space \begin{equation*} W^p(X\times V) = \aco{\varphi\in L^p(X\times V), v\cdot\nabla_x\varphi\in L^p(X\times V), \theta^{-\frac{1}{p}}\varphi\in L^p(X\times V)} \end{equation*} whose norm is $$\norm{\varphi}_{W^p(X\times V)}= \left[(p-1)\norm{\varphi}_p^p+\norm{v\cdot\nabla_x\varphi}_p^p +\norm{\theta^{-\frac{1}{p}}\varphi}_p^p\right]^{\frac{1}{p}}.$$ Then, the trace mappings $\gamma_+\, :\; W^p(X\times V) \longrightarrow L^p(\Gamma_+)$ and $\gamma_-\, :\; W^p(X\times V) \longrightarrow L^p(\Gamma_-)$ are continuous. Furthermore $\norm{\gamma_+\varphi}_{L^p(\Gamma_+)} \le \norm{\varphi}_{W^p(X\times V)}$ and $\norm{\gamma_-\varphi}_{L^p(\Gamma_-)} \le\norm{\varphi}_{W^p(X\times V)}.$ \end{The} \begin{thebibliography}{99} \bibitem{Boulanouar} \textsc{M. Boulanouar.} \textsl{New trace theorem for neutronic function spaces.} Trans. Theor. Stat. Phys., 38, 228--242, 2009. \bibitem{Cessenat1} M. Cessenat. \textsl{Th\'eor\emes de trace $L^p$ pour des espaces de fonctions de la neutronique}, C. R. Acad. Sc. Paris, t.299, pp. 831--834, 1984. \bibitem{Cessenat2} M. Cessenat. \textsl{Th\'eor\emes de trace pour des espaces de fonctions de la neutronique}, C. R. Acad. Sc. Paris. t.300, pp. 89--92, 1985. \end{thebibliography}
##### Generation of thrust and lift with airfoils in plunging and pitching motion

Manuel Moriche, Universidad Carlos III de Madrid, Bioengineering and Aerospace Department, Spain

Abstract: The aerodynamics of flapping wings has been the focus of attention for many researchers for the last decades. Classical unsteady models of Wagner and Theodorsen have been widely used, but they become inadequate at high amplitude/frequency motions and very low Reynolds number, where the flow is massively separated. This is precisely the regime where insects and small birds fly. Improving our understanding of this regime is a key factor to unlock the standardized design of bio-inspired Micro Air Vehicles (MAV), able to achieve high maneuverability and efficiency in their flapping flight. In this talk we present fully resolved Direct Numerical Simulations (DNS) of flow over moving airfoils (2D). The solver used is an in-house code that implements the Immersed Boundary method of Uhlmann1 . A combination of sinusoidal plunging and pitching motions is imposed in rigid airfoils. Starting from a thrust producing case (Reduced frequency, k = 1.41. Plunging amplitude h/c = 1. Mean pitching angle theta_m = 0º. Pitching amplitude theta_0 = 30◦ . Phase shift phi = 90◦ ), we increase the mean pitching angle (in order to produce lift) and vary the phase shift between pitching and plunging (to optimize the direction and magnitude of the net force on the airfoil). These cases will be discussed in terms of their lift coefficient, thrust coefficient and power efficiency.
##### Bistability in self-activating genes regulated by non-coding RNAs

Jesus Miro-Bueno, Research Institute of the IT4Innovations Centre of Excellence, Faculty of Philosophy and Science, Silesian University in Opava, , Czech Republic

Abstract: Non-coding RNA molecules are able to regulate gene expression and play an essential role in cells. On the other hand, bistability is an important behaviour of genetic networks. Here, we propose and study an ODE model in order to show how non-coding RNA can produce bistability in a simple way. The model comprises a single gene with positive feedback that is repressed by non-coding RNA molecules. We show how the values of all the reaction rates involved in the model are able to control the transitions between the high and low states. This new model can be interesting to clarify the role of non-coding RNA molecules in genetic networks. As well, these results can be interesting in synthetic biology for developing new genetic memories and biomolecular devices based on non-coding RNAs.
##### An efficient implicit direct forcing immersed boundary method for incompressible flows

shanggui, Univiersité de Technologie de Compiègne, , France

Abstract: We present a novel efficient implicit direct forcing immersed boundary method for incompressible flows with complex boundaries, based on the work of [M. Uhlmann, An immersed boundary method with direct forcing for the simulation of particulate flows, J. Comput. Phys., 209(2005) 448-476]. In Uhlmann's work, the calculation is performed on the Eulerian grid regardless the immersed object, with a fictitious force to mimic the presence of the physical boundaries. It is also known that Uhlmann's explicit method fails to accurately impose the non-slip boundary condition on the immersed interface. To improve the accuracy, iteration of the whole system is required, which could be extremely time-consuming. In the present work, an implicit treatment of the artificial force is preferred, but with an effective way of system iteration. The Navier-Stokes equation is resolved with the rotational incremental pressure-correction projection method of [J.L. Guermond and J. Shen, An overview of projection methods for incompressible flows, Comput. Methods Appl. Mech. Engrg., 195(2006) 6011-6045], thus a higher precision is reserved. Several numerical simulations are performed to verify the proposed method, which are in good agreement with those in the literature.
##### Evaluation of connectedness between the University courses of Physics and Chemistry basing on the graph model of inter-subject links

Tatiana Gnitetskaya, Far Eastern Federal University, School of natural sciences, Russian Federation

Elena Ivanova, Far Eastern Federal University, , Russian Federation

Abstract: An application of the graph model of inter-subject links to University courses of Physics and Chemistry is presented in this article. A part of inter-subject space with directions of inter-subject links from Physics to Chemistry in the group of physical concepts has been shown. The graph model of inter-subject links includes quantitative indicators. Its numerical values are given in the article. The degree of connectedness between the data of Physics and Chemistry courses is discussed for the courses considered. The effect of the courses placement within a curriculum on the value of their connectedness is shown. The placement of courses within a curriculum can provide the study of the courses at the same time or consecutive study, when one course precedes another.
##### SOLITONIC EFFECTS IN SUPERCONTINNUM GENERATION IN HIGHLY NONLINEAR FIBERS

Mário Ferreira, University of Aveiro, Physics, Portugal

Abstract: When ultrashort pulses with large enough power are launched into highly nonlinear fibers, soliton fission gives origin to multiple fundamental solitons of different widths and peak powers. Among the host of soliton related effects contributing to supercontinuum generation, there are two which become particularly important: the soliton self-frequency shift and the emission of dispersive radiation in the normal dispersion region. The peculiar dispersive characteristics of highly nonlinear fibers play a key role in these circumstances. If long duration pulses or a continuous wave are used to pump the supercontinuum, the modulation instability becomes the main initial mechanism, breaking the pump wave into multiple fundamental solitons. In this case, both the soliton self-frequency shift and inelastic collisions among such solitons play an important role. In this paper we discuss in detail the solitonic effects in the supercontinuum generation process for different pumping conditions.
##### COMPTON OPERATOR IN QUANTUM ELECTRODYNAMICS

Abstract: In the frame in the quantum electrodynamics exist four basic operators; the electron self-energy , vacuum polarization , vertex correction , and the Compton operator, denoted by C. The first three operators are very important by its relation with renormalized and Ward identity. However, the Compton operator has equal importance, but without divergence, and little attention has been given it. We have calculated the Compton operator and obtained the closed expression for it in the frame of dimensionally continuous integration and hypergeometric functions.
##### Representing Structure formation in the Early Universe as a result of non Linear Electrodynamics influencing scale factor size

Andrew Beckwith, Chongqing University department of physics( visitor), , China

Abstract: We find that having the scale factor close to zero due to a given magnetic field value in, an early universe magnetic field affects how we would interpret Mukhanov’s chapter on ‘self reproduction of the universe’ in in his reference “Physical foundations of cosmology” terms of production of inhomogeneity during inflation and its aftermath. The stronger an early universe magnetic field is, the greater the likelihood of production of about 20 new domains of size 1/ H, with H early universe Hubble’s constant, per Planck time interval in evolution. One final caveat to consider. What may happen is that the Camara (2004) density and Quintessential density (Corda et al.) are both simultaneously satisfied, which would put additional restrictions on the magnetic field which in turn affects structure formation. In time, once Eq.(16) of this paper is refined further, the author hopes that some of the issues raised by Kobayashi and Seto as to allowed inflation models may be addressed, once further refinement of these preliminary results commences
##### Numerical Investigation of Grain Coarsening and Coalescence Model

Aliki Muradova, Tecnical University of Crete, Mineral Resources Engineering, Greece

Dionisios Hristopulos, Technical University of Crete, Department of Mineral Resources Engineering, Greece

Abstract: A kinetic nonlinear model of mass transfer, grain coarsening and coales- cence with potential applications in sintering processes is considered. The model involves nonlinear differential equations that determine the transport of mass between grains. The rate of mass transfer is controlled by the acti- vation energy which depends on the degree of amorphization of the grains, leading to a nonlinear model of mass transfer and grain coarsening. The resulting dynamical system with random initial conditions (grain mass con- guration) is solved by means of the Runge-Kutta method. In order to incorporate coalescence of smaller grains with larger neighbors a computa- tional algorithm is adapted. An analysis of the two-grain system is carried out and the solution of the multi-grain system is studied. Numerical results are presented by examples.
##### Basic idea of Corbino-type single-electron transistor

Akira Suzuki, Tokyo University of Science, Physics, Japan

Abstract: We have formulated the transmission probability of an electron in a Corbino quantum disk by taking into account charging effect. The geometrical potential of the Corbino disk has a singularity at the centre of the disk. In order to avoid this singularity problem, we have to reformulate the Schroedinger equation in the Riemannian manifold. The Schroedinger equation describing the motion of the electron in the Corbino disk must be expressed by introducing a momentum operator reformed by the metric tensor. In order to obtain a Hermitian momentum operator, we must deform the Hilbert space by introducing a new wave function. This deformation leads to the extra potential term in the Schroedinger equation, which depends on the metric, i.e., the geometry of the disk. It should be noted that the charging energy of confining electrons in the Corbino disk should depend on the geometry of the disk. We discuss the quantum tunneling of an electron confined in the Corbino disk in order to investigate the effect of both geometrical potential and charging energy of confining electrons in the Corbino disk by using the Wentzel-Kramers-Brillouin (WKB) method. It is expected that the charging energy, which depends on the effective confining potential, plays an important role in the transmission probability. This suggests that the formulated transmission probability is applicable to the analysis of the single-electron transistor.
##### Mathematical modeling of a rigid circular membrane in a tri-material transversely isotropic full-space undergoing a prescribed displacement

Ladan Mohtati, School of Civil Engineering, College of Engineering, University of Tehran, , Iran (Islamic Republic of)

Mohammad Rahimian, School of Civil Engineering, College of Engineering, University of Tehran, , Iran, Islamic Republic Of

Morteza Eskandari-Ghadi, School of Civil Engineering, College of Engineering, University of Tehran, , Iran, Islamic Republic Of

Abstract: A horizontal layer between two different transversely isotropic half-spaces forms a tri-material full-space, which constitutes the domain for the boundary value problem involved in this paper. The axes of symmetries of all materials are assumed to be parallel to each other and normal to the interface of the middle layer and surrounding half-spaces. A mathematical formulation is presented to determine the response of a rigid circular membrane, which is laid down at an interface of the tri-material transversely isotropic full-space and is considered to be under a prescribed translation, which happens parallel to the plane of membrane. The governing equations are expressed in the cylindrical coordinate system, where the axis of the symmetry of the cylinder and the material are aligned. With the aid of a system of two scalar potential functions, the governing equations of motion can be uncoupled into two separated partial differential equations, which may be transformed to some ordinary differential equations by applying the Hankel integral transforms in the radial direction and Fourier series along the angular coordinate. After determining the unknown functions by imposing the relaxed boundary conditions, they are transformed to a set of four coupled integral equations, which are reduced to two coupled Fredholm-Volterra integral equations of the second kind, from which both displacement and the stress fields are computed. The proposed solutions are applied for a transversely isotropic half-space and the results match with the existing solutions. To confirm the accuracy of the numerical evaluation of the integrals involved, the numerical results are compared with the solutions exists for a transversely isotropic half-space. In order to investigate the degree of material anisotropy, some numerical evaluations are given for different combinations of transversely isotropic region. Final results show that the boundary conditions are completely satisfied and the verified numerical method leads to compatible solutions.
##### Quantum Toeplitz Matrix Appearing in Numerical Approximations of Quantum Integrals

S.M.Hashemiparast, K.N.Toosi University of Technology, Mathematics and Statistics, Iran (Islamic Republic of)

Abstract: Many complex systems in various area of science and technology can be described by supper statistics ,in model building situations the usual procedure is to select a model from a parametric family of distributions or generally the pathway models have been proposed which have connections with quantum integral and quantum fractional integral ,finally numerical approximations of quantum integral. In this paper we present a Quantum Toeplitz matrix for which the numerical q-calculation of the q- integrals based on specified nodes and weights are carried on with minimized error, in this q-calculation we consider the various properties of this Toeplitz matrix and the relations with the ordinary Toeplitz matrix, first the conditions for invetibility of this quantum Toeplitz matrix can be determined, then for having a unique solution, the conditions for the original system of equations (q=1 or ordinary Toeplitz matrix) is compared for the different values of q , so, some numerical examples are given for comparison with the exact solution and alternative numerical solution based on procedures using quantum Gauss -quantum Hermit quadrature rules, and also based on q-variation with iteration method , the related tables of approximations for all methods are presented .
##### Graph Model of Intradisciplinary Connections in Example of General Physics Course

Tatiana Gnitetskaya, Far Eastern Federal University, School of natural sciences, Russian Federation

Abstract: The model of an intradisciplinary connections was elaborated on the base of the theory of graphs. Every connection which appears in training content may be presented as oriented marked graphs. Each graph is a tree. In this paper we presented definition and model of intradisciplinary connections for example of physics course. The quantitative parameters of model are described in this paper. Quantitative method based on this model could be help to optimize a content of physics course. Furthermore using this model we can distinguish fundamental notions, laws and other elements of knowledge to separate group. This group is very important during the process of creating training course. Method of semantic structure attached to content of physics’ lections and physical problems was described.
##### Characterizing the non-stationary evolution of volume-prices evolution in the New York Stock Market

Paulo Rocha, Centro de Matemática e Aplicações Fundamentais, , Portugal

Frank Raischel, Instituto Dom Luiz, CGUL, University of Lisbon, 1749-016 Lisbon, Portugal, , Portugal

Pedro G. Lind, ForWind and Institute of Physics, University of Oldenburg, DE-26111 Oldenburg, Germany, , Portugal

Abstract: We introduce a framework for describing the stochastic evolution of the parameters defining volume-price distributions of the New York stock market. The volume-price data appears to follow a specific statistical pattern, other than the evolution of prices measured in similar studies. We argue that the inverse Gamma distribution fits well the volume-price evolution, with the two corresponding parameters varying in time. Assuming that the evolution of these parameters is governed by coupled Langevin equations, we derive the corresponding drift and diffusion coefficients, which then provide insight for understanding the mechanisms underlying the evolution of the stock market. The data analysed comprises volume-prices distributions extracted from yahoo available data with a sampling frequency of 10 minutes.
##### Multiple Time Series Ising Model for Financial Market Simulations

Tetsuya Takaishi, Hiroshima University of Economics, , Japan

Abstract: In this paper we propose an Ising model which simulates multiple financial time series. Our model introduces a global spin which couples to spins of other systems. Depending on the value of the global spin, correlations between time series will appear. Simulations from our model show that time series exhibit the volatility clustering that is often observed in real financial markets. We also find non-zero covariance of volatility time series. Thus our model can simulate stock markets where volatilities of stocks are mutually correlated.
##### Effects of Disorder on Critical Behavior of Generalized Conserved Lattice Gas

Meesoon Ha, Chosun University, Department of Physics Education, Korea, Republic of

##### Approximate analytical channel potential model of poly-Si thin film transistors operated in the strong inversion region under the high gate and low drain biases

Zhen Zhu, Suzhou Vocational University, , China

Junhao Chu, East China Normal University, , China

Abstract: An approximate analytical channel potential model of polycrystalline silicon thin film transistors operated in the strong inversion region under the high gate and low drain biases is proposed. Thus, the linear relationship between the channel potential and the drain voltage is derived in the strong inversion region under the above bias condition when the polysilicon layer is ultrathin. This model agrees with the two-dimensional-device simulation results under different gate voltages, different drain voltages and different channel lengths. By comparing the relative errors between the model and the simulation results, it presents that this model is more suitable under the higher gate voltage Vg or the lower drain voltage Vd, regardless of the channel length. And this approximate analytical model is helpful in solving the two-dimensional-device problem by one-dimensional Poisson’s equation since the drain bias is taken into account in the channel potential.
##### First-principle studies of phonons and thermal properties of AlN in wurtzite Structure

J.Q. Fu, Inner Mongolia University, Physical Science and Technology, China

Tie Lei Song, Department of Physics, School of Physical Science and Technology, Inner Mongolia University, , China

Xi Xia Liang, Inner Mongolia University, Department of Physics, School of Physical Science and Technology, , China

Guo Jun Zhao, Inner Mongolia University, Department of Physics, China

Abstract: Group-III nitrides have attracted considerable attention during the past decade due to technological applications in the optoelectronic and electronic devices. Aluminum nitride is one of the important materials among them. It has high thermal conductivity, high melting point, large bulk modulus, and large band gap. In this work, we calculate the band structure, density of states, phonon dispersion and thermodynamic properties of A1N in wurtzite structure based on the first-principle with the software QUANTUMESPRESSO. The local density approximation (LDA) and the generalized gradient approximation (GGA) exchange- correlation potentials are applied in the calculations. The results show that A1N belongs to direct band gap semiconductor. The phonon dispersion characteristics and thermodynamic properties are discussed in detail. The obtained results for the entropy S and the constant-volume specific heat CV as functions of the temperature T based on GGA and LDA are given. CV and S increase with the increasing T. Both the values of S and CV by LDA are slightly lower than those by GGA. The calculated values are in agreement with available experimental data.
##### On Modelling Quantum Potential Flow

KUO CHUNG HSUAN, Department of Aeronautics and astronautics, NCKU, , Afghanistan

Abstract: In this paper, a new hydrodynamic formulation of complex-valued quantum mechanics is derived to reveal a novel analogy between the probability flow and the potential flow on the complex plane. For a given complex-valued wavefunction , , we first define a complex potential function with and then prove that the streamline lines and the potential lines in the potential flow defined by are equivalent to the constant-probability lines and the constant-phase lines in the probability flow defined by . The discovered analogy is very useful in visualizing the unobservable probability flow on the complex plane by analogy with the 2D potential flow on the real plane, which can be visualized by using dye streaks in a fluid laboratory.
##### THE PROBLEM BY CHOICE OF PRECEDENT WITH EXAMPLE OF DATA ENCRYPTION ALGORITHMS

Ayman Iskakova, L.N. Gumilyov Eurasian National University, Department of Mechanical and Mathematics, Kazakhstan

Abstract: In this work the process of operational withdrawal of knowledge matrix on precedent described situational vectors realization of computing technique is presented. The new method of withdrawal precedent and its application in data encryption algorithms are presented.
##### Risk management mechanism of payment transactions

Ayman Iskakova, L.N. Gumilyov Eurasian National University, Department of Mechanical and Mathematics, Kazakhstan

Abstract: The mechanism of protection of smart card software from unauthorized copying based authentication based on symmetric has presented
##### Group theory Analysis of the Magnetic Structures in Rare Earth Iron Garnets revisited

Mahieddine LAHOUBI, Badji Mokhtar-Annaba University, Faculty of Sciences , Physics, Laboratory, L.P.S., Algeria

Abstract: The rare earth iron garnets with general formula REIG (where RE is a trivalent rare earth ion or the yttrium) discovered by Bertaut and Forrat [1] have been extensively studied during the last decades for the large application as well as fundamental works. Nevertheless, it appears up today that the magnetic and dielectric properties of these ferrimagnets are not achieved as indicated by the large magnetodielectric (MD) effects in low external magnetic fields and at low temperatures revealed by Hur et al. [2] and Song et al. [3] respectively in TbIG and DyIG. The magnetic ions are distributed over the three crystallographic sites of the paramagnetic space group (G) Ia-3d No. 230, RE ions in dodecahedral {24c} sites, and the iron Fe in both octahedral [16a] and tetrahedral (24d) sites respectively. According to the Néel theory of ferrimagnetism [4], the magnetic moments of the RE ions form below the Néel temperature (TN) (which is the same in all REIG (≈ 560 K) [5]) a collinear arrangement antiparallel to the resultant Fe magnetization along the [111] crystallographic direction which is the easy axis of magnetization in a given domain without external magnetic field. At low temperatures, due to the competition between three types of anisotropies (spin-orbit coupling, crystalline field and RE–Fe superexchange interactions) complicated noncollinear arrangement appears with the rhombohedral distortion from the space group cubic G to its highest subgroup (G’) R-3c No. 167. The ''representation analysis'' devised and developed by Bertaut for the determination of magnetic structures [6, 7] and the generalized identification method of Olbrychski [8, 9] of the irreducible representations (irreps) of G and G’ are applied here. These elegant group theoretical techniques are used in the interpretation of our recent neutron powder diffraction studies of the temperature evolution of the ''double umbrella'' magnetic structure observed in TbIG and DyIG [10-12] and the occurrence of the MD effects in these compounds is discussed. According to the concept of the single irrep developed by Izyumov [13] and to the fact that the third power of the three-dimensional irrep of G, (Г4g = T1g)3 contains always the identity representation [14], one can conclude that the magnetic anomaly at the so called low-temperature point of Belov (TB = 54 K) cannot be a second order phase transition, despite the low-T symmetry is a subgroup of the high-T symmetry. [1] E. F. Bertaut, and F. Forrat, C. R. Acad. Sc. Paris, vol. 242 (1956) 382. [2] N. Hur, S. Park, S. Guha, A. Borissov, V. Kiryukhin, and S.-W. Cheong, Appl. Phys. Lett., vol. 87, (2005) 042901. [3] K. M. Song, Y. A. Park, K. D. Lee, B. K. Yun, M. H. Jung, J. Cho, J. H. Jung, and N. Hur, Phys. Rev. B 83 (2011) 012404. [4] L. Néel, Ann. Phys., t. 3 (1948) pp. 137-198. [5] R. Pauthenet, Thesis Grenoble, France, n°. Ordre 81, pp.1-38, 1958. [6] E. F. Bertaut, Acta Cryst. A. 24 (1968) 217. [7] E. F. Bertaut, J. Phys. Colloque C1, suppl. n°. 2-3, vol. 32 (1971) C1- 462. [8] E. F. Bertaut, C. R. Acad. Sc. Paris, vol. 268 (1969) 281. [9] K. Olbrychski, Phys. Stat. Sol., 3, (1963) 2143. [10] M. Lahoubi, “Temperature evolution of the double umbrella magnetic structure in terbium iron garnet”, in Neutron Diffraction, (Ed) I. Khidirov, (2012) Chap. 10, pp. 203-230. [11] M. Lahoubi, Journal of Phys. Conf. Ser., 340 (2012) 012068. [12] M. Lahoubi, W. Younsi, M.-L. Soltani, and B. Ouladdiaf, Journal of Phys. Conf. Ser., 200 (2010) 082018. [13] Yu. A. Izyumov, Physics–Uspekhi, vol. 23, 356−274, 1980 [14] L. D. Landau, and E. M. Lifshitz, Statistical physics, Pergamon Press, 1969.
##### An Axisymmetric Contact Problem for an Elastic Layer Subjected to a Tensile Stress Applied Over an Annular Region

Kebli Belkacem, Ecole Nationale Polytechnique, Mechanical Engineering, Algeria

Abstract: The study is concerned with an axisymmetric contact problem for an elastic layer on which a tensile uniform stress is applied on the bottom over an annular hole while the lower surface is rigidly clamped. By using the Hankel integral transforms method we reduce the three-part mixed boundary value problem to a system of triple integral equations. With the help of the Gegenbauer formula and some integral representations of the Bessel function, we get an infinite system of algebraic equations for determining the unknown function. The expressions of the stress intensity factors are given analytically. Some quantities of physical interest are shown graphically followed by a discussion of the effect of the radii of the hole as well as the medium thickness on the layer deformation.
##### Precise Approximated Solution for the Bohm Sheath Potential

Pablo Martin, Antofagasta University, Physics, Chile

Luis Antonio Cortes Vega, Antofagasta University, Mathematics , Chile

Abstract: The Poisson equation for the plasma sheath potential near a wall, leads to a non-lineal differential equation, whose analytic solution is not know. The usual approximation taking only the first term does not give good accuracy. Other approximations taken two additional terms gives better accuracy , but it fails to give high accuracy in the intermediate region . Here we present an analytic approximations using quasi-rational multipoint approximation method (MPQA), which leads to new analytic approximate solution with much higher accuracy, and very precise results not only near wall and far away, but also in the transition region, Absolute values of the maximum per cent error as a function of the wall potential will be present. Other several figures showing thess new analytic solutions as a function of the relevant parameters for this problem. The advantages of the present solution compared with those of pervious works will be shown.
##### Reduction Solutions Related to the Residual Symmetry for a (2+1)-Dimensional Burgers Equation

Jianping FANG, Lishui University, , China

Abstract: The residual symmetry, coming from the standard truncated Painlev\'{e} expansion of a (2+1)-dimensional Burgers equation, is localized in the properly prolonged system with the Lie point symmetry vector. Some different transformation invariance are derived through the obtained symmetries. Further, reduction solution, especially interactive solution is obtained through a generalized tanh function expansion approach.
##### Effective mass Schrödinger equation with Thomas-Fermi potential

Gerardo Ovando, Universidad Autonoma Metropolitana, DCBI - Ciencias Básicas , Mexico

Jose Juan Peña Gil, Universidad Autonoma Metropoitana Azc., Ciencias Básicas, Mexico

Jesús Morales, Universidad Autónoma Metropolitana, DCBI- Ciencias Básicas, Mexico

Abstract: The exactly-solvable position-dependent mass Schrödinger equation (PDMSE) for the Thomas-Fermi potential is presented. This model has recent interest because its applications to δ-doped semiconductor structures. To that, the PDMSE is transformed into a standard Schrödinger equation with constant mass (CMSE) with the aim to obtain the transformation that should be used to find the exactly solvable CMSE. In that case, the potential associated to the PDMSE and the potential involved in the CMSE are connected through a Riccati-type relationship which includes a superpotential that determines the position-dependent mass distribution m(x) leading to exactly-solvable PDMSE. To attain the purpose of this work, we assume that the superpotential is constant to find m(x) and consequently an exactly-solvable PDMSE for the Thomas-Fermi potential. Orthogonal eigenfunctions and normalization constants are determined in closed form by developing the particular details involved in the Sturm-Liouville theory for the corresponding CMSE. It is proved that the eigenstates are those of definite parity with eigenvalues expressed through the order of the Bessel functions of the first kind. Beyond the case considered in this work, the approach is general and can be useful in the study of electronic properties of non-uniform materials in which the carrier effective mass depends on the position as well as in the search of new solvable potentials suitable for quantum systems.
##### A Simulation-Based Evaluation of Local Fluxes of Secondary Cosmic Rays at Tehran

Zahra bagheri, Research Institute for Astronomy and Astrophysics of Maragha, , Iran (Islamic Republic of)

Pantea Davoudifar, Research Institute for Astronomy and Astrophysics of Maragha, Astroparticle Physics, Iran (Islamic Republic of)

Abstract: Carrying segments, instruments and small detectors up to the upper atmosphere is a difficult and expensive task. So to evaluate the errors (hard and soft) caused by cosmic rays on onboard parts of a satellite, ordinarily ground-based experiments were being designed to evaluate these errors. In this paper taking the advantage of simulation, the local fluxes of secondary cosmic rays were calculated for Tehran (51°, 35 °). OMERE software were used to get the appropriate fluxes of cosmic rays at the first interaction level of the earth atmosphere with the height of 112.8 kilometers. Considering a flat surface on interaction level, appropriate area and energies were calculated and used to interpret the results of CORSIKA simulations. The local fluxes of secondaries of different types were calculated for Tehran.
##### Scattering of vortex-antivortex pair by single quantum vortex in a Bose-Einstein Condensate.

Lev A. Smirnov, Institute of Applied Physics of the Russian Academy of Sciences (IAP RAS), Nonlinear electrodynamics department, Russian Federation

Aleksander I. Smirnov, Institute of Applied Physics of the Russian Academy of Sciences (IAP RAS), , Russian Federation

Abstract: Dynamics of vortex structures and their interaction with each other in many respects determine key aspects of evolution for the ultracold Bose gas cloud with the repulsing interaction between atoms. Quantum vortices are associated with the breaking of the superfluidity mode and the transition of the Bose-Einstein condensate (BEC) to the turbulent state. Thus, it is important to make a maximal progress in solving the problem of describing different kinds of motion for interacting vortex formations. In the case when the distance between a vortex and an antivortex is substantially smaller than the characteristic scale of the medium inhomogenity the distributions of both the density and the velocities field in vortex pairs are similar to those that take place for the homogeneous flow-free condensate in the correspoding solitary solutions of the Gross-Pitayevsky (GP) equation. In this case one may introduce a concept of the “two-dimensional dark quasisolitons”. Such quasisolitons represent the dips of BEC concentration, propagating in the inhomogeneous Bose gas at subsonic velocities with acceleration and, generally speaking, along the curved paths. We have developed the asymptotic theory describing behaviour of two-dimensional quasisolitons structures in a smoothly inhomogeneous condensate. Using this theory in the case of the quasisolitons motion in the smoothly inhomogeneous small flows of BEC we have succeeded in detailed analizing and explaining peculiarities of scattering for vortex-antivortex pairs on a single vortex. The results of the direct numerical simulation performed within the framework of the GP equation have demonstrated a good agreement with the developed theory.
##### Analytic approximations to the energy eingenvalues of the quadratic Zeeman effect in 2-D for hydrogenlike atoms

Pablo Martin, Antofagasta University, Physics, Chile

Flavio Citti, Simon Bolivar University, Physics, Venezuela

Enrique Castro, Simon Bolivar University, Physics, Venezuela

Luis Antonio Cortes Vega, Antofagasta University, Mathematics , Chile

Abstract: Analytic energy eigenvalues approximations have been found for the 2-D quadratic Zeeman effect in hydrogenlike atoms. Here a new technique has been used, which is an extension of the two-point quasirational approximation (TPQA) technique, previously used for the same problem . In our previous calculations, two expansions were used, one for small values of the magnetic field and other for large values, that is, expansions around zero and the inﬁnity were determined. Later a bridge between both expansions was built using simple auxiliary functions as well as rational ones. In the present paper expansions around intermediate points were also determined, and used. In this way, the corresponding analytic bridge function is determined, taken care also of these new expansions, so the technique is a multiple point quasirational approximation (MPQA). The calculations are carried out for the ground state 1s and the excited states 2p- and 3d- . These results are now much better than those obtained in previous works . We will also show the advantage of our approximants with respect to other approximations.
##### Bound state solutions of D-dimensional Feynman propagator for the q-deformed Woods–Saxon potential

Ahmed Diaf, Laboratory of Energy and Smart Systems, Khemis Miliana University, Sciences and Technology, Algeria

Abstract: Using the space-time transformations, approximate analytical solutions of the D-dimensional propagator in the presence of q-deformed Woods–Saxon potential are obtained. The analytical expression of the energy eigenvalues is given for various quantum numbers and the corresponding normalized eigenfunctions are obtained in terms of hypergeometric function. Our results are compared with those given by The Nikivorov-Uvarov method.
##### Designing a Supply Chain Network on the Basis of Variational inequality

Abstract: Increasing developments and changes in the field of production and business resulted in rising and development of concepts including supply chain and supply chain management. Supply change management from 190s onward developed increasingly among scientific societies, companies and industries and aimed at establishing cooperation among suppliers, lowering the costs, satisfying the needs of customers, increasing purchasing ability, profits and competitive advantages. The method applied in this study is Modified Image Method that is a common method to solve variational inequality. It assumes that K is a possible space and F is a function that acts in monotonous and Leap Sheets conditions. The method to solve inequality is proved and to apply it generally it was used on different issues with specified capacities. The result of proving the method of solution showed that it is an appropriate and reliable method to solve problems of chain supply management.
##### High Performance Computing of Meshless Time Domain Method on Multi-GPU Cluster

Soichiro Ikuno, Tokyo University of Technology, School of Computer Science, Japan

Yuta Hirokawa, University of Tsukuba, Graduate School of Systems and Information Engineering, Japan

Taku Itoh, Tokyo University of Technology, , Japan

Susumu Nakata, Ritsumeikan University, College of Information Science and Engineering, Japan

Abstract: High performance computing of Meshless Time Domain Method (MTDM) on Multi-GPU cluster is numerically investigated, and the electromagnetic wave propagation simulation is achieved by using parallelized MTDM on Multi-GPU. Generally, Finite Difference Time Domain (FDTD) Method is adopted for the simulation of the electromagnetic wave propagation phenomena. FDTD has great advantages in terms of parallelization and the scheme of FDTD is a transparent algorithm. However, a numerical domain must be divided in to rectangle meshes, and it is difficult to treat the problem in the complexed shaped domain. On the other hand, a meshless approach does not require meshes of a geometrical structure, and various meshless approaches such as the Radial Point Interpolation Method (RPIM) have been developed. In particular, meshless approaches based on RPIM are applied to time dependent problems, and the method is called Meshless Time Domain Method. The Graphics Processing Unit (GPU) is one of the most progressive device in recent years, and various researches of General Purpose computing on GPU (GPGPU) have been proposed aggressively. The results of computation show that the execution time of the time evolution calculation on GPU with various parallelization techniques is about 54 time faster than that of serial CPU.
##### Organization model for Mobile Wireless Sensor Networks (MWSN) inspired in Artificial Bee Colony (ABC)

Guilherme Roberto, São Paulo State University (UNESP), , Brazil

Luis Maschi, Federal Institute of São Paulo (IFSP) - Campus Catanduva, Informática, Brazil

Daniel Fernando Pigatto, Institute of Mathematics and Computer Sciences (ICMC) / University of São Paulo (USP), , Brazil

Carlos Montez, Federal University of Santa Catarina (UFSC), , Brazil

Leandro Neves, São Paulo State University, DCCE, Brazil

Kalinka Branco, Institute of Mathematics and Computer Sciences (ICMC) / University of São Paulo (USP), , Brazil

Alex Sandro Pinto, Universidade Federal de Santa Catarina, Câmpus de Blumenau, Brazil

Abstract: Wireless sensor networks (WSN) are a special type of computer networks, which basically aim at monitoring and tracking a specific area, typically with military, environmental or industrial purposes. These networks are considered a type of Ad Hoc network, and are composed of small computational nodes usually with limited processing capacity, capable of environment sensing and messages routing through radio frequency. A major problem in this type of networks is the energy consumption required for nodes operations, usually derived from batteries. The exchange of messages between nodes and data collection stations is the largest energy consumer within a typical WSN. Mobile wireless sensor networks (MWSN) are an interesting alternative to solve the problem of communication range, which usually is subjected to the limitations of a WSN. However, due to the mobility of nodes in a MWSN, these networks require frequent reformulations in route for forwarding messages to the base. Considering all the aspects mentioned above, the purpose of this study is to find a self-organizing model for MWSN based on bee colonies in order to reduce the number of messages transmitted among nodes, and thus reduce the overall consumption energy while maintaining the efficiency of message delivery. The results obtained in this article are originated from simulations carried out with SINALGO software, which demonstrates the effectiveness of the proposed approach. The BeeAODV (Bee Ad-Hoc On Demand Distance Vector) proposed in this paper allows to considerably reduce message exchanges whether compared to AODV (Ad-Hoc On Demand Distance Vector).
##### First-principle studies of lattice and electronic structure of BexZn1-xO

X.Lei, Inner Mongolia University, Physical Science and Technology, China

Guo Jun Zhao, Inner Mongolia University, Department of Physics, China

Xi Xia Liang, Inner Mongolia University, Department of Physics, School of Physical Science and Technology, , China

Tie Lei Song, Department of Physics, School of Physical Science and Technology, Inner Mongolia University, , China

Abstract: With the development of optoelectronics industry, II-VI wide-gap semiconductor zinc oxide has become a great importance to the photoelectric material, and it was used widely in ultraviolet light photoelectric devices. The causes of study ZnO alloy is due to its easy to doping ternary alloy formation mechanism and it is easy to meet the needs of a variety of optoelectronic devices. In this paper, beryllium doped zinc oxide(BexZn1-xO) ternary alloy was studied based on the density functional theory. After lattice structure and atomic position optimization for different doping concentration BexZn1-xO (x=0, 0.125, 0.25, 0.375, 0.5, 0.625, 0.75, 0.875 and 1, respectively), the lattice parameters for the stable structure and the band gap of BexZn1-xO are obtained. It turned out that the lattice constants a and lattice constants c decrease linearly as Be doping concentration increases, compliance with Vegard's law, the lattice parameters of the BexZn1-xO ternary alloys are consistent with experimental results. Band gap increases with increasing Be content, and the band gap values are fixed. Furthermore, the density of states of BexZn1-xO is calculated, the results shown that the valence band maximum (VBM) occupied by O2p states, and the conduction band minimum (CBM) occupied by Zn4s states. As Be doping, Zn 4s states dominate the conduction band and the conduction band bottom position constantly moving to higher energy region. Density of states strength Zn 4s states with decreasing proportion of Zn is constantly reduced.
##### Mathematical modeling of high-strain-rate deformation, kinetics of alpha-epsilon-phase transition, as well as nucleation, growth and recompaction of spall damage for recovered samples of quenched 30KhGSA steel under explosive loading

Natalia Zhilyaeva, Russian Federal Nuclear Center-Zababakhin Russian Science Research Institute of Technical Physics, , Russian Federation

Abstract: Advanced models for shear and spall strengths of the quenched 30KhGSA steel (HRc35...40) are described to take into account the reversing solid-state alpha-epsilon-phase transition, as well as the setup and results of earlier experiments with registration under two explosive loading modes of the two-, and three-wave configurations of stress waves in samples having different initial thickness. Consideration is given to the main theoretical predictions and simulation results for explosive experiments wherein wave profiles in the quenched 30KhGSA steel were registered with the help of streak-camera using the optical lever method and with the help of laser-interferometry using the Fabry-Perot and PDV methods. Comparing the calculated wave profiles with the measured ones, as well as comparing the predicted locations of the spall and shear micro-, meso-, and macrodamage zones with those observed on micro-, and macrophotographs of recovered samples cross-sections indicates that our calculations adequately describe wave processes in samples of the quenched 30KhGSA steel taking into account the elastic-plastic properties, direct and reverse alpha-epsilon-phase transition in stress waves, as well as kinetics of nucleation, growth, and recompaction of spall and shear damages.
##### DFT/TDDFT CALCULATIONS OF GROUND AND EXCITED STATES OF CIS AND TRANS 1, 3 INDANEDIOL

Moussaoui Yahia, Universite des sciences et de la technologie houari boumediene, Chimie physique, Algeria

Souhila Rezzouk, Université des Sciences et de la Technologie Houari Boumediene, Department of Chimie Physique, Algeria

Abstract: The system 1, 3-indanediol is presented under two isomers of cis and trans configurations the cis 1, 3-indanediol was used as a building unit in the synthesis of a new aromatic polyester obtained by poly-condensation. The ground state of cis and trans 1,3-indandiol were optimized by using the B3LYP DFT approach in combination with the basis set 6-31G*. The cis form have been found more stable than the trans configuration. This observation is in good agreement with experimental results. The excited states of our system were computed using Time dependent-Density Functional Theory (TD-DFT). References: [1] H. Guemmour, F. Carrière, A. Benaboura, Polym. Bull. 1-6 (2001) 46.. [2] E. Runge, E.K.U. Gross, Phys. Rev. Lett. 52, 997–1000 (1984). [3] E.K.U. Gross, C.A. Ullrich, U.J. Gossmann, Density Functional Theory, NATO ASI Series, Plenum, New York, pp. 149–171(1994). [4] R. van Leeuwen, Int. J. Mod. Phys. B 15, 1969 (2001). [5] M.A.L. Marques, E.K.U. Gross, Annu. Rev. Phys. Chem. 55, 427–455 (2004).
##### AB INITIO STUDY OF NEW COMPLEXES OF IMIDAZOLE AND SUBSTITUTED IMIDAZOLE WITH Zn ATOMS

Nait Achour Madjid, Universite des Sciences et de la Technologie Houari Boumediene, Chimie Physique, Algeria

Abstract: Imidazole and its derivatives are involved in important biological processes where they play a variety of roles [1]. Imidazole is also well known to adsorb and form strong complexes with transition metal ions[2]. Bonding between metal ions and imidazole ligands is very common in nature and many examples of metals in metalloproteins coordinated by imidazole have been reported [3]. Among these metals are found Zn2+ ions which are key structural components in a large number of metalloproteins. In the present work, a series of new compounds formed by imidazole ring and its derivatives with several Zn atoms, in which the first atom of Zn is formally in the +2 oxidation state and the others atoms are in their neutral state, are reported and studied theoretically. Thereby, both Density Functional Theory (DFT) and Moller Plesset (MP2) calculations were performed for these complexes. The elucidation of the structures of these complexes was made using the 6-311++G** function basis set. We evaluate the binding energy taking basis set superposition error (BSSE) and zero point vibrational energy (ZPVE) correction into account, as we calculate the IR spectrum, heats of formation and thermodynamic properties for all the obtained complexes. Harmonic vibrational frequencies calculations confirm that these elucidated structures are stable. All the calculations were performed using Gaussian 03 program. [1] (a) H.C. Freeman, in G.L. Eichhom (ed), Inorganic Biochemistry, Elsevier, New York, 1973, Ch.4 (b) R.J.Sundberg and R.B.Martin, Chem. Rev.,74(1974)471. [2]J.Hedin, D.Issaksson, M.Andersson, M.Nyden,J.of Colloid and Intrface Science 336(2009) 388-392. [3]J.M.Guss and H.C.Freeman, J.Mol.Biol.,169 (1983) 521.
##### Arbitrary l state solutions of the Feynman propagator with the Deng-Fan molecular potential

Ahmed Diaf, Laboratory of Energy and Smart Systems, Khemis Miliana University, Sciences and Technology, Algeria

Abstract: The bound state solutions of the Feynman propagator with the rotating Deng-Fan molecular potential are presented approximately. An approximation of the centrifugal potential is used and nonlinear space-time transformations are applied. A relation between the original path integral and the Green function of a new quantum soluble system is derived. The energy spectrum and the normalized eigenfunctions are both obtained for the application of this technique to the Deng-Fan molecular potential. Our results are in very good agreement with those found by using numerical and other methods.
##### Searching for AdS3 waves and Asymptotically Lifshitz black holes in R3-NMG

Giorgos Anastasiou, Universidad Andres Bello, Physics Department, Chile

Abstract: In this paper we consider the structure of the AdS3 vacua in R3 expansion of New Massive Gravity (R3-NMG). We obtain the degeneracies of the AdS3 vacua at several points of the parametric space. Additionally, following a speci c analysis we show that AdS3 wave solutions are present. Using these wave solutions, we single out two special points of the parametric space for which logarithmic terms appear in the solutions. The rst one is a point at which the e ective mass of the wave pro le, which is interpreted as a scalar mode, completely saturates the Breitenlohner- Freedman bound of the AdS3 space in which the wave is propagating. The second special point is a point at which the central charge of the theory vanishes. Furthermore, we investigate the possibility of asymptotically Lifshitz black hole solutions to be present in the three-dimensional R3-NMG. We derive analytically the Lifshitz vacua considering speci c relations between the mass parameters of R3-NMG. A certain polynomial equation arises at the rst special point where solutions with logarithmic fallo in the AdS3 space appear. Solving this polynomial equation, we obtain the values of the dynamical exponent z which correspond to possible asymptotically Lifshitz black hole solutions. However, it is shown that asymptotically Lifshitz black hole solutions do not exist in the three-dimensional R3-NMG for a speci c ansatz of the black hole metric.
##### Monte-Carlo investigation of proton pencil beam effects in the rectangular target, with using GEANT4 code

Seyed Ali Mahdipour, Hakim Sabzevari University, , Iran (Islamic Republic of)

Abstract: Radiotherapy with hadron beams like proton have been used for treatment of different cancers for many years. Determination of the particle’s deposited dose in the different tissues of budy is very important in radiotherapy and medical physics. GEANT4 Monte Carlo code reused for considering this effect. Simulation results help us to determine the physical parameters exactly. In this article we have calculate the total depth and lateral dose profiles for the 110 - 220 MeV proton pencil beams in the arbitrary rectangular target made of skin, adipose and tumor tissues. The material of tissues are taken from compositions of the ICRU 46. The results of this Monte Carlo code has a good agreement experiment data. moreover 2D depth dose profiles for the proton beam has been calculated in the phantom.
##### Meteorological time series forecasting based on MLP modeling using heterogeneous transfer functions

Voyant Cyril, University of Corsica, , France

Abstract: In this paper we propose to study four meteorological and seasonal time series coupled with multi-layer perceptron (MLP) modeling. We choose to combine two transfer functions in one heterogeneous function and to use a temporal indicator (time index as input) in order to take into account the seasonal aspect of the studied time series. The results of the prediction concern two years of never-used data (8 years for the learning step) and show that this methodology can improve the accuracy of meteorological data estimation compared to a classical MLP modeling with homogenous transfer function. The time index is the tool giving the best results for all cases, while in 75% of the cases it is for the use of the heterogeneous transfer function. For the time series with low variation coefficient, the lowest nRMSE is under 10% whereas for the other (solar irradiation and wind speed) the nRMSE is close to 30%.
##### AB INITIO STUDY OF NEW COMPLEXES OF IMIDAZOLE AND SUBSTITUTED IMIDAZOLE WITH Zn ATOMS

Nait Achour Madjid, Universite des Sciences et de la Technologie Houari Boumediene, Chimie Physique, Algeria

Abstract: Imidazole and its derivatives are involved in important biological processes where they play a variety of roles [1]. Imidazole is also well known to adsorb and form strong complexes with transition metal ions[2]. Bonding between metal ions and imidazole ligands is very common in nature and many examples of metals in metalloproteins coordinated by imidazole have been reported [3]. Among these metals are found Zn2+ ions which are key structural components in a large number of metalloproteins. In the present work, a series of new compounds formed by imidazole ring and its derivatives with several Zn atoms, in which the first atom of Zn is formally in the +2 oxidation state and the others atoms are in their neutral state, are reported and studied theoretically. Thereby, both Density Functional Theory (DFT) and Moller Plesset (MP2) calculations were performed for these complexes. The elucidation of the structures of these complexes was made using the 6-311++G** function basis set. We evaluate the binding energy taking basis set superposition error (BSSE) and zero point vibrational energy (ZPVE) correction into account, as we calculate the IR spectrum, heats of formation and thermodynamic properties for all the obtained complexes. Harmonic vibrational frequencies calculations confirm that these elucidated structures are stable. All the calculations were performed using Gaussian 03 program. [1] (a) H.C. Freeman, in G.L. Eichhom (ed), Inorganic Biochemistry, Elsevier, New York, 1973, Ch.4 (b) R.J.Sundberg and R.B.Martin, Chem. Rev.,74(1974)471. [2]J.Hedin, D.Issaksson, M.Andersson, M.Nyden,J.of Colloid and Intrface Science 336(2009) 388-392. [3]J.M.Guss and H.C.Freeman, J.Mol.Biol.,169 (1983) 521.
##### Activation energy of disordering lattices superionic crystals LnF3 (Ln = La, Ce, Pr): experimental results and quantum-chemical calculations

Sirojiddin Mirzaev, Instituteof Ionic, Plazma and Lazer Technology, , Uzbekistan

Valeriy Krivorotov, Instituteof Ionic, Plazma and Lazer Technology, , Uzbekistan

Georgiy Nujdov, Institute of Ionic, Plazma and Lazer Technology , , Uzbekistan

Abstract: The lattice disordering in superionic crystal LaF3 investigated by inelastic light scattering spectroscopy and quantum-chemical calculations energy process. Based on analysis of the line widths of Raman scattering in the crystal LaF3 found that the values of the activation energy Ea = 0.18 eV F1 sublattice disordering in the insulating phase and 0.04 eV for the superconductive phase. Quantum-chemical calculations of the potential relief in the LaF3 cluster ions of 1200 proved that the energy Ea of defect vacancy- interstitial fluorine ion equal 0.16 eV when moving to an interstitial single ion F1 ( insulating phase ) and 0.03 eV for the movement of ions in the superionic phase F1 . On the basis of quasi-elastic light scattering in LaF3 shown that potential barriers Ed that impede the movement of ions F1, decrease from 0.28 eV ( insulating phase ) to 0.07 eV in the superionic phase . The results lead to an important conclusion for the physics of superionic conductors : Ed and Ea values ​​obtained for the insulating phase ( T
##### Mathematical Tools for Discovery of Nanoporous Materials for Energy Applications

Maciej Haranczyk, Lawrence Berkeley National Laboratory, Computational Research Division, United States

Richard Martin, Lawrence Berkeley National Laboratory, Computational Research Division, United States

Abstract: Porous materials such as zeolites and metal organic frameworks have been of growing importance as materials for energy-related applications such as CO2 capture, hydrogen and methane storage, and as catalysts. The current state-of-the-art molecular simulations allow for accurate in silico prediction of materials’ properties but the computational cost of such calculations prohibits their application in the characterization of very large sets of structures, which would be required to perform brute-force screening. Our work focuses on the development of novel methodology to harness this complexity of the materials space. In particular, we have been developing algorithms and tools for enumeration and characterization of porous material databases as well as efficient screening approaches. The latter include similarity-based, optimization-based approaches as well as application of structure-property relationships. Our methodology represents a “soup” of mathematical methods. We have used Voronoi tessellation-based techniques to enable high-throughput structure characterization and comparison, PDE-based techniques to predict guest-molecule accessibility, and continues and discrete optimization to design materials. The resulting hybrid material discovery suite requires expensive characterization only for carefully selected and statistically relevant subset of a database, therefore enabling discoveries at a minimal computational cost. Our presentation will give an overview of recent developments as well as highlight few interesting applications.
##### Firewalls: arising from vaccum?

Luis Cabarique, Universidad Nacional de Colombia, Departamento de Física, Colombia

Robel Arenas, Universidad Nacional de Colombia, Observatorio Astronómico Nacional, Colombia

Abstract: The appearance of a firewall near a black hole has caused great controversy. In this paper we discuss the black hole firewall centering in the vaccum properties of each corresponding observer in (1+1)-dimensions. From this properties we found that the horizon appears naturally as a massless radiating Shell, but it no necessarily is the same firewall proposed as a possible solution to an apparent inconsistency in black hole complementarity.
##### Numerical study of the propagation of train of chirped vector solitons in birefringent optical fibers with variable coefficients

Siham Aziez , Université Hadj-Lakhdar de Batna, Département de Science de la Matière, Faculté des Sciences, Algeria

Derradji Bahloul, Université El Hadj Lakhdar de Batna, Sciences de la Matière, Algeria

Abstract: We study numerically in this work the propagation of train of chirped vector solitons in birefringent optical fibers with variable coefficients using the compact split step Padé scheme (CSSPS) [5]. Generally, the propagation of managed vector solitons in birefringent optical fibers is governed by the coupled NLS equations with variable coefficients which are not integrable, but may possess exact solitary-wave solutions (vector solitons) in some particular cases.
##### Generalized Radon Inversion from an Efficient Number of Moment Based Projections Applied on Computerized Tomography

Spiros Chountasis, Independent Power Transmission Operator, Department of Systems & Infrastructure, Greece

Abstract: This paper presents a new method for image computerized tomography that is based on the usage of a generalized inversion that is rotation-invariant utilizing Radon and time-frequency transforms. The proposed technique employs the Fourier and Haar coefficients for spectral and spatial moment based image analysis, respectively. It provides a new approach to the problem of tomographic image reconstruction where an X-ray image obtained from a set of line projections. The experimental evaluation of the scheme leads to the conclusion that the increased selectivity of the method provides a faster and more robust approach to the problem. The method has been tested by reconstructing a commonly used tomography image at the presence of blur caused by uniform linear motions. The noise during normal transmission process is another issue that is considered in the current work.
##### Keto-enol tautomerism in 5- and 6-azauracils in water solution

Nadezhda Markova - Petrova, Institute of Organic Chemistry with Centre of Phytochemistry, Bulgarian Academy of Sciences, , Bulgaria

Abstract: Ab initio quantum chemical investigations on the tautomeric equilibrium in 5- and 6-azauracils in water were performed. The solvent effects were considered by explicit inclusion of two pairs of water molecules, which model a first hydration shell around the solute. The effects of the water as solvent were introduced at two different levels - using solute-solvent clusters (four water molecules) and using the same clusters embedded in an external continuum. Full geometry optimizations of these complexes were carried out at MP2/6-31+G(d,p) and C-PCM/MP2/6-31+G(d,p). Single point calculations were performed at CCSD(T)/6-31+G(d,p)//MP2/6-31+G(d,p) computational level to obtain accurate energies. According to our calculations, hydrated azauracils should exist in three forms: mainly dioxo form and two hydroxy forms. The calculated proton transfer activation energies for the tautomeric reactions of 5-azauracil and 6-azauracil show different pictures for these two compounds. According to C-PCM/MP2/6-31+G(d,p) data water-assisted proton transfer in 5-azauracil occurs through two parallel reactions A→B and A→D. Tautomeric equilibrium in 6-azauracil in water could occur by two contiguous reactions, A→D and D→C. The proton transfer investigated reactions in 5- and 6-azauracils involve concerted atomic movement.
##### Stochastic Multiscale Modeling of Physical and Biological Systems

George Karniadakis, Brown University, Applied Mathematiacs, United States

Abstract: TBA
##### The Enigma of Dark Energy

Ofer Lahav, University College London, , United Kingdom

Abstract: There is strong observational evidence that our Universe is flat and it consists of three main ingredients: ordinary matter, Dark Matter and Dark Energy. Dark Energy is commonly interpreted as the cause of the observed acceleration of the cosmic expansion. The lecture will review the chequered history of Dark Energy, the current observations and new surveys such as the “Dark Energy Survey”, and mathematical and statistical approaches.
##### Structure and Dynamics of Multi-scale Interacting Diseases

Yamir Moreno, Institute BIFI, University of Zaragoza, , Spain

Abstract: TBA
##### Hybrid Ab Initio/Molecular Mechanical Simulations of Solutions - Perspectives for Complex Chemistry, Proton Transfer Reactions and Interfaces

Thomas Hofer, University of Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, Austria

Abstract: Following the continous improvmement of computational resources, theoretical chemistry has proven to be an increasingly valuable approach for investigations of a large number of chemical phenomena, especially in case of liquids and solutions, which are ont one hand the most important but at the same time also one of the most challenging states of matter to describe. Hybrid quantum mechanical/molecular mechanical (QM/MM) simulation techniques are still considered as one of the most promising approaches for investigations of chemical system. These methods separate the system into two regions. While interactions in the chemical most relevant part are accounted for via quantum mechanics, empirical molecular mechanical potentials are sufficiently accurate to describe the remaining part of the system. Especially for the determination of experimentally elusive properties such as ultrafast phenomena or single ion properties QM/MM techniques proved to be particularly useful alternative approaches. Details of advanced QM/MM technqiues are presented and simulation results for various systems ranging from coordination complexes in solution to the description of proton transfer events will be presented. A particularly challenging application is the description of solid interfaces, which is a prerequisite to study of adsorption and related phenomena such as catalysis. The combination of the QM/MM methodology with a periodic quantum mechanical description of the system will presented. Exemplary results of a QM/MM simulation study of the MgO(100)/water interface demonstrate the capabilities of this approach.
##### Geometrodynamics in Cosmology

Spyros Basilakos, Academy of Athens, RCAAM, Greece

Abstract: In this review talk I will discuss some ideas based on Geometrodynamics towards explaining the accelerated expansion of the Universe. In particular, I will describe some basic modified gravity models namely f(R), Dvali, Gabadadze, and Porrati(DGP) braneworld, Gauss-Bonnet and Finsler, and I will discuss their Cosmological implications.
##### Polarons on nonlinear lattice in the Su–Schrieffer–Heeger approximation: Exact solution and multipeaked polarons.

George Vinogradov, Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Statistical physics, Russian Federation

Abstract: We investigate the polaron dynamics on the nonlinear lattice with the cubic nonlinearity. The electron-phonon interaction is accounted in the Su-Schrieffer-Heeger approximation. An exact analytical solution is obtained in the continuum approximation. The numerical simulation agrees with analytics very well. Moreover, colliding polarons recover their shapes and velocities after the elastic collision suggesting that the solution belongs to the exactly integrable system. When the continuum approximation is invalid (parameters of nonlinearity and electron-phonon interaction are not small), a new family of stable multipeaked polarons is found. These polarons are formed by the coupled solitons hold together by the electron-phonone interaction.

#### International Conference on Mathematical Modeling in Physical Sciences

August 28-31, 2014