Depleted-Wave Electric-Field Solver for Pulsed Bessel-Gaussian Type-II Frequency Doubling in KTP Crystals

Depleted-Wave Electric-Field Solver for Pulsed Bessel-Gaussian Type-II Frequency Doubling in KTP Crystals | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Depleted-Wave Electric-Field Solver for Pulsed Bessel-Gaussian Type-II Frequency Doubling in KTP Crystals Mohammad Sabaeian, Alireza Motazedian, Mostafa M. Rezaee, Fatemeh Sedaghat Jalil-Abadi, and 1 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-9203079/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract This toolkit addresses the computational modeling of three-dimensional time-dependent electric-field distributions in type-II pulsed Bessel-Gaussian second-harmonic generation within potassium titanyl phosphate (KTP) crystals, where depleted-wave dynamics dominate nonlinear optical conversion. The toolkit solves three coupled nonlinear wave equations for ordinary and extraordinary fundamental waves and extraordinary second-harmonic waves using finite difference methods in cylindrical coordinates. It accepts user-defined parameters including pulse energy, beam spot size, crystal dimensions, and optical properties as inputs, and produces spatiotemporal electric-field amplitude distributions and conversion efficiency profiles as outputs. A unified Fortran codebase provides reproducible simulation pipelines, parametric sweep capabilities for examining interaction length dependencies, and computational optimizations that enable execution on standard desktop systems despite the fine radial meshes required for Bessel-Gaussian beam fluctuations. The implementation reproduces previously published quantitative predictions, specifically confirming that for ωf ≈ 80 μm spot sizes and 0.8 J pulse energies, complete energy exchange occurs over approximately 5 mm, which validates the necessity of depleted-wave formalism for crystals exceeding this interaction length. The toolkit is available as an open-source GitHub repository and is released as version v1.0.1 under the MIT License, with a permanently archived distribution identified by DOI 10.5281/zenodo.17362587. Optics/Lasers Bessel-Gaussian beams second-harmonic generation type-II phase matching depleted-wave model KTP pulsed SHG Full Text Additional Declarations The authors declare no competing interests. Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-9203079","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":610862634,"identity":"aa936c7a-6db3-4b59-a3e5-da0770d4307d","order_by":0,"name":"Mohammad Sabaeian","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAx0lEQVRIiWNgGAWjYFACHiA2YGDgh/HZiNYi2UCaFpCuA8Q6S7797MEPHwrs8ozPH34mwVBjx8AnTUCzwZm8ZMkZBsnFZjfSzCQYjiUzsPElENDCkGMgzWPAnLjtBg+bBAPbAQY2HkIO639j/PuPQX3i5v4zQC3/iNDCcCPHTJrB4HDiBoYcNgnGNiK0GNx4l2bZY3A8ccaNNGOLxL5kHiIclnv4xo8/1Yn9/Ycf3vjwzU5OvoeQw1BAAjyaRsEoGAWjYBRQBAA5hjmolQ8MQwAAAABJRU5ErkJggg==","orcid":"","institution":"Department of Physics, Shahid Chamran University of Ahvaz, Ahvaz, Khuzestan, Iran","correspondingAuthor":true,"prefix":"","firstName":"Mohammad","middleName":"","lastName":"Sabaeian","suffix":""},{"id":610862635,"identity":"5571ca1e-c825-4b9b-898a-b9d57b394ba0","order_by":1,"name":"Alireza Motazedian","email":"","orcid":"","institution":"Department of Physics, University of New Hampshire, Durham, NH, USA","correspondingAuthor":false,"prefix":"","firstName":"Alireza","middleName":"","lastName":"Motazedian","suffix":""},{"id":610862636,"identity":"94b85c39-5bb5-4e0b-8971-8b617cf5e5a0","order_by":2,"name":"Mostafa M. Rezaee","email":"","orcid":"","institution":"Department of Data Science, Bowling Green State University, Bowling Green, OH, USA","correspondingAuthor":false,"prefix":"","firstName":"Mostafa","middleName":"M.","lastName":"Rezaee","suffix":""},{"id":610862637,"identity":"9236ba09-2748-4029-b455-a9478d5f0e21","order_by":3,"name":"Fatemeh Sedaghat Jalil-Abadi","email":"","orcid":"","institution":"Department of Energy Engineering and Physics, Amirkabir University of Technology, Tehran, Iran","correspondingAuthor":false,"prefix":"","firstName":"Fatemeh","middleName":"Sedaghat","lastName":"Jalil-Abadi","suffix":""},{"id":610862638,"identity":"d73e012b-b279-4a59-9d63-ab7a3743b5ad","order_by":4,"name":"Mohammad Ghadri","email":"","orcid":"https://orcid.org/0009-0004-4513-3306","institution":"MIAE Department, Concordia University, Montreal, QC, Canada","correspondingAuthor":false,"prefix":"","firstName":"Mohammad","middleName":"","lastName":"Ghadri","suffix":""}],"badges":[],"createdAt":"2026-03-23 16:47:59","currentVersionCode":1,"declarations":{"humanSubjects":false,"vertebrateSubjects":false,"conflictsOfInterestStatement":false,"humanSubjectEthicalGuidelines":false,"humanSubjectConsent":false,"humanSubjectClinicalTrial":false,"humanSubjectCaseReport":false,"vertebrateSubjectEthicalGuidelines":false},"doi":"10.21203/rs.3.rs-9203079/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-9203079/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":105354195,"identity":"d9091e3b-f695-475a-bd46-cbf69c6a8143","added_by":"auto","created_at":"2026-03-25 06:27:20","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":293640,"visible":true,"origin":"","legend":"","description":"","filename":"main.pdf","url":"https://assets-eu.researchsquare.com/files/rs-9203079/v1_covered_1eaf6c86-2ce3-4af8-855f-6497a388047f.pdf"}],"financialInterests":"The authors declare no competing interests.","formattedTitle":"\u003cp\u003eDepleted-Wave Electric-Field Solver for Pulsed Bessel-Gaussian Type-II Frequency Doubling in KTP Crystals\u003c/p\u003e","fulltext":[],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":false,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":true,"highlight":"","institution":"Shahid Chamran University of Ahvaz","isAcceptedByJournal":false,"isAuthorSuppliedPdf":true,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":true,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Bessel-Gaussian beams, second-harmonic generation, type-II phase matching, depleted-wave model, KTP, pulsed SHG","lastPublishedDoi":"10.21203/rs.3.rs-9203079/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-9203079/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eThis toolkit addresses the computational modeling of three-dimensional time-dependent electric-field distributions in type-II pulsed Bessel-Gaussian second-harmonic generation within potassium titanyl phosphate (KTP) crystals, where depleted-wave dynamics dominate nonlinear optical conversion. The toolkit solves three coupled nonlinear wave equations for ordinary and extraordinary fundamental waves and extraordinary second-harmonic waves using finite difference methods in cylindrical coordinates. It accepts user-defined parameters including pulse energy, beam spot size, crystal dimensions, and optical properties as inputs, and produces spatiotemporal electric-field amplitude distributions and conversion efficiency profiles as outputs. A unified Fortran codebase provides reproducible simulation pipelines, parametric sweep capabilities for examining interaction length dependencies, and computational optimizations that enable execution on standard desktop systems despite the fine radial meshes required for Bessel-Gaussian beam fluctuations. The implementation reproduces previously published quantitative predictions, specifically confirming that for ωf ≈ 80 μm spot sizes and 0.8 J pulse energies, complete energy exchange occurs over approximately 5 mm, which validates the necessity of depleted-wave formalism for crystals exceeding this interaction length. The toolkit is available as an open-source GitHub repository and is released as version v1.0.1 under the MIT License, with a permanently archived distribution identified by DOI 10.5281/zenodo.17362587.\u003c/p\u003e","manuscriptTitle":"Depleted-Wave Electric-Field Solver for Pulsed Bessel-Gaussian Type-II Frequency Doubling in KTP Crystals","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2026-03-25 06:25:03","doi":"10.21203/rs.3.rs-9203079/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"1087aade-baf8-4c8b-8a9d-468a95fd269c","owner":[],"postedDate":"March 25th, 2026","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":64997604,"name":"Optics/Lasers"}],"tags":[],"updatedAt":"2026-03-25T06:25:03+00:00","versionOfRecord":[],"versionCreatedAt":"2026-03-25 06:25:03","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-9203079","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-9203079","identity":"rs-9203079","version":["v1"]},"buildId":"XKTyCvWXoU3ODBz1xrDgd","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}