We study ground states in the nonlinear Schrödinger equation (NLS) with an isotropic harmonic potential, in energy-critical and energy-supercritical cases. In both cases, we prove existence of a family of ground states parametrized by their amplitude, together with the corresponding values of the spectral parameter. Moreover, we derive asymptotic behavior of the spectral parameter when the amplitude of ground states tends to infinity. We show that in the energy-supercritical case the family of ground states converges to a limiting singular solution and the spectral parameter converges to a nonzero limit, where the convergence is oscillatory for smaller dimensions, and monotone for larger dimensions. In the energy-critical case, we show that the spectral parameter converges to zero, with a specific leading-order term behavior depending on the spatial dimension. Furthermore, we study the Morse index of the ground states in the energy-supercritical case. We show that in the case of monotone behavior of the spectral parameter, that is for large values of the dimension, the Morse index of the ground state is finite and independent of its amplitude. Moreover, we show that it asymptotically equals to the Morse index of the limiting singular solution. This result suggests how to estimate the Morse index of the ground state numerically. / Dissertation / Doctor of Philosophy (PhD)

The purpose of this thesis is to study the effects of gain and damping on a nonlinear waveguide array with a strong disorder that is driven in the first site, and try to find regimes which have stable stationary solutions. This has been done with a modified DNLS (Discrete nonlinear Schrödinger equation). Stable stationary solutions were mainly found when the damping was stronger than the gain, but some stable stationary regimes were also found when the gain was stronger than the damping.

This thesis deals with problems arising in the study of nonlinear partial differential equations arising from many-body problems. It is divided into two parts: The first part concerns the derivation of a nonlinear diffusion equation from a microscopic stochastic process. We give a new method to show that in the hydrodynamic limit, the particle densities of a one-dimensional zero range process on a periodic lattice converge to the solution of a nonlinear diffusion equation. This method allows for the first time an explicit uniform-in-time bound on the rate of convergence in the hydrodynamic limit. We also discuss how to extend this method to the multi-dimensional case. Furthermore we present an argument, which seems to be new in the context of hydrodynamic limits, how to deduce the convergence of the microscopic entropy and Fisher information towards the corresponding macroscopic quantities from the validity of the hydrodynamic limit and the initial convergence of the entropy. The second part deals with problems arising in the analysis of nonlinear Schrödinger equations of Gross-Pitaevskii type. First, we consider the Cauchy problem for (energy-subcritical) nonlinear Schrödinger equations with sub-quadratic external potentials and an additional angular momentum rotation term. This equation is a well-known model for superfluid quantum gases in rotating traps. We prove global existence (in the energy space) for defocusing nonlinearities without any restriction on the rotation frequency, generalizing earlier results given in the literature. Moreover, we find that the rotation term has a considerable influence in proving finite time blow-up in the focusing case. Finally, a mathematical framework for optimal bilinear control of nonlinear Schrödinger equations arising in the description of Bose-Einstein condensates is presented. The obtained results generalize earlier efforts found in the literature in several aspects. In particular, the cost induced by the physical work load over the control process is taken into account rather then often used L^2- or H^1-norms for the cost of the control action. We prove well-posedness of the problem and existence of an optimal control. In addition, the first order optimality system is rigorously derived. Also a numerical solution method is proposed, which is based on a Newton type iteration, and used to solve several coherent quantum control problems.

京都大学 / 0048 / 新制・課程博士 / 博士(情報学) / 甲第22863号 / 情博第742号 / 新制||情||127(附属図書館) / 京都大学大学院情報学研究科数理工学専攻 / (主査)教授 矢ヶ崎 一幸, 教授 中村 佳正, 准教授 柴山 允瑠, 教授 國府 寛司 / 学位規則第4条第1項該当 / Doctor of Informatics / Kyoto University / DFAM

Submitted by Luciana Ferreira (lucgeral@gmail.com) on 2018-07-17T13:07:11Z No. of bitstreams: 2 Tese - Luciano Calaça Alves - 2018.pdf: 4699371 bytes, checksum: 706846314ebb74648be890b03e18d4cc (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Approved for entry into archive by Luciana Ferreira (lucgeral@gmail.com) on 2018-07-17T13:55:46Z (GMT) No. of bitstreams: 2 Tese - Luciano Calaça Alves - 2018.pdf: 4699371 bytes, checksum: 706846314ebb74648be890b03e18d4cc (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Made available in DSpace on 2018-07-17T13:55:46Z (GMT). No. of bitstreams: 2 Tese - Luciano Calaça Alves - 2018.pdf: 4699371 bytes, checksum: 706846314ebb74648be890b03e18d4cc (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Previous issue date: 2018-06-07 / Fundação de Amparo à Pesquisa do Estado de Goiás - FAPEG / This work presents the study of solitary wave solutions, known as solitons, in non-linear and non- homogeneous media using non-linear Schrödinger equations. Three cases are studied: first considering a logarithmic nonlinear term; second with saturation effect and finally including effects of high orders (Raman scattering). Solutions are modulated by three different types of potential. First, linear in the spatial and oscillatory coordinate in the temporal coordinate. The second, quadratic in the spatial and oscillatory in the temporal coordinates. Finally, it is also modulated using a mixed potential, which is the junction of the two potentials presented above. After including inomogeneities in linear and nonlinear coefficients, the similarity transformation technique is used to convert the non-linear, non-autonomous equation into an autonomous one that will be solved analytically. This field of study has potential applications in crystals, optical fibers and in Bose- Einstein condensates, also serving to understand the fundamentals related to this state of matter. The stability of the solutions are checked by numerical simulations. / Este trabalho apresenta o estudo de soluções de ondas solitárias, conhecidas como sólitons, em meios não lineares e não homogêneos por meio de equações não lineares de Schrödinger. São estudados três casos: primeiro considerando um termo não linear do tipo logarítmico; segundo com efeito de saturação e por último incluindo efeitos de altas ordens (espalhamento Raman). As soluções são moduladas por três tipos diferentes de potencial. O primeiro, linear na coordenada espacial e oscilatório na coordenada temporal. O segundo, quadrático na coordenada espacial e oscilatório na temporal. Por fim, modula-se também utilizando um potencial misto, que é a junção dos dois potenciais apresentados anteriormente. Depois de incluir heterogeneidades nos coeficientes lineares e não lineares, é utilizada a técnica da transformação de similaridade para converter a equação não linear, não autônoma em uma autônoma que será resolvida analiticamente. Esse campo de estudo tem potenciais aplicações em cristais, fibras ópticas e em condensados de Bose-Einstein, servindo também para o entendimento dos fundamentos relacionados a esse estado da matéria. A estabilidade das soluções são checadas por meio de simulações numéricas.

Optical fibers play a vital role in modern telecommunication systems and networks. An optical fiber link imposes some linear and nonlinear distortions on the propagating light-wave signal due to the inherent dispersive nature and nonlinear behavior of the fiber. These distortions impede the increasing demand for higher data rate transmission over longer distances. Developing efficient and computationally non-expensive digital signal processing (DSP) techniques to effectively compensate for the fiber impairments is therefore essential and of preeminent importance. This thesis proposes two DSP-based approaches for mitigating the induced distortions in short-reach and long-haul fiber-optic communication systems. The first approach introduces a powerful digital nonlinear feed-forward equalizer (NFFE), exploiting multilayer artificial neural network (ANN). The proposed ANN-NFFE mitigates nonlinear impairments of short-haul optical fiber communication systems, arising due to the nonlinearity introduced by direct photo-detection. In a direct detection system, the detection process is nonlinear due to the fact that the photo-current is proportional to the absolute square of the electric field intensity. The proposed equalizer provides the most efficient computational cost with high equalization performance. Its performance is comparable to the benchmark compensation performance achieved by maximum-likelihood sequence estimator. The equalizer trains an ANN to act as a nonlinear filter whose impulse response removes the intersymbol interference (ISI) distortions of the optical channel. Owing to the proposed extensive training of the equalizer, it achieves the ultimate performance limit of any feed-forward equalizer. The performance and efficiency of the equalizer are investigated by applying it to various practical short-reach fiber-optic transmission system scenarios. These scenarios are extracted from practical metro/media access networks and data center applications. The obtained results show that the ANN-NFFE compensates for the received BER degradation and significantly increases the tolerance to the chromatic dispersion distortion. The second approach is devoted for blindly combating impairments of long-haul fiber-optic systems and networks. A novel adjoint sensitivity analysis (ASA) approach for the nonlinear Schrödinger equation (NLSE) is proposed. The NLSE describes the light-wave propagation in optical fiber communication systems. The proposed ASA approach significantly accelerates the sensitivity calculations in any fiber-optic design problem. Using only one extra adjoint system simulation, all the sensitivities of a general objective function with respect to all fiber design parameters are estimated. We provide a full description of the solution to the derived adjoint problem. The accuracy and efficiency of our proposed algorithm are investigated through a comparison with the accurate but computationally expensive central finite-differences (CFD) approach. Numerical simulation results show that the proposed ASA algorithm has the same accuracy as the CFD approach but with a much lower computational cost. Moreover, we propose an efficient, robust, and accelerated adaptive digital back propagation (A-DBP) method based on adjoint optimization technique. Provided that the total transmission distance is known, the proposed A-DBP algorithm blindly compensates for the linear and nonlinear distortions of point-to-point long-reach optical fiber transmission systems or multi-point optical fiber transmission networks, without knowing the launch power and channel parameters. The NLSE-based ASA approach is extended for the sensitivity analysis of general multi-span DBP model. A modified split-step Fourier scheme method is introduced to solve the adjoint problem, and a complete analysis of its computational complexity is studied. An adjoint-based optimization (ABO) technique is introduced to significantly accelerate the parameters extraction of the A-DBP. The ABO algorithm utilizes a sequential quadratic programming (SQP) technique coupled with the extended ASA algorithm to rapidly solve the A-DBP training problem and optimize the design parameters using minimum overhead of extra system simulations. Regardless of the number of A-DBP design parameters, the derivatives of the training objective function with respect to all parameters are estimated using only one extra adjoint system simulation per optimization iterate. This is contrasted with the traditional finite-difference (FD)-based optimization methods whose sensitivity analysis calculations cost per iterate scales linearly with the number of parameters. The robustness, performance, and efficiency of the proposed A-DBP algorithm are demonstrated through applying it to mitigate the distortions of a 4-span optical fiber communication system scenario. Our results show that the proposed A-DBP achieves the optimal compensation performance obtained using an ideal fine-mesh DBP scheme utilizing the correct channel parameters. Compared to A-DBPs trained using SQP algorithms based on forward, backward, and central FD approaches, the proposed ABO algorithm trains the A-DBP with 2.02 times faster than the backward/forward FD-based optimizers, and with 3.63 times faster than the more accurate CFD-based optimizer. The achieved gain further increases as the number of design parameters increases. A coarse-mesh A-DBP with less number of spans is also adopted to significantly reduce the computational complexity, achieving compensation performance higher than that obtained using the coarse-mesh DBP with full number of spans. / Thesis / Doctor of Philosophy (PhD) / This thesis proposes two powerful and computationally efficient digital signal processing (DSP)-based techniques, namely, artificial neural network nonlinear feed forward equalizer (ANN-NFFE) and adaptive digital back propagation (A-DBP) equalizer, for mitigating the induced distortions in short-reach and long-haul fiber-optic communication systems, respectively. The ANN-NFFE combats nonlinear impairments of direct-detected short-haul optical fiber communication systems, achieving compensation performance comparable to the benchmark performance obtained using maximum-likelihood sequence estimator with much lower computational cost. A novel adjoint sensitivity analysis (ASA) approach is proposed to significantly accelerate sensitivity analyses of fiber-optic design problems. The A-DBP exploits a gradient-based optimization method coupled with the ASA algorithm to blindly compensate for the distortions of coherent-detected fiber-optic communication systems and networks, utilizing the minimum possible overhead of performed system simulations. The robustness and efficiency of the proposed equalizers are demonstrated using numerical simulations of varied examples extracted from practical optical fiber communication systems scenarios.

We clear up two issues regarding the eigenvalue problem for the Manakov system; these problems relate directly to the existence of the soliton [sic] effect in fiber optic cables. The first issue is a bound on the eigenvalues of the Manakov system: if the parameter ξ is an eigenvalue, then it must lie in a certain region in the complex plane. The second issue has to do with a chirped Manakov system. We show that if a system is chirped too much, the soliton effect disappears. While this has been known for some time experimentally, there has not yet been a theoretical result along these lines for the Manakov system. / Ph. D.

Cette thèse étudie l'influence du vent sur l'instabilité modulationnelle. Une première partie unifie les travaux de Segur et al. qui intègrent la dissipation et ceux de Leblanc qui prennent en compte le vent. Une équation non linéaire de Schrödinger est établie avec un terme additionnel linéaire résultant de la compétition entre le vent et la dissipation. La dissipation est traduite par le modèle de Lundgren et l'effet du vent se manifeste par l'intermédiaire de la pression atmosphérique selon le modèle de Miles. La profondeur est finie. Une étude de stabilité de l'onde de Stokes est détaillée, et des simulations numériques sont menées pour illustrer les résultats. Des expérimentations sont menées pour apporter une validation qualitative à ces travaux. Cette première partie a été validée par une publication au Journal of Fluid Mechanics (2010). La deuxième partie étudie l'influence du vent sur l'instabilité modulationnelle par l'intermédiaire de la vorticité qu'il crée en surface. Le modèle est simplifié par l'hypothèse d'un écoulement unidirectionnel et d'une vorticité constante. La profondeur est encore supposée finie. Une équation non linéaire de Schrödinger est établie, qui prend en compte cette vorticité constante. La stabilité de l'onde de Stokes est alors étudiée en détail(diagramme d'instabilité en fonction de la vorticité et de la profondeur, bande d'instabilité, taux d'instabilité, etc.). Il est démontré qu'une vorticité négative, au delà d'un certain seuil, supprime l'instabilité modulationnelle indépendamment de la profondeur. Cette deuxième partie a été soumise pour publication au journal Physics of Fluids. / This thesis manuscript treats about the influence of wind on modulational instability. A first part merges the works of Segur at al. which take into account viscous dissipation and Leblanc's work which deals with wind. A nonlinear Schrödinger equation is derived, with a forcing linear term which represents the result of the balance between wind forcing and dissipation. Visous dissipation is represented by Lundgren's model and the effect of wind is integrated into atmospheric pressure following Miles' model. Depth is finite. The stability of Stokes's waves is investigated, and numerical simulations are presented to illustrate the results. Some experimentations are done to confirm qualitatively these works. This first part was validated by a publication in the Journal of Fluid Mechanics~(2010). The second part studies the influence of the wind on the modulational instability by the intermediary of the vorticity whom it creates on the water at the surface. The model is simplified by the hypothesis of an unidirectional flow and a constant vorticity. The depth is still supposed finite. A non linear Schrödinger equation is derived, which takes into account this constant vorticity. The stability of the Stokes' wave is studied then in detail (instability diagram function of vorticity and depth, instability bandwidth, instability rate, etc.). It is demonstrated that a negative vorticity, beyond a certain threshold, eliminates the modulational instability independently of the depth. This second part has been submitted for publication in the journal Physics of Fluids.

Ce travail cherche a comprendre comment l'ajout d'une géométrie non euclidienne dans un problème de Schrödinger non linéaire influe sur l'existence et l'unicité des solutions explosives de masse critique. On s'inspire pour beaucoup des travaux de Merle et Raphaël sur la méthode de modulation des paramètres d'invariance géométrique pour une EDP qui possède de bonnes lois de conservations. On s'appuie ici plus particulièrement sur un article de Raphaël et Szeftel qui prouve l'existence et l'unicité d'une solution de masse critique en dimension 2 pour l'équation de Schrödinger non linéaire avec potentiel d'inhomogénéité devant la non-linéarité, et qui explose par ailleurs au maximum de l'inhomogénéité. Dans un premier temps, il s'agit de reprendre la méthode dans son ensemble afin de l'adapter à des cas où le Laplacien n'est plus plat, et est remplacé par un opérateur de type Laplace-Beltrami ou Laplacien généralisé. Ayant mis en avant le rôle de la courbure au point d'explosion, en termes de conditions sur les dérivées de termes métriques, on reprend dans un deuxième temps l'étude dans le cas plus général d'une variété riemannienne. Grâce à un ansatz sur la solution qui intègre maintenant la transformation induite par la métrique, on est capable d'énoncer un résultat d'existence et d'unicité en termes de conditions géométriques sur la variété elle même. Par soucis de simplicité, on se limite néanmoins au rôle local de la métrique, en la supposant globalement définie dans une certaine carte, et asymptotiquement équivalente a la métrique euclidienne. / The present work aims at investigating the effects of a non-euclidean geometry on existence and uniqueness results for critical blow up NLS solutions. We will use many ideas from the works of Merle and Raphaël, particularly ideas from modulation theory which describes a solution in terms of geometric invariants parameters. We will rely more specically on a paper from Raphaël and Szeftel for existence and uniqueness of a critical mass blow up solution in dimension two tothe nonlinear Schrödinger equation with inhomogeneous potential acting on the nonlinearity, and which blows up where the inhomogeneity reaches its maximum. At first, we consider a generalized Laplacian operator and deploy the classical ansatz method to point out difficulties inherited from the non-flat metric terms, and in particular the key role played by the curvature at the blow-up point. In a second part, we reproduce the method when modifying the geometrical ansatz on which the parametrix is constructed, and investigate more precisely what is needed for existence and then uniqueness when dealing with a Laplace-Beltrami operator associated to a riemannian manifold. For simplicity, we shall only consider the role of g locally around the blow up point we are constructing, by assuming g is globally defined in some map, and asymptotically equals the usual euclidean metric.

<p>Numerical methods and algorithms for solving of partial differential equations, especially parallel algorithms, are an important research topic, given the very broad applicability range in all areas of science. Rapid advances of computer technology open up new possibilities for development of faster algorithms and numerical simulations of higher resolution. This is achieved through paralleliza-tion at different levels that&nbsp; practically all current computers support.</p><p>In this thesis we develop parallel algorithms for solving one kind of partial differential equations known as nonlinear Schr&ouml;dinger equation (NLSE) with a convolution integral kernel. Equations of this type arise in many fields of physics such as nonlinear optics, plasma physics and physics of ultracold atoms, as well as economics and quantitative&nbsp; finance. We focus on a special type of NLSE, the dipolar Gross-Pitaevskii equation (GPE), which characterizes the behavior of ultracold atoms in the state of Bose-Einstein condensation.</p><p>We present novel parallel algorithms for numerically solving GPE for a wide range of modern parallel computing platforms, from shared memory systems and dedicated hardware accelerators in the form of graphics processing units (GPUs), to&nbsp;&nbsp; heterogeneous computer clusters. For shared memory systems, we provide an algorithm and implementation targeting multi-core processors us-ing OpenMP. We also extend the algorithm to GPUs using CUDA toolkit and combine the OpenMP and CUDA approaches into a hybrid, heterogeneous al-gorithm that is capable of utilizing all&nbsp; available resources on a single computer. Given the inherent memory limitation a single&nbsp; computer has, we develop a distributed memory algorithm based on Message Passing Interface (MPI) and previous shared memory approaches. To maximize the performance of hybrid implementations, we optimize the parameters governing the distribution of data&nbsp; and workload using a genetic algorithm. Visualization of the increased volume of output data, enabled by the efficiency of newly developed algorithms, represents a challenge in itself. To address this, we integrate the implementations with the state-of-the-art visualization tool (VisIt), and use it to study two use-cases which demonstrate how the developed programs can be applied to simulate real-world systems.</p> / <p>Numerički metodi i algoritmi za re&scaron;avanje parcijalnih diferencijalnih jednačina, naročito paralelni algoritmi, predstavljaju izuzetno značajnu oblast istraživanja, uzimajući u obzir veoma &scaron;iroku primenljivost u svim oblastima nauke. Veliki napredak informacione tehnologije otvara nove mogućnosti za razvoj bržih al-goritama i&nbsp; numeričkih simulacija visoke rezolucije. Ovo se ostvaruje kroz para-lelizaciju na različitim nivoima koju poseduju praktično svi moderni računari. U ovoj tezi razvijeni su paralelni algoritmi za re&scaron;avanje jedne vrste parcijalnih diferencijalnih jednačina poznate kao nelinearna &Scaron;redingerova jednačina sa inte-gralnim konvolucionim kernelom. Jednačine ovog tipa se javljaju u raznim oblas-tima fizike poput nelinearne optike, fizike plazme i fizike ultrahladnih atoma, kao i u ekonomiji i kvantitativnim finansijama. Teza se bavi posebnim oblikom nelinearne &Scaron;redingerove jednačine, Gros-Pitaevski jednačinom sa dipol-dipol in-terakcionim članom, koja karakteri&scaron;e pona&scaron;anje ultrahladnih atoma u stanju Boze-Ajn&scaron;tajn kondenzacije.<br />U tezi su predstavljeni novi paralelni algoritmi za numeričko re&scaron;avanje Gros-Pitaevski jednačine za &scaron;irok spektar modernih računarskih platformi, od sis-tema sa deljenom memorijom i specijalizovanih hardverskih akceleratora u ob-liku grafičkih procesora, do heterogenih računarskih klastera. Za sisteme sa deljenom memorijom, razvijen je&nbsp; algoritam i implementacija namenjena vi&scaron;e-jezgarnim centralnim procesorima&nbsp; kori&scaron;ćenjem OpenMP tehnologije. Ovaj al-goritam je pro&scaron;iren tako da radi i u&nbsp; okruženju grafičkih procesora kori&scaron;ćenjem CUDA alata, a takođe je razvijen i&nbsp; predstavljen hibridni, heterogeni algoritam koji kombinuje OpenMP i CUDA pristupe i koji je u stanju da iskoristi sve raspoložive resurse jednog računara.<br />Imajući u vidu inherentna ograničenja raspoložive memorije koju pojedinačan računar poseduje, razvijen je i algoritam za sisteme sa distribuiranom memorijom zasnovan na Message Passing Interface tehnologiji i prethodnim algoritmima za sisteme sa deljenom memorijom. Da bi se maksimalizovale performanse razvijenih hibridnih implementacija, parametri koji određuju raspodelu podataka i računskog opterećenja su optimizovani kori&scaron;ćenjem genetskog algoritma. Poseban izazov je vizualizacija povećane količine izlaznih podataka, koji nastaju kao rezultat efikasnosti novorazvijenih algoritama. Ovo je u tezi re&scaron;eno kroz inte-graciju implementacija sa najsavremenijim alatom za vizualizaciju (VisIt), &scaron;to je omogućilo proučavanje dva primera koji pokazuju kako razvijeni programi mogu da se iskoriste za simulacije realnih sistema.</p>