ndispers phase matching

%23%20%2F%2F%2F%20script%0A%23%20requires-python%20%3D%20%22%3E%3D3.11%22%0A%23%20dependencies%20%3D%20%5B%22ndispers%3D%3D0.18.0%22%2C%20%22numpy%22%2C%20%22plotly%22%5D%0A%23%20%2F%2F%2F%0A%22%22%22Phase-matching%20calculator%20for%20sum-frequency%20generation.%0A%0APhase-matching%20angles%2C%20acceptance%20bandwidths%20(wavelength%2C%20angle%2C%20temperature)%2C%0Awalk-off%20and%20effective%20nonlinearity.%20Runs%20in%20the%20browser%20via%20marimo's%20WASM%0Aexport%3B%20see%20apps%2FREADME.md.%0A%22%22%22%0A%0Aimport%20marimo%0A%0A__generated_with%20%3D%20%220.24.0%22%0Aapp%20%3D%20marimo.App(width%3D%22medium%22%2C%20app_title%3D%22ndispers%20phase%20matching%22)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%0A%20%20%20%20return%20(mo%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_()%3A%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20import%20plotly.graph_objects%20as%20go%0A%20%20%20%20from%20plotly.subplots%20import%20make_subplots%0A%20%20%20%20from%20ndispers.helper%20import%20brentq%0A%20%20%20%20import%20ndispers%20as%20nd%0A%20%20%20%20import%20ndispers.media.crystals%20as%20_C%0A%0A%20%20%20%20%23%20Parameterisations%20kept%20in%20the%20package%20only%20to%20reproduce%20earlier%0A%20%20%20%20%23%20calculations%20are%20not%20offered%20in%20the%20apps%3A%20every%20LBO%20set%20except%20KK2018%2C%0A%20%20%20%20%23%20and%20the%20deprecated%20AGS_Takaoka1999%20(see%20docs%2Fvalidation.md%2C%20%22Module%20status%22).%0A%20%20%20%20_HIDDEN%20%3D%20lambda%20name%3A%20((name.startswith(%22LBO_%22)%20and%20not%20name.startswith(%22LBO_KK2018%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20or%20name%20%3D%3D%20%22AGS_Takaoka1999%22)%0A%20%20%20%20%23%20only%20media%20that%20can%20phase-match%3A%20one%20of%20the%20two%20angles%20must%20be%20free%0A%20%20%20%20XTALS%20%3D%20%7B%7D%0A%20%20%20%20for%20_name%20in%20sorted(n%20for%20n%20in%20dir(_C)%20if%20not%20n.startswith(%22_%22))%3A%0A%20%20%20%20%20%20%20%20_cls%20%3D%20getattr(_C%2C%20_name)%0A%20%20%20%20%20%20%20%20if%20not%20isinstance(_cls%2C%20type)%20or%20_HIDDEN(_name)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20continue%0A%20%20%20%20%20%20%20%20_x%20%3D%20_cls()%0A%20%20%20%20%20%20%20%20if%20%22var%22%20in%20(_x.theta_rad%2C%20_x.phi_rad)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20XTALS%5B_name%5D%20%3D%20_cls%0A%20%20%20%20return%20XTALS%2C%20brentq%2C%20go%2C%20make_subplots%2C%20nd%2C%20np%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(XTALS%2C%20mo%2C%20nd)%3A%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20f%22%22%22%0A%20%20%20%20%23%20Phase-matching%20calculator%0A%0A%20%20%20%20Collinear%20three-wave%20mixing%20%E2%80%94%20sum-frequency%20generation%2C%20or%20difference-frequency%0A%20%20%20%20generation%20%2F%20optical%20parametric%20amplification%20and%20oscillation%20%E2%80%94%20in%20%7Blen(XTALS)%7D%0A%20%20%20%20birefringent%20crystals%3A%20phase-matching%20angles%2C%20tuning%20curves%2C%20acceptance%0A%20%20%20%20bandwidths%2C%20walk-off%20and%20effective%20nonlinearity%2C%20computed%20in%20your%20browser.%20Powered%20by%0A%20%20%20%20%5Bndispers%20%7Bnd.__version__%7D%5D(https%3A%2F%2Fgithub.com%2Fakihiko-shimura%2Fndispers).%0A%20%20%20%20%22%22%22%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20process%20%3D%20mo.ui.radio(%0A%20%20%20%20%20%20%20%20options%3D%5B%22SFG%20%2F%20SHG%22%2C%20%22DFG%20%2F%20OPA%20%2F%20OPO%22%5D%2C%20value%3D%22SFG%20%2F%20SHG%22%2C%20label%3D%22Process%22%2C%20inline%3DTrue)%0A%20%20%20%20process%0A%20%20%20%20return%20(process%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(XTALS%2C%20mo%2C%20process)%3A%0A%20%20%20%20IS_DFG%20%3D%20process.value.startswith(%22DFG%22)%0A%20%20%20%20xtal%20%3D%20mo.ui.dropdown(options%3Dlist(XTALS)%2C%20value%3D%22BetaBBO_Eimerl1987%22%2C%20label%3D%22Crystal%22)%0A%20%20%20%20T%20%3D%20mo.ui.number(-50%2C%20300%2C%20value%3D25%2C%20step%3D1%2C%20label%3D%22Temperature%20(%C2%B0C)%22)%0A%20%20%20%20L%20%3D%20mo.ui.number(0.1%2C%20100.0%2C%20value%3D10.0%2C%20step%3D0.1%2C%20label%3D%22Crystal%20length%20(mm)%22)%0A%20%20%20%20%23%20SFG%3A%20the%20two%20inputs%0A%20%20%20%20shg%20%3D%20mo.ui.checkbox(value%3DTrue%2C%20label%3D%22Second%20harmonic%20(%CE%BB%E2%82%81%20%3D%20%CE%BB%E2%82%82)%22)%0A%20%20%20%20wl1%20%3D%20mo.ui.number(200.0%2C%204000.0%2C%20value%3D1064.0%2C%20step%3D0.1%2C%20label%3D%22%CE%BB%E2%82%81%20(nm)%22)%0A%20%20%20%20wl2%20%3D%20mo.ui.number(200.0%2C%204000.0%2C%20value%3D1064.0%2C%20step%3D0.1%2C%20label%3D%22%CE%BB%E2%82%82%20(nm)%22)%0A%20%20%20%20%23%20DFG%20%2F%20OPA%20%2F%20OPO%3A%20the%20pump%20and%20one%20of%20the%20two%20generated%20waves%0A%20%20%20%20wlp%20%3D%20mo.ui.number(200.0%2C%204000.0%2C%20value%3D800.0%2C%20step%3D0.1%2C%20label%3D%22%CE%BB_pump%20(nm)%22)%0A%20%20%20%20given%20%3D%20mo.ui.radio(options%3D%5B%22signal%22%2C%20%22idler%22%5D%2C%20value%3D%22signal%22%2C%20label%3D%22given%20wave%22%2C%20inline%3DTrue)%0A%20%20%20%20wlg%20%3D%20mo.ui.number(200.0%2C%2020000.0%2C%20value%3D1300.0%2C%20step%3D0.1%2C%20label%3D%22%CE%BB_signal%20or%20%CE%BB_idler%20(nm)%22)%0A%20%20%20%20return%20IS_DFG%2C%20L%2C%20T%2C%20given%2C%20shg%2C%20wl1%2C%20wl2%2C%20wlg%2C%20wlp%2C%20xtal%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(IS_DFG%2C%20L%2C%20T%2C%20given%2C%20mo%2C%20shg%2C%20wl1%2C%20wl2%2C%20wlg%2C%20wlp%2C%20xtal)%3A%0A%20%20%20%20if%20IS_DFG%3A%0A%20%20%20%20%20%20%20%20_mid%20%3D%20mo.vstack(%5Bwlp%2C%20wlg%2C%20given%5D%2C%20gap%3D0.25)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20_mid%20%3D%20mo.vstack(%5Bwl1%2C%20wl2%20if%20not%20shg.value%20else%20mo.md(%22%3Csub%3E%CE%BB%E2%82%82%20%3D%20%CE%BB%E2%82%81%3C%2Fsub%3E%22)%5D)%0A%20%20%20%20mo.hstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.vstack(%5Bxtal%5D%20%2B%20(%5B%5D%20if%20IS_DFG%20else%20%5Bshg%5D)%2C%20gap%3D0.25)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_mid%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.vstack(%5BT%2C%20L%5D%2C%20gap%3D0.25)%2C%0A%20%20%20%20%20%20%20%20%5D%2C%0A%20%20%20%20%20%20%20%20justify%3D%22start%22%2C%0A%20%20%20%20%20%20%20%20gap%3D0.75%2C%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(IS_DFG%2C%20T%2C%20XTALS%2C%20given%2C%20mo%2C%20shg%2C%20wl1%2C%20wl2%2C%20wlg%2C%20wlp%2C%20xtal)%3A%0A%20%20%20%20x%20%3D%20XTALS%5Bxtal.value%5D()%0A%20%20%20%20%23%20ndispers%20works%20in%20%C2%B5m%3B%20the%20interface%20is%20in%20nm.%20Internally%20everything%20is%0A%20%20%20%20%23%20the%20SFG%20of%20waves%201%20and%202%20into%203%3B%20for%20DFG%20%2F%20OPA%20%2F%20OPO%20those%20are%0A%20%20%20%20%23%20(signal%2C%20idler%2C%20pump)%2C%20so%20every%20cell%20below%20serves%20both%20processes.%0A%20%20%20%20if%20IS_DFG%3A%0A%20%20%20%20%20%20%20%20mo.stop(wlg.value%20%3C%3D%20wlp.value%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22*The%20signal%20or%20idler%20must%20be%20longer%20than%20the%20pump%20%E2%80%94%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%221%2F%CE%BB_p%20%3D%201%2F%CE%BB_s%20%2B%201%2F%CE%BB_i.*%22))%0A%20%20%20%20%20%20%20%20WL3%20%3D%20wlp.value%20*%201e-3%0A%20%20%20%20%20%20%20%20_g%20%3D%20wlg.value%20*%201e-3%0A%20%20%20%20%20%20%20%20_other%20%3D%20x.wl_idler(WL3%2C%20_g)%0A%20%20%20%20%20%20%20%20WL1%2C%20WL2%20%3D%20(_g%2C%20_other)%20if%20given.value%20%3D%3D%20%22signal%22%20else%20(_other%2C%20_g)%0A%20%20%20%20%20%20%20%20WAVES%20%3D%20(%22signal%22%2C%20%22idler%22%2C%20%22pump%22)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20WL1%20%3D%20wl1.value%20*%201e-3%0A%20%20%20%20%20%20%20%20WL2%20%3D%20WL1%20if%20shg.value%20else%20wl2.value%20*%201e-3%0A%20%20%20%20%20%20%20%20WL3%20%3D%201.0%20%2F%20(1.0%20%2F%20WL1%20%2B%201.0%20%2F%20WL2)%0A%20%20%20%20%20%20%20%20WAVES%20%3D%20(%22%CE%BB%E2%82%81%22%2C%20%22%CE%BB%E2%82%82%22%2C%20%22%CE%BB%E2%82%83%22)%0A%20%20%20%20T_C%20%3D%20T.value%0A%20%20%20%20ANGLE_KEY%20%3D%20%22theta%22%20if%20x.theta_rad%20%3D%3D%20%22var%22%20else%20%22phi%22%0A%20%20%20%20ANGLE_SYM%20%3D%20%22%CE%B8%22%20if%20ANGLE_KEY%20%3D%3D%20%22theta%22%20else%20%22%CF%86%22%20%20%20%20%20%20%23%20for%20titles%20and%20labels%0A%20%20%20%20return%20ANGLE_KEY%2C%20ANGLE_SYM%2C%20T_C%2C%20WAVES%2C%20WL1%2C%20WL2%2C%20WL3%2C%20x%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo%2C%20x)%3A%0A%20%20%20%20%23%20the%20docstring%20is%20preformatted%2C%20but%20its%20reference%20lines%20are%20long%3B%20render%20it%0A%20%20%20%20%23%20in%20a%20wrapping%20block%20rather%20than%20a%20code%20fence%20so%20it%20does%20not%20scroll%20sideways%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20'%3Cdiv%20style%3D%22white-space%3Apre-wrap%3B%20overflow-wrap%3Aanywhere%3B%20'%0A%20%20%20%20%20%20%20%20'font-family%3Aui-monospace%2Cmonospace%3B%20font-size%3A0.82em%3B%20line-height%3A1.5%22%3E'%0A%20%20%20%20%20%20%20%20%2B%20(x.__doc__%20or%20%22%22).strip().replace(%22%26%22%2C%20%22%26amp%3B%22).replace(%22%3C%22%2C%20%22%26lt%3B%22)%0A%20%20%20%20%20%20%20%20%2B%20%22%3C%2Fdiv%3E%22%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(HAS_PM%2C%20POLS%2C%20mo%2C%20x)%3A%0A%20%20%20%20if%20HAS_PM%3A%0A%20%20%20%20%20%20%20%20_md%20%3D%20%22%5Cn%5Cn%22.join(%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22%24%24n_%7B%7B%7B_p%7D%7D%7D%20%3D%20%7Bx.n_latex(_p)%7D%24%24%22%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_p%20in%20dict.fromkeys(POLS)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.accordion(%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Sellmeier%20equation%20as%20evaluated%22%3A%20mo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22What%20ndispers%20differentiates%20and%20evaluates%2C%20with%20the%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22coefficients%20already%20substituted%20%E2%80%94%20not%20a%20quotation%20from%20the%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22paper.%20**%CE%BB%20is%20in%20%C2%B5m%20and%20T%20in%20%C2%B0C%20in%20this%20expression**%2C%20unlike%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22the%20nm%20inputs%20above.%5Cn%5Cn%22%20%2B%20_md%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%7D)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.md(%22%22)%0A%20%20%20%20_out%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(ANGLE_KEY%2C%20ANGLE_SYM%2C%20IS_DFG%2C%20T_C%2C%20WAVES%2C%20WL1%2C%20WL2%2C%20WL3%2C%20mo%2C%20x)%3A%0A%20%20%20%20%23%20Type%20I%3A%20both%20inputs%20share%20a%20polarization.%20Type%20II%3A%20they%20differ.%0A%20%20%20%20TYPE_OF%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22ooe%22%3A%20%22I%20(negative)%22%2C%20%22eeo%22%3A%20%22I%20(positive)%22%2C%0A%20%20%20%20%20%20%20%20%22oee%22%3A%20%22II%20(negative)%22%2C%20%22eoe%22%3A%20%22II%20(negative)%22%2C%0A%20%20%20%20%20%20%20%20%22eoo%22%3A%20%22II%20(positive)%22%2C%20%22oeo%22%3A%20%22II%20(positive)%22%2C%0A%20%20%20%20%7D%0A%20%20%20%20_pm%20%3D%20x.pmAngles_sfg(WL1%2C%20WL2%2C%20T_C%2C%20deg%3DTrue)%0A%20%20%20%20SOLUTIONS%20%3D%20%7Bk%3A%20_pm%5Bk%5D%5BANGLE_KEY%5D%5B0%5D%20for%20k%20in%20TYPE_OF%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20isinstance(_pm.get(k)%2C%20dict)%20and%20_pm%5Bk%5D.get(ANGLE_KEY)%7D%0A%0A%20%20%20%20_rows%20%3D%20%22%5Cn%22.join(%0A%20%20%20%20%20%20%20%20f%22%7C%20%60%7Bk%7D%60%20%7C%20Type%20%7BTYPE_OF%5Bk%5D%7D%20%7C%20%22%0A%20%20%20%20%20%20%20%20%2B%20(f%22**%7BSOLUTIONS%5Bk%5D%3A.3f%7D%C2%B0**%20%7C%20%7BSOLUTIONS%5Bk%5D%20*%203.141592653589793%20%2F%20180%3A.5f%7D%20rad%22%0A%20%20%20%20%20%20%20%20%20%20%20if%20k%20in%20SOLUTIONS%20else%20%22%E2%80%94%20none%20%7C%20%E2%80%94%22)%0A%20%20%20%20%20%20%20%20%2B%20%22%20%7C%22%0A%20%20%20%20%20%20%20%20for%20k%20in%20TYPE_OF%0A%20%20%20%20)%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20f%22%22%22%0A%20%20%20%20%23%23%20Phase-matching%20angles%0A%0A%20%20%20%20%7B(f%22pump%20%7BWL3%20*%201e3%3A.2f%7D%20nm%20%E2%86%92%20signal%20%7BWL1%20*%201e3%3A.2f%7D%20nm%20%2B%20**idler%20%7BWL2%20*%201e3%3A.2f%7D%20nm**%22%0A%20%20%20%20%20%20if%20IS_DFG%20else%0A%20%20%20%20%20%20f%22%CE%BB%E2%82%81%20%3D%20%7BWL1%20*%201e3%3A.2f%7D%20nm%2C%20%CE%BB%E2%82%82%20%3D%20%7BWL2%20*%201e3%3A.2f%7D%20nm%20%E2%86%92%20**%CE%BB%E2%82%83%20%3D%20%7BWL3%20*%201e3%3A.2f%7D%20nm**%22)%7D%0A%20%20%20%20%26nbsp%3B%C2%B7%26nbsp%3B%20the%20varying%20angle%20of%20this%20crystal%20is%20**%7BANGLE_SYM%7D**%0A%0A%20%20%20%20%7C%20polarizations%20(%7B%22%2C%20%22.join(WAVES)%20if%20IS_DFG%20else%20%221%2C2%2C3%22%7D)%20%7C%20type%20%7C%20angle%20%7C%20%7C%0A%20%20%20%20%7C---%7C---%7C---%7C---%7C%0A%20%20%20%20%7B_rows%7D%0A%20%20%20%20%22%22%22%0A%20%20%20%20)%0A%20%20%20%20return%20SOLUTIONS%2C%20TYPE_OF%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(SOLUTIONS%2C%20T_C%2C%20WL1%2C%20mo%2C%20x)%3A%0A%20%20%20%20pmtype%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20options%3Dlist(SOLUTIONS)%20or%20%5B%22ooe%22%5D%2C%0A%20%20%20%20%20%20%20%20value%3D(list(SOLUTIONS)%20or%20%5B%22ooe%22%5D)%5B0%5D%2C%0A%20%20%20%20%20%20%20%20label%3D%22Interaction%20to%20analyse%22%2C%0A%20%20%20%20)%0A%20%20%20%20%23%20a%20medium%20with%20an%20extraordinary-ray%20Sellmeier%20equation%20only%20(SLN)%20cannot%0A%20%20%20%20%23%20birefringently%20phase-match%20at%20all%3A%20say%20that%2C%20rather%20than%20%22change%20%CE%BB%22%0A%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20x.n(WL1%2C%200.0%2C%20T_C%2C%20pol%3D%22o%22)%0A%20%20%20%20%20%20%20%20_has_o%20%3D%20True%0A%20%20%20%20except%20ValueError%3A%0A%20%20%20%20%20%20%20%20_has_o%20%3D%20False%0A%20%20%20%20_none%20%3D%20(%0A%20%20%20%20%20%20%20%20%22*No%20birefringent%20phase%20matching%3A%20this%20medium%20has%20a%20Sellmeier%20equation%20for%20the%20%22%0A%20%20%20%20%20%20%20%20%22extraordinary%20ray%20only%2C%20and%20every%20birefringent%20interaction%20needs%20an%20ordinary%20wave.%20%22%0A%20%20%20%20%20%20%20%20%22It%20is%20meant%20for%20**quasi-phase%20matching**%20(all%20waves%20extraordinary%2C%20d33)%20%E2%80%94%20see%20the%20%22%0A%20%20%20%20%20%20%20%20%22QPM%20section%20below.*%22%20if%20not%20_has_o%20else%0A%20%20%20%20%20%20%20%20%22*No%20phase-matching%20solution%20%E2%80%94%20change%20%CE%BB%2C%20T%20or%20crystal.*%22)%0A%20%20%20%20pmtype%20if%20SOLUTIONS%20else%20mo.md(_none)%0A%20%20%20%20return%20(pmtype%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(ANGLE_SYM%2C%20SOLUTIONS%2C%20T_C%2C%20np%2C%20pmtype%2C%20xtal)%3A%0A%20%20%20%20HAS_PM%20%3D%20bool(SOLUTIONS)%20and%20pmtype.value%20in%20SOLUTIONS%0A%20%20%20%20ANGLE_PM%20%3D%20np.radians(SOLUTIONS%5Bpmtype.value%5D)%20if%20HAS_PM%20else%20np.nan%0A%20%20%20%20POLS%20%3D%20tuple(pmtype.value)%20if%20HAS_PM%20else%20(%22o%22%2C%20%22o%22%2C%20%22e%22)%0A%20%20%20%20%23%20the%20state%20every%20plot%20below%20was%20computed%20at%0A%20%20%20%20PLOT_TITLE%20%3D%20(%0A%20%20%20%20%20%20%20%20f%22%7Bxtal.value%7D%20%20%20T%20%3D%20%7BT_C%3Ag%7D%20%C2%B0C%2C%20%7BANGLE_SYM%7D%20%3D%20%22%0A%20%20%20%20%20%20%20%20f%22%7BSOLUTIONS%5Bpmtype.value%5D%3A.3f%7D%C2%B0%2C%20pol%20%3D%20%7Bpmtype.value%7D%22%0A%20%20%20%20%20%20%20%20if%20HAS_PM%20else%20%22%22%0A%20%20%20%20)%0A%20%20%20%20PLOT_LAYOUT%20%3D%20dict(%0A%20%20%20%20%20%20%20%20title%3Ddict(font_size%3D11%2C%20x%3D0%2C%20xanchor%3D%22left%22%2C%20y%3D0.97%2C%20yanchor%3D%22top%22)%2C%0A%20%20%20%20%20%20%20%20width%3D560%2C%20margin%3Ddict(l%3D60%2C%20r%3D20%2C%20t%3D40%2C%20b%3D45)%2C%20font_size%3D11%2C%0A%20%20%20%20%20%20%20%20template%3D%22plotly_white%22%2C%20showlegend%3DFalse%2C%20modebar%3Ddict(orientation%3D%22v%22)%2C%0A%20%20%20%20)%0A%20%20%20%20return%20ANGLE_PM%2C%20HAS_PM%2C%20PLOT_LAYOUT%2C%20PLOT_TITLE%2C%20POLS%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20ANGLE_PM%2C%0A%20%20%20%20HAS_PM%2C%0A%20%20%20%20L%2C%0A%20%20%20%20PLOT_LAYOUT%2C%0A%20%20%20%20PLOT_TITLE%2C%0A%20%20%20%20POLS%2C%0A%20%20%20%20T_C%2C%0A%20%20%20%20WL1%2C%0A%20%20%20%20WL2%2C%0A%20%20%20%20make_subplots%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20np%2C%0A%20%20%20%20x%2C%0A)%3A%0A%20%20%20%20if%20HAS_PM%3A%0A%20%20%20%20%20%20%20%20_w%20%3D%20np.radians(0.15)%0A%20%20%20%20%20%20%20%20_a%20%3D%20np.linspace(ANGLE_PM%20-%20_w%2C%20ANGLE_PM%20%2B%20_w%2C%20400)%0A%20%20%20%20%20%20%20%20_dk%20%3D%20x.dk_sfg(WL1%2C%20WL2%2C%20_a%2C%20T_C%2C%20*POLS)%0A%20%20%20%20%20%20%20%20_sinc%20%3D%20x.pmFactor_sfg(WL1%2C%20WL2%2C%20_a%2C%20T_C%2C%20*POLS%2C%20L.value)%0A%0A%20%20%20%20%20%20%20%20_deg%20%3D%20np.degrees(_a)%0A%20%20%20%20%20%20%20%20_fig%20%3D%20make_subplots(rows%3D2%2C%20cols%3D1%2C%20shared_xaxes%3DTrue%2C%20vertical_spacing%3D0.06)%0A%20%20%20%20%20%20%20%20_fig.add_scatter(x%3D_deg%2C%20y%3D_dk%2C%20mode%3D%22lines%22%2C%20row%3D1%2C%20col%3D1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20hovertemplate%3D%22%25%7Bx%3A.4f%7D%C2%B0%3Cbr%3E%CE%94k%20%25%7By%3A.5g%7D%20rad%2F%C2%B5m%3Cextra%3E%3C%2Fextra%3E%22)%0A%20%20%20%20%20%20%20%20_fig.add_scatter(x%3D_deg%2C%20y%3D_sinc%2C%20mode%3D%22lines%22%2C%20row%3D2%2C%20col%3D1%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20hovertemplate%3D%22%25%7Bx%3A.4f%7D%C2%B0%3Cbr%3Esinc%C2%B2%20%25%7By%3A.4f%7D%3Cextra%3E%3C%2Fextra%3E%22)%0A%20%20%20%20%20%20%20%20_fig.add_hline(y%3D0%2C%20line%3Ddict(color%3D%22gray%22%2C%20width%3D0.8)%2C%20row%3D1%2C%20col%3D1)%0A%20%20%20%20%20%20%20%20_fig.add_vline(x%3Dnp.degrees(ANGLE_PM)%2C%20line%3Ddict(color%3D%22red%22%2C%20width%3D0.8%2C%20dash%3D%22dash%22))%0A%20%20%20%20%20%20%20%20_fig.update_yaxes(title_text%3D%22%CE%94k%20(rad%2F%C2%B5m)%22%2C%20row%3D1%2C%20col%3D1)%0A%20%20%20%20%20%20%20%20_fig.update_yaxes(title_text%3D%22sinc%C2%B2(%CE%94kL%2F2)%22%2C%20row%3D2%2C%20col%3D1)%0A%20%20%20%20%20%20%20%20_fig.update_xaxes(title_text%3D%22angle%20(deg)%22%2C%20row%3D2%2C%20col%3D1)%0A%20%20%20%20%20%20%20%20_fig.update_layout(**PLOT_LAYOUT%2C%20title_text%3DPLOT_TITLE%2C%20height%3D360)%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.vstack(%5Bmo.md(%22%23%23%20Phase%20mismatch%20and%20conversion%22)%2C%20_fig%5D)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.md(%22%22)%0A%20%20%20%20_out%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(np)%3A%0A%20%20%20%20def%20zero_curve(xs%2C%20z)%3A%0A%20%20%20%20%20%20%20%20%22%22%22(x%2C%20fractional%20index)%20points%20where%20z%20changes%20sign%20along%20its%20second%0A%20%20%20%20%20%20%20%20axis%2C%20the%20crossing%20interpolated%20linearly%20-%20a%20zero%20locus%20as%20a%20sparse%0A%20%20%20%20%20%20%20%20point%20set%20(a%20contour%20trace%20would%20ship%20the%20whole%20grid%20to%20the%20browser).%22%22%22%0A%20%20%20%20%20%20%20%20_px%2C%20_py%20%3D%20%5B%5D%2C%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20_i%2C%20_row%20in%20enumerate(z)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_ok%20%3D%20np.isfinite(_row)%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_j%20in%20np.nonzero(np.diff(np.signbit(_row))%20%26%20_ok%5B%3A-1%5D%20%26%20_ok%5B1%3A%5D)%5B0%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_f%20%3D%20_row%5B_j%5D%20%2F%20(_row%5B_j%5D%20-%20_row%5B_j%20%2B%201%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_px.append(xs%5B_i%5D)%3B%20_py.append(_j%20%2B%20_f)%0A%20%20%20%20%20%20%20%20return%20np.array(_px)%2C%20np.array(_py)%0A%0A%20%20%20%20return%20(zero_curve%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20ANGLE_PM%2C%0A%20%20%20%20ANGLE_SYM%2C%0A%20%20%20%20HAS_PM%2C%0A%20%20%20%20IS_DFG%2C%0A%20%20%20%20PLOT_LAYOUT%2C%0A%20%20%20%20POLS%2C%0A%20%20%20%20T_C%2C%0A%20%20%20%20WL1%2C%0A%20%20%20%20WL2%2C%0A%20%20%20%20WL3%2C%0A%20%20%20%20make_subplots%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20np%2C%0A%20%20%20%20x%2C%0A%20%20%20%20zero_curve%2C%0A)%3A%0A%20%20%20%20%23%20OPO%20%2F%20OPA%20tuning%20curves%3A%20the%20zero%20contour%20of%20dk%20over%20(angle%2C%20signal%0A%20%20%20%20%23%20frequency)%20and%20over%20(temperature%2C%20signal%20frequency)%20-%20one%20vectorised%0A%20%20%20%20%23%20evaluation%20each%2C%20every%20branch%20drawn%2C%20no%20root%20finder%20needed%20for%20a%20plot.%0A%20%20%20%20if%20IS_DFG%20and%20HAS_PM%3A%0A%20%20%20%20%20%20%20%20import%20re%20as%20_re%0A%20%20%20%20%20%20%20%20%23%20longest%20idler%20to%20draw%3A%20the%20crystal's%20transparency%20edge%20(the%20Sellmeier%0A%20%20%20%20%20%20%20%20%23%20validity%20range%20is%20usually%20much%20narrower%20than%20where%20an%20OPO%20actually%20tunes)%0A%20%20%20%20%20%20%20%20_doc%20%3D%20type(x).__doc__%20or%20%22%22%0A%20%20%20%20%20%20%20%20_sec%20%3D%20(_re.search(r%22Transparency%20range%5Cs*%3A(%5B%5E%5Cn%5D*)%22%2C%20_doc)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20or%20_re.search(r%22Validity%20range%5Cs*%5Cn%5Cs*-%2B%5Cs*%5Cn(.%7B0%2C400%7D)%22%2C%20_doc%2C%20_re.S))%0A%20%20%20%20%20%20%20%20_m%20%3D%20_sec%20and%20_re.search(r%22(%5Cd%2B%5C.%3F%5Cd*)%5Cs*(%3F%3A%C2%B5m%7Cum)%3F%5Cs*(%3F%3Ato%7C-%7C%E2%80%93)%5Cs*(%5Cd%2B%5C.%3F%5Cd*)%22%2C%20_sec.group(1))%0A%20%20%20%20%20%20%20%20_wl_i_max%20%3D%20min(float(_m.group(2))%2C%2020.0)%20if%20_m%20else%205.0%0A%20%20%20%20%20%20%20%20_wl_i_max%20%3D%20max(_wl_i_max%2C%201.2%20*%20WL2)%20%20%20%20%20%20%20%20%20%20%23%20the%20current%20idler%20must%20be%20inside%0A%20%20%20%20%20%20%20%20_nu_p%20%3D%201.0%20%2F%20WL3%0A%20%20%20%20%20%20%20%20_nu_s%20%3D%20np.linspace(_nu_p%20%2F%202%2C%20_nu_p%20-%201.0%20%2F%20_wl_i_max%2C%20400)%20%20%20%23%20even%20in%20frequency%0A%20%20%20%20%20%20%20%20_wl_s%20%3D%201.0%20%2F%20_nu_s%0A%20%20%20%20%20%20%20%20_wl_i%20%3D%201.0%20%2F%20(_nu_p%20-%20_nu_s)%0A%20%20%20%20%20%20%20%20_ang%20%3D%20np.linspace(0.0%2C%2090.0%2C%20361)%0A%20%20%20%20%20%20%20%20_T%20%3D%20np.linspace(-50.0%2C%20300.0%2C%20141)%0A%20%20%20%20%20%20%20%20with%20np.errstate(all%3D%22ignore%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_dk_ang%20%3D%20np.asarray(x.dk_dfg(WL3%2C%20_wl_s%5BNone%2C%20%3A%5D%2C%20np.radians(_ang)%5B%3A%2C%20None%5D%2C%20T_C%2C%20*POLS)%2C%20float)%0A%20%20%20%20%20%20%20%20%20%20%20%20_dk_T%20%3D%20np.asarray(x.dk_dfg(WL3%2C%20_wl_s%5BNone%2C%20%3A%5D%2C%20ANGLE_PM%2C%20_T%5B%3A%2C%20None%5D%2C%20*POLS)%2C%20float)%0A%20%20%20%20%20%20%20%20_fig%20%3D%20make_subplots(rows%3D1%2C%20cols%3D2%2C%20horizontal_spacing%3D0.12%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20subplot_titles%3D(f%22angle%20tuning%20at%20T%20%3D%20%7BT_C%3Ag%7D%20%C2%B0C%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22temperature%20tuning%20at%20%7BANGLE_SYM%7D%20%3D%20%7Bnp.degrees(ANGLE_PM)%3A.3f%7D%C2%B0%22))%0A%20%20%20%20%20%20%20%20_idx%20%3D%20np.arange(_wl_s.size)%0A%20%20%20%20%20%20%20%20for%20_col%2C%20_xs%2C%20_dk%20in%20((1%2C%20_ang%2C%20_dk_ang)%2C%20(2%2C%20_T%2C%20_dk_T))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_px%2C%20_pj%20%3D%20zero_curve(_xs%2C%20_dk)%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_y%2C%20_name%2C%20_color%20in%20((_wl_s%20*%201e3%2C%20%22signal%22%2C%20%22%231f77b4%22)%2C%20(_wl_i%20*%201e3%2C%20%22idler%22%2C%20%22%23d62728%22))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_fig.add_scatter(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20x%3D_px%2C%20y%3Dnp.interp(_pj%2C%20_idx%2C%20_y)%20if%20_px.size%20else%20%5B%5D%2C%20mode%3D%22markers%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20marker%3Ddict(color%3D_color%2C%20size%3D3)%2C%20row%3D1%2C%20col%3D_col%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20name%3D_name%2C%20showlegend%3D(_col%20%3D%3D%201)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20hovertemplate%3Df%22%25%7B%7Bx%3A.3f%7D%7D%3Cbr%3E%7B_name%7D%20%25%7B%7By%3A.1f%7D%7D%20nm%3Cextra%3E%3C%2Fextra%3E%22)%0A%20%20%20%20%20%20%20%20_fig.add_scatter(x%3D%5Bnp.degrees(ANGLE_PM)%5D%2C%20y%3D%5BWL1%20*%201e3%5D%2C%20mode%3D%22markers%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20marker%3Ddict(color%3D%22%231f77b4%22%2C%20size%3D8)%2C%20row%3D1%2C%20col%3D1%2C%20showlegend%3DFalse%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20hovertemplate%3D%22current%20signal%3Cextra%3E%3C%2Fextra%3E%22)%0A%20%20%20%20%20%20%20%20_fig.add_scatter(x%3D%5Bnp.degrees(ANGLE_PM)%5D%2C%20y%3D%5BWL2%20*%201e3%5D%2C%20mode%3D%22markers%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20marker%3Ddict(color%3D%22%23d62728%22%2C%20size%3D8)%2C%20row%3D1%2C%20col%3D1%2C%20showlegend%3DFalse%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20hovertemplate%3D%22current%20idler%3Cextra%3E%3C%2Fextra%3E%22)%0A%20%20%20%20%20%20%20%20_fig.add_scatter(x%3D%5BT_C%5D%2C%20y%3D%5BWL1%20*%201e3%5D%2C%20mode%3D%22markers%22%2C%20marker%3Ddict(color%3D%22%231f77b4%22%2C%20size%3D8)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20row%3D1%2C%20col%3D2%2C%20showlegend%3DFalse%2C%20hovertemplate%3D%22current%20signal%3Cextra%3E%3C%2Fextra%3E%22)%0A%20%20%20%20%20%20%20%20_fig.add_scatter(x%3D%5BT_C%5D%2C%20y%3D%5BWL2%20*%201e3%5D%2C%20mode%3D%22markers%22%2C%20marker%3Ddict(color%3D%22%23d62728%22%2C%20size%3D8)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20row%3D1%2C%20col%3D2%2C%20showlegend%3DFalse%2C%20hovertemplate%3D%22current%20idler%3Cextra%3E%3C%2Fextra%3E%22)%0A%20%20%20%20%20%20%20%20_fig.update_xaxes(title_text%3Df%22%7BANGLE_SYM%7D%20(deg)%22%2C%20range%3D%5B0%2C%2090%5D%2C%20row%3D1%2C%20col%3D1)%0A%20%20%20%20%20%20%20%20_fig.update_xaxes(title_text%3D%22T%20(%C2%B0C)%22%2C%20range%3D%5B-50%2C%20300%5D%2C%20row%3D1%2C%20col%3D2)%0A%20%20%20%20%20%20%20%20_fig.update_yaxes(title_text%3D%22wavelength%20(nm)%22%2C%20range%3D%5BWL3%20*%201e3%2C%20_wl_i_max%20*%201e3%5D%2C%20row%3D1%2C%20col%3D1)%0A%20%20%20%20%20%20%20%20_fig.update_yaxes(range%3D%5BWL3%20*%201e3%2C%20_wl_i_max%20*%201e3%5D%2C%20row%3D1%2C%20col%3D2)%0A%20%20%20%20%20%20%20%20_fig.update_layout(**%7B**PLOT_LAYOUT%2C%20%22showlegend%22%3A%20True%2C%20%22width%22%3A%20760%7D%2C%20height%3D340%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20legend%3Ddict(x%3D0.01%2C%20y%3D0.99%2C%20bgcolor%3D%22rgba(255%2C255%2C255%2C0.6)%22))%0A%20%20%20%20%20%20%20%20_fig.update_annotations(font_size%3D11)%20%20%20%20%20%20%20%20%20%20%23%20subplot%20titles%20at%20axis-label%20size%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.vstack(%5Bmo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%23%23%20Tuning%20curves%5Cn%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22Signal%20and%20idler%20that%20phase-match%20for%20pump%20%7BWL3%20*%201e3%3A.2f%7D%20nm%20and%20polarizations%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22%60%7B''.join(POLS)%7D%60%2C%20as%20the%20angle%20or%20the%20temperature%20is%20tuned%3B%20the%20dots%20are%20the%20present%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%22operating%20point.%20Curves%20are%20drawn%20up%20to%20the%20crystal's%20transparency%20edge%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22(%7B_wl_i_max%3Ag%7D%20%C2%B5m)%3B%20outside%20the%20Sellmeier%20validity%20range%20they%20extrapolate%20the%20fit.%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%22A%20medium%20without%20a%20temperature%20term%20gives%20a%20flat%20temperature%20curve.%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_fig%5D)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.md(%22%22)%0A%20%20%20%20_out%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(T_C%2C%20WL1%2C%20WL2%2C%20mo%2C%20np%2C%20x)%3A%0A%20%20%20%20%23%20the%20poling%20period%20that%20phase-matches%20the%20present%20wavelengths%20-%20the%0A%20%20%20%20%23%20default%20for%20the%20temperature-tuning%20panel%2C%20editable%20to%20%22my%20crystal%20has%20%CE%9B%20%3D%2030%20%C2%B5m%22%0A%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20QPM_NOW%20%3D%20float(x.qpm_period_sfg(WL1%2C%20WL2%2C%20np.pi%20%2F%202%2C%20T_C%2C%20%22e%22%2C%20%22e%22%2C%20%22e%22))%0A%20%20%20%20except%20ValueError%3A%0A%20%20%20%20%20%20%20%20QPM_NOW%20%3D%20np.nan%0A%20%20%20%20qpm_period%20%3D%20mo.ui.number(%0A%20%20%20%20%20%20%20%200.1%2C%2010000.0%2C%20value%3Dround(QPM_NOW%2C%202)%20if%20np.isfinite(QPM_NOW)%20else%2030.0%2C%20step%3D0.01%2C%0A%20%20%20%20%20%20%20%20label%3D%22Poling%20period%20%CE%9B%20(%C2%B5m)%20for%20the%20temperature%20panel%22)%0A%20%20%20%20return%20QPM_NOW%2C%20qpm_period%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20IS_DFG%2C%0A%20%20%20%20PLOT_LAYOUT%2C%0A%20%20%20%20QPM_NOW%2C%0A%20%20%20%20T_C%2C%0A%20%20%20%20WL1%2C%0A%20%20%20%20WL2%2C%0A%20%20%20%20WL3%2C%0A%20%20%20%20make_subplots%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20np%2C%0A%20%20%20%20qpm_period%2C%0A%20%20%20%20shg%2C%0A%20%20%20%20x%2C%0A%20%20%20%20zero_curve%2C%0A)%3A%0A%20%20%20%20%23%20Quasi-phase-matching%3A%20the%20first-order%20poling%20period%20for%20all%20waves%0A%20%20%20%20%23%20extraordinary%20along%20a%20principal%20axis%20(d33%20of%20LiNbO3%2C%20LiTaO3%2C%20KTP)%2C%20and%0A%20%20%20%20%23%20what%20a%20fixed%20period%20selects%20as%20the%20temperature%20is%20tuned.%20Drawn%20for%20every%0A%20%20%20%20%23%20crystal%20with%20an%20e-ray%3B%20meaningful%20for%20the%20ferroelectrics.%0A%20%20%20%20if%20np.isfinite(QPM_NOW)%3A%0A%20%20%20%20%20%20%20%20_T%20%3D%20np.linspace(-50.0%2C%20300.0%2C%20141)%0A%20%20%20%20%20%20%20%20_Lam%20%3D%20qpm_period.value%0A%20%20%20%20%20%20%20%20with%20np.errstate(all%3D%22ignore%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20IS_DFG%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20sweep%20the%20signal%20at%20fixed%20pump%3B%20the%20idler%20follows%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_nu_p%20%3D%201.0%20%2F%20WL3%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_nu%20%3D%20np.linspace(_nu_p%20%2F%202%2C%20_nu_p%20-%201.0%20%2F%20min(5.5%2C%2012%20*%20WL3)%2C%20300)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_wl_x%20%3D%201.0%20%2F%20_nu%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_wl_other%20%3D%201.0%20%2F%20(_nu_p%20-%20_nu)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_per%20%3D%20np.asarray(x.qpm_period_dfg(WL3%2C%20_wl_x%2C%20np.pi%20%2F%202%2C%20T_C%2C%20%22e%22%2C%20%22e%22%2C%20%22e%22)%2C%20float)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_res%20%3D%20np.abs(np.asarray(x.dk_dfg(WL3%2C%20_wl_x%5BNone%2C%20%3A%5D%2C%20np.pi%20%2F%202%2C%20_T%5B%3A%2C%20None%5D%2C%20%22e%22%2C%20%22e%22%2C%20%22e%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20float))%20-%202%20*%20np.pi%20%2F%20_Lam%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_xlab%2C%20_curves%20%3D%20%22signal%20(nm)%22%2C%20((%22signal%22%2C%20_wl_x%2C%20%22%231f77b4%22)%2C%20(%22idler%22%2C%20_wl_other%2C%20%22%23d62728%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_now%20%3D%20(WL1%2C%20WL2)%0A%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20sweep%20%CE%BB%E2%82%81%20(for%20SHG%20%CE%BB%E2%82%82%20%3D%20%CE%BB%E2%82%81%20follows%3B%20otherwise%20%CE%BB%E2%82%82%20is%20held)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_wl_x%20%3D%20np.linspace(WL1%20%2F%201.8%2C%20WL1%20*%201.8%2C%20300)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_wl_2%20%3D%20_wl_x%20if%20shg.value%20else%20np.full_like(_wl_x%2C%20WL2)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_per%20%3D%20np.asarray(x.qpm_period_sfg(_wl_x%2C%20_wl_2%2C%20np.pi%20%2F%202%2C%20T_C%2C%20%22e%22%2C%20%22e%22%2C%20%22e%22)%2C%20float)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_res%20%3D%20np.abs(np.asarray(x.dk_sfg(_wl_x%5BNone%2C%20%3A%5D%2C%20_wl_2%5BNone%2C%20%3A%5D%2C%20np.pi%20%2F%202%2C%20_T%5B%3A%2C%20None%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22e%22%2C%20%22e%22%2C%20%22e%22)%2C%20float))%20-%202%20*%20np.pi%20%2F%20_Lam%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_xlab%20%3D%20%22%CE%BB%E2%82%81%20(nm)%22%20%2B%20(%22%22%20if%20shg.value%20else%20f%22%2C%20%CE%BB%E2%82%82%20%3D%20%7BWL2%20*%201e3%3Ag%7D%20nm%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_curves%20%3D%20((%22%CE%BB%E2%82%81%22%2C%20_wl_x%2C%20%22%231f77b4%22)%2C)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_now%20%3D%20(WL1%2C)%0A%20%20%20%20%20%20%20%20_fig%20%3D%20make_subplots(rows%3D1%2C%20cols%3D2%2C%20horizontal_spacing%3D0.12%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20subplot_titles%3D(f%22period%20for%20each%20wavelength%2C%20T%20%3D%20%7BT_C%3Ag%7D%20%C2%B0C%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22temperature%20tuning%20at%20%CE%9B%20%3D%20%7B_Lam%3Ag%7D%20%C2%B5m%22))%0A%20%20%20%20%20%20%20%20_fig.add_scatter(x%3D_wl_x%20*%201e3%2C%20y%3D_per%2C%20mode%3D%22lines%22%2C%20row%3D1%2C%20col%3D1%2C%20line%3Ddict(color%3D%22%231f77b4%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20showlegend%3DFalse%2C%20hovertemplate%3D%22%25%7Bx%3A.1f%7D%20nm%3Cbr%3E%CE%9B%20%3D%20%25%7By%3A.2f%7D%20%C2%B5m%3Cextra%3E%3C%2Fextra%3E%22)%0A%20%20%20%20%20%20%20%20_fig.add_scatter(x%3D%5BWL1%20*%201e3%5D%2C%20y%3D%5BQPM_NOW%5D%2C%20mode%3D%22markers%22%2C%20marker%3Ddict(size%3D8%2C%20color%3D%22%23d62728%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20row%3D1%2C%20col%3D1%2C%20showlegend%3DFalse%2C%20hovertemplate%3D%22present%20wavelengths%3Cextra%3E%3C%2Fextra%3E%22)%0A%20%20%20%20%20%20%20%20_px%2C%20_pj%20%3D%20zero_curve(_T%2C%20_res)%0A%20%20%20%20%20%20%20%20_idx%20%3D%20np.arange(_wl_x.size)%0A%20%20%20%20%20%20%20%20for%20_name%2C%20_y%2C%20_color%20in%20_curves%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_fig.add_scatter(x%3D_px%2C%20y%3Dnp.interp(_pj%2C%20_idx%2C%20_y%20*%201e3)%20if%20_px.size%20else%20%5B%5D%2C%20mode%3D%22markers%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20marker%3Ddict(color%3D_color%2C%20size%3D3)%2C%20row%3D1%2C%20col%3D2%2C%20name%3D_name%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20hovertemplate%3Df%22%25%7B%7Bx%3A.1f%7D%7D%20%C2%B0C%3Cbr%3E%7B_name%7D%20%25%7B%7By%3A.1f%7D%7D%20nm%3Cextra%3E%3C%2Fextra%3E%22)%0A%20%20%20%20%20%20%20%20if%20abs(_Lam%20-%20QPM_NOW)%20%3C%200.005%20%2B%201e-9%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_fig.add_scatter(x%3D%5BT_C%5D%20*%20len(_now)%2C%20y%3D%5Bw%20*%201e3%20for%20w%20in%20_now%5D%2C%20mode%3D%22markers%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20marker%3Ddict(size%3D8%2C%20color%3D%5Bc%20for%20_%2C%20_%2C%20c%20in%20_curves%5D%5B%3Alen(_now)%5D)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20row%3D1%2C%20col%3D2%2C%20showlegend%3DFalse%2C%20hovertemplate%3D%22present%20operating%20point%3Cextra%3E%3C%2Fextra%3E%22)%0A%20%20%20%20%20%20%20%20_fig.update_xaxes(title_text%3D_xlab%2C%20row%3D1%2C%20col%3D1)%0A%20%20%20%20%20%20%20%20_fig.update_yaxes(title_text%3D%22poling%20period%20%CE%9B%20(%C2%B5m)%22%2C%20row%3D1%2C%20col%3D1)%0A%20%20%20%20%20%20%20%20_fig.update_xaxes(title_text%3D%22T%20(%C2%B0C)%22%2C%20range%3D%5B-50%2C%20300%5D%2C%20row%3D1%2C%20col%3D2)%0A%20%20%20%20%20%20%20%20_ymax%20%3D%20max(w.max()%20for%20_%2C%20w%2C%20_%20in%20_curves)%20*%201e3%0A%20%20%20%20%20%20%20%20_ymin%20%3D%20(WL3%20if%20IS_DFG%20else%20_wl_x.min())%20*%201e3%0A%20%20%20%20%20%20%20%20_fig.update_yaxes(title_text%3D%22wavelength%20(nm)%22%2C%20range%3D%5B_ymin%2C%20_ymax%5D%2C%20row%3D1%2C%20col%3D2)%0A%20%20%20%20%20%20%20%20_fig.update_layout(**%7B**PLOT_LAYOUT%2C%20%22showlegend%22%3A%20True%2C%20%22width%22%3A%20760%7D%2C%20height%3D320%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20title_text%3D%22QPM%2C%20e%2Be%E2%86%92e%20along%20the%20axis%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20(f%22%2C%20pump%20%7BWL3%20*%201e3%3A.2f%7D%20nm%22%20if%20IS_DFG%20else%20%22%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20legend%3Ddict(x%3D0.52%2C%20y%3D0.99%2C%20bgcolor%3D%22rgba(255%2C255%2C255%2C0.6)%22))%0A%20%20%20%20%20%20%20%20_fig.update_annotations(font_size%3D11)%0A%20%20%20%20%20%20%20%20_flat%20%3D%20float(np.ptp(_res%2C%20axis%3D0).max())%20%3C%201e-12%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.vstack(%5Bmo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%23%23%20Quasi-phase-matching%5Cn%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%22First-order%20poling%20period%20for%20all%20three%20waves%20extraordinary%2C%20propagating%20along%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%22a%20principal%20axis%20(the%20PPLN%20%2F%20PPSLT%20%2F%20PPKTP%20geometry%2C%20using%20d33)%3A%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22**%CE%9B%20%3D%20%7BQPM_NOW%3A.2f%7D%20%C2%B5m**%20for%20the%20present%20wavelengths.%20The%20right%20panel%20is%20what%20a%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%22crystal%20of%20the%20period%20entered%20below%20phase-matches%20as%20it%20is%20heated%20%E2%80%94%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%2B%20(%22the%20temperature-tuning%20curve%20of%20a%20PPLN%20OPO.%22%20if%20IS_DFG%20else%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22the%20fundamental%20a%20PPLN%20doubler%20(or%20SFG%20mixer%2C%20%CE%BB%E2%82%82%20held)%20accepts%20at%20each%20temperature.%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%2B%20(%22%20**Flat%20here%3A%20%CE%94k%20does%20not%20change%20with%20temperature%20for%20this%20medium%20%E2%80%94%20either%20its%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Sellmeier%20equation%20has%20no%20temperature%20term%2C%20or%20its%20dn%2FdT%20is%20the%20same%20for%20all%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22wavelengths%20and%20cancels%20between%20the%20three%20waves.%20SLN%20and%20SLT%20carry%20a%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22wavelength-dependent%20dn%2FdT.**%22%20if%20_flat%20else%20%22%22))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20qpm_period%2C%20_fig%5D)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.md(%22%22)%0A%20%20%20%20_out%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(brentq%2C%20np)%3A%0A%20%20%20%20def%20acceptance(factor%2C%20x0%2C%20span0%2C%20threshold%3D0.5)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Full%20width%20of%20factor(x)%20at%20%60threshold%60%2C%20centred%20on%20the%20peak%20at%20x0.%0A%0A%20%20%20%20%20%20%20%20factor%20is%20sinc%5E2%2C%20so%20it%20falls%20monotonically%20from%201%20either%20side%20of%20x0%0A%20%20%20%20%20%20%20%20until%20its%20first%20zero%3A%20bracket%20that%20and%20solve%2C%20rather%20than%20stepping.%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20def%20f(u)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20factor(u)%20-%20threshold%0A%0A%20%20%20%20%20%20%20%20widths%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20_sign%20in%20(%2B1%2C%20-1)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20span%20%3D%20span0%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_%20in%20range(80)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20f(x0%20%2B%20_sign%20*%20span)%20%3C%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20break%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20span%20*%3D%201.6%0A%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20return%20np.nan%0A%20%20%20%20%20%20%20%20%20%20%20%20widths.append(abs(brentq(f%2C%20x0%2C%20x0%20%2B%20_sign%20*%20span%2C%20xtol%3Dspan0%20*%201e-8)%20-%20x0))%0A%20%20%20%20%20%20%20%20return%20sum(widths)%0A%0A%20%20%20%20return%20(acceptance%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20threshold%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%200.1%2C%200.9%2C%20value%3D0.5%2C%20step%3D0.05%2C%20label%3D%22Acceptance%20threshold%22%2C%20show_value%3DTrue%0A%20%20%20%20)%0A%20%20%20%20%23%20scale%20from%20ndispers'%20native%20unit%20(%C2%B5m%20for%20wavelength%2C%20rad%20for%20angle)%0A%20%20%20%20WL_UNITS%20%3D%20%7B%22pm%22%3A%201e6%2C%20%22nm%22%3A%201e3%2C%20%22%C2%B5m%22%3A%201.0%7D%0A%20%20%20%20ANG_UNITS%20%3D%20%7B%22mdeg%22%3A%20180e3%20%2F%203.141592653589793%2C%20%22mrad%22%3A%201e3%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22%C2%B5rad%22%3A%201e6%2C%20%22deg%22%3A%20180%20%2F%203.141592653589793%7D%0A%20%20%20%20wl_unit%20%3D%20mo.ui.dropdown(options%3Dlist(WL_UNITS)%2C%20value%3D%22pm%22%2C%20label%3D%22Wavelength%20unit%22)%0A%20%20%20%20ang_unit%20%3D%20mo.ui.dropdown(options%3Dlist(ANG_UNITS)%2C%20value%3D%22mrad%22%2C%20label%3D%22Angle%20unit%22)%0A%20%20%20%20mo.hstack(%5Bthreshold%2C%20wl_unit%2C%20ang_unit%5D%2C%20justify%3D%22start%22%2C%20gap%3D0.75)%0A%20%20%20%20return%20ANG_UNITS%2C%20WL_UNITS%2C%20ang_unit%2C%20threshold%2C%20wl_unit%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20ANGLE_PM%2C%0A%20%20%20%20ANG_UNITS%2C%0A%20%20%20%20HAS_PM%2C%0A%20%20%20%20IS_DFG%2C%0A%20%20%20%20L%2C%0A%20%20%20%20POLS%2C%0A%20%20%20%20T_C%2C%0A%20%20%20%20WL1%2C%0A%20%20%20%20WL2%2C%0A%20%20%20%20WL3%2C%0A%20%20%20%20WL_UNITS%2C%0A%20%20%20%20acceptance%2C%0A%20%20%20%20ang_unit%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20np%2C%0A%20%20%20%20shg%2C%0A%20%20%20%20threshold%2C%0A%20%20%20%20wl_unit%2C%0A%20%20%20%20x%2C%0A)%3A%0A%20%20%20%20if%20HAS_PM%3A%0A%20%20%20%20%20%20%20%20_th%20%3D%20threshold.value%0A%20%20%20%20%20%20%20%20_L%20%3D%20L.value%0A%20%20%20%20%20%20%20%20_au%2C%20_as%20%3D%20ang_unit.value%2C%20ANG_UNITS%5Bang_unit.value%5D%0A%20%20%20%20%20%20%20%20_wu%2C%20_ws%20%3D%20wl_unit.value%2C%20WL_UNITS%5Bwl_unit.value%5D%0A%0A%20%20%20%20%20%20%20%20_d_ang%20%3D%20acceptance(%0A%20%20%20%20%20%20%20%20%20%20%20%20lambda%20a%3A%20float(x.pmFactor_sfg(WL1%2C%20WL2%2C%20a%2C%20T_C%2C%20*POLS%2C%20_L))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ANGLE_PM%2C%201e-6%2C%20_th)%0A%20%20%20%20%20%20%20%20_d_T%20%3D%20acceptance(%0A%20%20%20%20%20%20%20%20%20%20%20%20lambda%20t%3A%20float(x.pmFactor_sfg(WL1%2C%20WL2%2C%20ANGLE_PM%2C%20t%2C%20*POLS%2C%20_L))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20T_C%2C%200.01%2C%20_th)%0A%20%20%20%20%20%20%20%20if%20IS_DFG%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20pump%20fixed%3A%20tuning%20the%20signal%20moves%20the%20idler%20the%20other%20way%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20(energy%20conservation)%20-%20the%20parametric%20gain%20bandwidth%2C%20whose%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20leading%20term%20is%20the%20signal-idler%20group-velocity%20mismatch%0A%20%20%20%20%20%20%20%20%20%20%20%20_d_sig%20%3D%20acceptance(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lambda%20w%3A%20float(x.pmFactor_dfg(WL3%2C%20w%2C%20ANGLE_PM%2C%20T_C%2C%20*POLS%2C%20_L))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20WL1%2C%201e-6%2C%20_th)%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20signal%20fixed%2C%20pump%20tuned%20(idler%20follows)%3A%20the%20pump%20bandwidth%0A%20%20%20%20%20%20%20%20%20%20%20%20_d_pump%20%3D%20acceptance(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lambda%20w%3A%20float(x.pmFactor_dfg(w%2C%20WL1%2C%20ANGLE_PM%2C%20T_C%2C%20*POLS%2C%20_L))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20WL3%2C%201e-6%2C%20_th)%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20%CE%BB%E2%82%82%20held%20at%20its%20stated%20value%2C%20only%20%CE%BB%E2%82%81%20tuned%20%E2%80%94%20the%20SNLO%20convention%0A%20%20%20%20%20%20%20%20%20%20%20%20_d_wl1%20%3D%20acceptance(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lambda%20w%3A%20float(x.pmFactor_sfg(w%2C%20WL2%2C%20ANGLE_PM%2C%20T_C%2C%20*POLS%2C%20_L))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20WL1%2C%201e-6%2C%20_th)%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20both%20tuned%20together%3B%20meaningful%20for%20SHG%20of%20a%20broadband%20source%0A%20%20%20%20%20%20%20%20%20%20%20%20_d_wboth%20%3D%20acceptance(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lambda%20w%3A%20float(x.pmFactor_sfg(w%2C%20w%2C%20ANGLE_PM%2C%20T_C%2C%20*POLS%2C%20_L))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20WL1%2C%201e-6%2C%20_th)%20if%20shg.value%20else%20np.nan%0A%0A%20%20%20%20%20%20%20%20def%20_row(label%2C%20value%2C%20scale%2C%20unit)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%22%22value%20is%20in%20ndispers'%20native%20unit%3B%20L-product%20is%20per%20centimetre.%22%22%22%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20not%20np.isfinite(value)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20return%20f%22%7C%20%7Blabel%7D%20%7C%20%E2%80%94%20%7C%20%E2%80%94%20%7C%22%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20(f%22%7C%20%7Blabel%7D%20%7C%20**%7Bvalue%20*%20scale%3A.4g%7D**%20%7Bunit%7D%20%7C%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22**%7Bvalue%20*%20scale%20*%20_L%20*%200.1%3A.4g%7D**%20%7Bunit%7D%C2%B7cm%20%7C%22)%0A%0A%20%20%20%20%20%20%20%20%23%20(label%2C%20width%20in%20the%20chosen%20unit%2C%20unit)%20-%20the%20report%20cell%20reads%20this%20too%0A%20%20%20%20%20%20%20%20ACCEPT%20%3D%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20(%22angle%22%2C%20_d_ang%20*%20_as%2C%20_au)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20(%22temperature%22%2C%20_d_T%2C%20%22%C2%B0C%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20if%20IS_DFG%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20ACCEPT.append((%22signal%20(pump%20fixed%2C%20idler%20follows)%22%2C%20_d_sig%20*%20_ws%2C%20_wu))%0A%20%20%20%20%20%20%20%20%20%20%20%20ACCEPT.append((%22pump%20(signal%20fixed%2C%20idler%20follows)%22%2C%20_d_pump%20*%20_ws%2C%20_wu))%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20ACCEPT.append((%22%CE%BB%E2%82%81%20only%2C%20%CE%BB%E2%82%82%20fixed%22%2C%20_d_wl1%20*%20_ws%2C%20_wu))%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20shg.value%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ACCEPT.append((%22%CE%BB%E2%82%81%20and%20%CE%BB%E2%82%82%20tuned%20together%22%2C%20_d_wboth%20*%20_ws%2C%20_wu))%0A%20%20%20%20%20%20%20%20_rows%20%3D%20%5B_row(_lbl%2C%20_v%2C%201%2C%20_u)%20for%20_lbl%2C%20_v%2C%20_u%20in%20ACCEPT%5D%0A%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22%22%22%0A%20%20%20%20%23%23%20Acceptance%20bandwidths%0A%0A%20%20%20%20Full%20width%20at%20%7B_th%3A.2f%7D%20of%20sinc%C2%B2(%CE%94kL%2F2)%20for%20L%20%3D%20%7B_L%3Ag%7D%20mm%2C%20and%20the%0A%20%20%20%20length%20product%2C%20which%20is%20what%20is%20usually%20tabulated.%0A%0A%20%20%20%20%7C%20%7C%20full%20width%20at%20L%20%3D%20%7B_L%3Ag%7D%20mm%20%7C%20%C3%97%20L%20%7C%0A%20%20%20%20%7C---%7C---%7C---%7C%0A%20%20%20%20%22%22%22%20%2B%20%22%5Cn%22.join(%22%20%20%20%20%22%20%2B%20r%20for%20r%20in%20_rows)%20%2B%20%22%22%22%0A%0A%20%20%20%20%22%22%22%20%2B%20(%0A%20%20%20%20%22%22%22The%20**signal**%20row%20is%20the%20parametric%20gain%20bandwidth%3A%20the%20pump%20is%20fixed%20and%0A%20%20%20%20the%20idler%20follows%20the%20signal%20through%201%2F%CE%BB_i%20%3D%201%2F%CE%BB_p%20%E2%88%92%201%2F%CE%BB_s%2C%20so%20to%20first%20order%0A%20%20%20%20it%20is%20set%20by%20the%20signal%E2%80%93idler%20group-velocity%20mismatch%20%E2%80%94%20it%20opens%20up%20towards%0A%20%20%20%20degeneracy%2C%20and%20fully%20in%20a%20noncollinear%20geometry.%20This%20is%20a%20different%20quantity%0A%20%20%20%20from%20SNLO's%20SFG%20acceptance%2C%20which%20tunes%20one%20input%20with%20the%20other%20fixed.%0A%20%20%20%20%22%22%22%20if%20IS_DFG%20else%0A%20%20%20%20%22%22%22The%20wavelength%20row%20to%20compare%20against%20tabulated%20values%20is%20**%CE%BB%E2%82%81%20only**%20%E2%80%94%0A%20%20%20%20that%20is%20what%20SNLO%20reports.%20Tuning%20%CE%BB%E2%82%81%20and%20%CE%BB%E2%82%82%20together%20makes%20%CE%94k%20move%20twice%0A%20%20%20%20as%20fast%2C%20so%20that%20width%20is%20exactly%20half%3B%20it%20is%20the%20relevant%20one%20when%20a%0A%20%20%20%20single%20broadband%20beam%20supplies%20both%20photons.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20ACCEPT%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.md(%22%22)%0A%20%20%20%20_out%0A%20%20%20%20return%20(ACCEPT%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(ANGLE_PM%2C%20HAS_PM%2C%20POLS%2C%20T_C%2C%20WAVES%2C%20WL1%2C%20WL2%2C%20WL3%2C%20mo%2C%20np%2C%20x)%3A%0A%20%20%20%20if%20HAS_PM%3A%0A%20%20%20%20%20%20%20%20_rows%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20_lbl%2C%20_wl%2C%20_p%20in%20zip(WAVES%2C%20(WL1%2C%20WL2%2C%20WL3)%2C%20POLS)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_rho%20%3D%20float(x.woa_theta(_wl%2C%20ANGLE_PM%2C%20T_C%2C%20pol%3D_p))%0A%20%20%20%20%20%20%20%20%20%20%20%20_n%20%3D%20float(x.n(_wl%2C%20ANGLE_PM%2C%20T_C%2C%20pol%3D_p))%0A%20%20%20%20%20%20%20%20%20%20%20%20_rows.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%7C%20%7B_lbl%7D%20%3D%20%7B_wl%20*%201e3%3A.2f%7D%20nm%20%7C%20%7B_p%7D%20%7C%20%7B_n%3A.5f%7D%20%7C%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%7Bnp.degrees(_rho)%3A.4f%7D%C2%B0%20%7C%20%7B_rho%20*%201e3%3A.4g%7D%20mrad%20%7C%22%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%23%23%20Refractive%20indices%20and%20walk-off%5Cn%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%22at%20the%20phase-matching%20angle.%20The%20ordinary%20ray%20of%20a%20uniaxial%20crystal%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%22has%20no%20walk-off%2C%20and%20the%20three%20indices%20coinciding%20is%20the%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%22phase-matching%20condition%20itself.%5Cn%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%7C%20wave%20%7C%20pol%20%7C%20n%20%7C%20walk-off%20%7C%20%7C%5Cn%7C---%7C---%7C---%7C---%7C---%7C%5Cn%22%20%2B%20%22%5Cn%22.join(_rows)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.md(%22%22)%0A%20%20%20%20_out%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(ANGLE_PM%2C%20HAS_PM%2C%20POLS%2C%20T_C%2C%20WL1%2C%20WL2%2C%20np%2C%20x)%3A%0A%20%20%20%20%23%20For%20a%20uniaxial%20crystal%20the%20refractive%20index%20does%20not%20depend%20on%20%CF%86%20(ndispers%0A%20%20%20%20%23%20reports%20phi_rad%20as%20'arb')%2C%20but%20%7Cd_eff%7C%20does%20through%20sin%203%CF%86%2C%20cos%202%CF%86%2C%20...%2C%0A%20%20%20%20%23%20so%20the%20azimuthal%20cut%20is%20a%20free%20choice.%20A%20principal%20plane%20of%20a%20biaxial%0A%20%20%20%20%23%20crystal%20fixes%20one%20angle%3B%20deff_sfg%20takes%20None%20for%20that%20one.%0A%20%20%20%20HAS_DEFF%20%3D%20hasattr(x%2C%20%22deff_sfg%22)%0A%20%20%20%20PHI_FREE%20%3D%20HAS_DEFF%20and%20x.phi_rad%20%3D%3D%20%22arb%22%0A%0A%20%20%20%20def%20deff_args(angle%2C%20phi)%3A%0A%20%20%20%20%20%20%20%20%22%22%22(theta%2C%20phi)%20for%20deff_sfg%3A%20the%20swept%20angle%20plus%20the%20fixed%2Fchosen%20one.%22%22%22%0A%20%20%20%20%20%20%20%20if%20x.theta_rad%20%3D%3D%20%22var%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20(angle%2C%20phi%20if%20PHI_FREE%20else%20None)%0A%20%20%20%20%20%20%20%20return%20(None%2C%20angle)%0A%0A%20%20%20%20if%20HAS_PM%20and%20PHI_FREE%3A%0A%20%20%20%20%20%20%20%20_p%20%3D%20np.linspace(0%2C%20180%2C%201801)%0A%20%20%20%20%20%20%20%20_d%20%3D%20np.abs(x.deff_sfg(WL1%2C%20WL2%2C%20ANGLE_PM%2C%20np.radians(_p)%2C%20T_C%2C%20*POLS))%0A%20%20%20%20%20%20%20%20_peak%20%3D%20_d.max()%0A%20%20%20%20%20%20%20%20_best%20%3D%20_p%5B_d%20%3E%20_peak%20*%200.9999%5D%0A%20%20%20%20%20%20%20%20PHI_OPT%20%3D%20%5B_best%5B0%5D%5D%0A%20%20%20%20%20%20%20%20for%20_v%20in%20_best%5B1%3A%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_v%20-%20PHI_OPT%5B-1%5D%20%3E%201.0%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20PHI_OPT.append(_v)%0A%20%20%20%20%20%20%20%20DEFF_MAX%20%3D%20_peak%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20PHI_OPT%2C%20DEFF_MAX%20%3D%20%5B%5D%2C%20np.nan%0A%20%20%20%20return%20DEFF_MAX%2C%20HAS_DEFF%2C%20PHI_FREE%2C%20PHI_OPT%2C%20deff_args%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(PHI_FREE%2C%20PHI_OPT%2C%20mo)%3A%0A%20%20%20%20phi_cut%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%200%2C%20180%2C%20value%3Dint(round(PHI_OPT%5B-1%5D))%20if%20PHI_OPT%20else%200%2C%20step%3D1%2C%0A%20%20%20%20%20%20%20%20label%3D%22Azimuthal%20cut%20%CF%86%20(deg)%22%2C%20show_value%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20phi_cut%20if%20PHI_FREE%20else%20mo.md(%22%22)%0A%20%20%20%20return%20(phi_cut%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20ANGLE_PM%2C%0A%20%20%20%20DEFF_MAX%2C%0A%20%20%20%20HAS_DEFF%2C%0A%20%20%20%20HAS_PM%2C%0A%20%20%20%20PHI_FREE%2C%0A%20%20%20%20PHI_OPT%2C%0A%20%20%20%20PLOT_LAYOUT%2C%0A%20%20%20%20PLOT_TITLE%2C%0A%20%20%20%20POLS%2C%0A%20%20%20%20T_C%2C%0A%20%20%20%20WL1%2C%0A%20%20%20%20WL2%2C%0A%20%20%20%20deff_args%2C%0A%20%20%20%20go%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20np%2C%0A%20%20%20%20phi_cut%2C%0A%20%20%20%20x%2C%0A)%3A%0A%20%20%20%20%23%20closed%20forms%20per%20point%20group%20(walk-off%20angles%20%CF%81_m%20of%20the%20e-waves%20included)%3B%0A%20%20%20%20%23%20the%20library%20evaluates%20the%20tensor%20contraction%20directly%2C%20these%20are%20for%20reading%0A%20%20%20%20_EXPR%20%3D%20%7B%0A%20%20%20%20%20%20%20%20(%22Uniax_3m%22%2C%20%22ooe%22)%3A%20r%22d_%7B31%7D%5Csin(%5Ctheta%2B%5Crho_3)%20-%20d_%7B22%7D%5Ccos(%5Ctheta%2B%5Crho_3)%5Csin%203%5Cphi%22%2C%0A%20%20%20%20%20%20%20%20(%22Uniax_3m%22%2C%20%22oee%22)%3A%20r%22d_%7B22%7D%5Ccos(%5Ctheta%2B%5Crho_2)%5Ccos(%5Ctheta%2B%5Crho_3)%5Ccos%203%5Cphi%22%2C%0A%20%20%20%20%20%20%20%20(%22Uniax_3m%22%2C%20%22eoe%22)%3A%20r%22d_%7B22%7D%5Ccos(%5Ctheta%2B%5Crho_1)%5Ccos(%5Ctheta%2B%5Crho_3)%5Ccos%203%5Cphi%22%2C%0A%20%20%20%20%20%20%20%20(%22Uniax_32%22%2C%20%22ooe%22)%3A%20r%22d_%7B11%7D%5Ccos(%5Ctheta%2B%5Crho_3)%5Ccos%203%5Cphi%22%2C%0A%20%20%20%20%20%20%20%20(%22Uniax_32%22%2C%20%22oee%22)%3A%20r%22d_%7B11%7D%5Ccos(%5Ctheta%2B%5Crho_2)%5Ccos(%5Ctheta%2B%5Crho_3)%5Csin%203%5Cphi%22%2C%0A%20%20%20%20%20%20%20%20(%22Uniax_32%22%2C%20%22eoe%22)%3A%20r%22d_%7B11%7D%5Ccos(%5Ctheta%2B%5Crho_1)%5Ccos(%5Ctheta%2B%5Crho_3)%5Csin%203%5Cphi%22%2C%0A%20%20%20%20%20%20%20%20(%22Uniax_42m%22%2C%20%22ooe%22)%3A%20r%22d_%7B36%7D%5Csin(%5Ctheta%2B%5Crho_3)%5Csin%202%5Cphi%22%2C%0A%20%20%20%20%20%20%20%20(%22Uniax_42m%22%2C%20%22oee%22)%3A%20r%22d_%7B36%7D%5Csin(2%5Ctheta%2B%5Crho_2%2B%5Crho_3)%5Ccos%202%5Cphi%22%2C%0A%20%20%20%20%20%20%20%20(%22Uniax_42m%22%2C%20%22eoe%22)%3A%20r%22d_%7B36%7D%5Csin(2%5Ctheta%2B%5Crho_1%2B%5Crho_3)%5Ccos%202%5Cphi%22%2C%0A%20%20%20%20%20%20%20%20(%22Uniax_4mm%22%2C%20%22ooe%22)%3A%20r%22d_%7B31%7D%5Csin(%5Ctheta%2B%5Crho_3)%22%2C%0A%20%20%20%20%7D%0A%20%20%20%20_group%20%3D%20type(x).__mro__%5B1%5D.__name__%0A%20%20%20%20_pols%20%3D%20%22%22.join(POLS)%0A%20%20%20%20if%20HAS_PM%20and%20HAS_DEFF%3A%0A%20%20%20%20%20%20%20%20_phi%20%3D%20np.radians(phi_cut.value)%20if%20PHI_FREE%20else%20None%0A%20%20%20%20%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_at%20%3D%20abs(float(x.deff_sfg(WL1%2C%20WL2%2C%20*deff_args(ANGLE_PM%2C%20_phi)%2C%20T_C%2C%20*POLS)))%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20point%20group%202%20(BiBO)%3A%20theta%20and%20pi%20-%20theta%20phase-match%20alike%20but%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20d_eff%20differs%2C%20so%20sweep%20the%20full%20half-circle%20there%0A%20%20%20%20%20%20%20%20%20%20%20%20_amax%20%3D%20np.pi%20if%20_group%20%3D%3D%20%22Biax_2%22%20else%20np.pi%20%2F%202%0A%20%20%20%20%20%20%20%20%20%20%20%20_a%20%3D%20np.linspace(0.02%2C%20_amax%20-%200.02%2C%20200%20if%20_amax%20%3C%202%20else%20400)%0A%20%20%20%20%20%20%20%20%20%20%20%20_curve%20%3D%20np.abs(x.deff_sfg(WL1%2C%20WL2%2C%20*deff_args(_a%2C%20_phi)%2C%20T_C%2C%20*POLS))%0A%20%20%20%20%20%20%20%20%20%20%20%20_fig%20%3D%20go.Figure(go.Scatter(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20x%3Dnp.degrees(_a)%2C%20y%3D_curve%2C%20mode%3D%22lines%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20hovertemplate%3D%22%25%7Bx%3A.2f%7D%C2%B0%3Cbr%3E%7Cd_eff%7C%20%25%7By%3A.4g%7D%20pm%2FV%3Cextra%3E%3C%2Fextra%3E%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20_fig.add_vline(x%3Dnp.degrees(ANGLE_PM)%2C%20line%3Ddict(color%3D%22red%22%2C%20width%3D0.8%2C%20dash%3D%22dash%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_group%20%3D%3D%20%22Biax_2%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_fig.add_vline(x%3D180%20-%20np.degrees(ANGLE_PM)%2C%20line%3Ddict(color%3D%22red%22%2C%20width%3D0.8%2C%20dash%3D%22dot%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20_fig.update_layout(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20**PLOT_LAYOUT%2C%20height%3D260%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20title_text%3DPLOT_TITLE%20%2B%20(f%22%2C%20%CF%86%20%3D%20%7Bphi_cut.value%3A.0f%7D%C2%B0%22%20if%20PHI_FREE%20else%20%22%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20xaxis_title%3D%22angle%20(deg)%22%2C%20yaxis_title%3D%22%7Cd_eff%7C%20(pm%2FV)%22)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20_forbidden%20%3D%20float(np.max(_curve))%20%3C%201e-9%20*%20max(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20abs(x.d_sfg(il%2C%20WL1%2C%20WL2%2C%20T_C))%20for%20il%20in%20x._d_ref)%0A%20%20%20%20%20%20%20%20%20%20%20%20_opt%20%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22%20%26nbsp%3B%C2%B7%26nbsp%3B%20%7Cd_eff%7C%20peaks%20at%20**%CF%86%20%3D%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20%22%2C%20%22.join(f%22%7Bv%3A.0f%7D%C2%B0%22%20for%20v%20in%20PHI_OPT)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20f%22**%20with%20%7BDEFF_MAX%3A.4g%7D%20pm%2FV%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20PHI_OPT%20else%20%22%22%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20_expr%20%3D%20_EXPR.get((_group%2C%20_pols))%0A%20%20%20%20%20%20%20%20%20%20%20%20_formula%20%3D%20(f%22%24%24%7Cd_%5C%5Cmathrm%7B%7Beff%7D%7D%7C%20%3D%20%5C%5Cleft%7C%7B_expr%7D%5C%5Cright%7C%24%24%5Cn%5Cn%22%20if%20_expr%20else%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22d_eff%20is%20the%20contraction%20of%20the%20d%20tensor%20with%20the%20field%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22directions%20of%20the%20three%20waves%2C%20walk-off%20included.%5Cn%5Cn%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20the%20tensor%20components%2C%20at%20their%20reference%20wavelengths%20and%20scaled%20to%20these%0A%20%20%20%20%20%20%20%20%20%20%20%20_rows%20%3D%20%22%5Cn%22.join(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%7C%20%7Bil%7D%20%7C%20%7Bd0%3Ag%7D%20%7C%20%7Bwl1%20*%201e3%3Ag%7D%20%2B%20%7Bwl2%20*%201e3%3Ag%7D%20%7C%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%7Bx.d_sfg(il%2C%20WL1%2C%20WL2%2C%20T_C)%3A.4g%7D%20%7C%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20for%20il%2C%20(d0%2C%20wl1%2C%20wl2)%20in%20x._d_ref.items())%0A%20%20%20%20%20%20%20%20%20%20%20%20_out%20%3D%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22%23%23%20Effective%20nonlinearity%5Cn%5Cn%22%20%2B%20_formula%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20f%22At%20the%20phase-matching%20angle%7Bf'%20and%20%CF%86%20%3D%20%7Bphi_cut.value%3A.0f%7D%C2%B0'%20if%20PHI_FREE%20else%20''%7D%3A%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22**%7B_at%3A.4g%7D%20pm%2FV**%7B_opt%7D.%5Cn%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20(%22This%20interaction%20phase-matches%20but%20its%20d_eff%20vanishes%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22identically%3A%20the%20tensor%20has%20no%20component%20that%20couples%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22these%20polarizations%20in%20the%20%7Bx.plane%7D%20plane%2C%20so%20no%20light%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22is%20generated.%20Choose%20another%20interaction%20or%20plane.%5Cn%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20_forbidden%20else%20%22%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20(%22In%20a%20monoclinic%20crystal%20the%20mirror%20angle%20180%C2%B0%20%E2%88%92%20%CE%B8%20phase-matches%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22too%20but%20with%20a%20different%20d_eff%3A%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22**%7Babs(float(x.deff_sfg(WL1%2C%20WL2%2C%20*deff_args(np.pi%20-%20ANGLE_PM%2C%20_phi)%2C%20T_C%2C%20*POLS)))%3A.4g%7D%20pm%2FV**%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22(dotted%20line).%5Cn%5Cn%22%20if%20_group%20%3D%3D%20%22Biax_2%22%20else%20%22%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20%22%7C%20component%20%7C%20reference%20(pm%2FV)%20%7C%20measured%20at%20(nm)%20%7C%20scaled%20to%20these%20wavelengths%20%7C%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22%7C---%7C---%7C---%7C---%7C%5Cn%22%20%2B%20_rows%20%2B%20%22%5Cn%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%3Csub%3EReference%20values%20and%20how%20far%20to%20trust%20the%20Miller-rule%20scaling%3A%20%7Bx._d_note%7D%3C%2Fsub%3E%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_fig%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D)%0A%20%20%20%20%20%20%20%20except%20ValueError%20as%20_e%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_out%20%3D%20mo.md(f%22%23%23%20Effective%20nonlinearity%5Cn%5Cn*Not%20available%20for%20%60%7B_pols%7D%60%3A%20%7B_e%7D*%22)%0A%20%20%20%20elif%20HAS_PM%3A%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%23%23%20Effective%20nonlinearity%5Cn%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22*%60%7Btype(x).__name__%7D%60%20is%20centrosymmetric%3A%20no%20second-order%20nonlinearity.*%22%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.md(%22%22)%0A%20%20%20%20_out%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20ACCEPT%2C%0A%20%20%20%20ANGLE_KEY%2C%0A%20%20%20%20ANGLE_PM%2C%0A%20%20%20%20DEFF_MAX%2C%0A%20%20%20%20HAS_DEFF%2C%0A%20%20%20%20HAS_PM%2C%0A%20%20%20%20IS_DFG%2C%0A%20%20%20%20L%2C%0A%20%20%20%20PHI_FREE%2C%0A%20%20%20%20PHI_OPT%2C%0A%20%20%20%20POLS%2C%0A%20%20%20%20SOLUTIONS%2C%0A%20%20%20%20TYPE_OF%2C%0A%20%20%20%20T_C%2C%0A%20%20%20%20WL1%2C%0A%20%20%20%20WL2%2C%0A%20%20%20%20WL3%2C%0A%20%20%20%20deff_args%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20nd%2C%0A%20%20%20%20np%2C%0A%20%20%20%20phi_cut%2C%0A%20%20%20%20shg%2C%0A%20%20%20%20threshold%2C%0A%20%20%20%20x%2C%0A%20%20%20%20xtal%2C%0A)%3A%0A%20%20%20%20%23%20A%20plain-text%20report%20of%20everything%20above%2C%20in%20the%20spirit%20of%20SNLO's%20QMIX%0A%20%20%20%20%23%20output%20but%20carrying%20its%20own%20provenance%3A%20which%20parameterisation%2C%20which%0A%20%20%20%20%23%20papers%2C%20which%20version%20-%20so%20the%20numbers%20can%20be%20traced%20later.%0A%20%20%20%20import%20datetime%20as%20_dt%0A%0A%20%20%20%20def%20_section(doc%2C%20head)%3A%0A%20%20%20%20%20%20%20%20%22%22%22The%20lines%20of%20a%20docstring%20section%20(%60head%60%20underlined%20with%20---).%22%22%22%0A%20%20%20%20%20%20%20%20_lines%20%3D%20%5Bl.strip()%20for%20l%20in%20(doc%20or%20%22%22).splitlines()%5D%0A%20%20%20%20%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_i%20%3D%20_lines.index(head)%0A%20%20%20%20%20%20%20%20except%20ValueError%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20%5B%5D%0A%20%20%20%20%20%20%20%20_body%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20l%20in%20_lines%5B_i%20%2B%202%3A%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20l%20and%20set(l)%20%3C%3D%20%7B%22-%22%7D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_body.pop()%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20the%20previous%20line%20was%20the%20next%20heading%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20break%0A%20%20%20%20%20%20%20%20%20%20%20%20_body.append(l)%0A%20%20%20%20%20%20%20%20while%20_body%20and%20not%20_body%5B-1%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_body.pop()%0A%20%20%20%20%20%20%20%20return%20_body%0A%0A%20%20%20%20if%20HAS_PM%3A%0A%20%20%20%20%20%20%20%20_deg%20%3D%20np.degrees(ANGLE_PM)%0A%20%20%20%20%20%20%20%20_cut%20%3D%20f%22theta%20%3D%20%7B_deg%3A.3f%7D%20deg%22%20if%20ANGLE_KEY%20%3D%3D%20%22theta%22%20else%20f%22phi%20%3D%20%7B_deg%3A.3f%7D%20deg%22%0A%20%20%20%20%20%20%20%20if%20x.theta_rad%20!%3D%20%22var%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_cut%20%3D%20f%22theta%20%3D%20%7Bnp.degrees(x.theta_rad)%3Ag%7D%20deg%2C%20%22%20%2B%20_cut%0A%20%20%20%20%20%20%20%20elif%20PHI_FREE%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_cut%20%2B%3D%20f%22%2C%20phi%20%3D%20%7Bphi_cut.value%3Ag%7D%20deg%22%0A%20%20%20%20%20%20%20%20elif%20not%20isinstance(x.phi_rad%2C%20str)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_cut%20%2B%3D%20f%22%2C%20phi%20%3D%20%7Bnp.degrees(x.phi_rad)%3Ag%7D%20deg%22%0A%20%20%20%20%20%20%20%20_cut_note%20%3D%20%22%20%20%20(internal%20angles%3B%20phi%20is%20a%20free%20choice%20here)%22%20if%20PHI_FREE%20else%20%22%20%20%20(internal%20angles)%22%0A%20%20%20%20%20%20%20%20_pols%20%3D%20%22%22.join(POLS)%0A%20%20%20%20%20%20%20%20_type%20%3D%20TYPE_OF%5B_pols%5D%0A%0A%20%20%20%20%20%20%20%20_waves%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20_wn%20%3D%20(%22signal%22%2C%20%22idler%22%2C%20%22pump%22)%20if%20IS_DFG%20else%20(%22wl1%22%2C%20%22wl2%22%2C%20%22wl3%22)%0A%20%20%20%20%20%20%20%20for%20_lbl%2C%20_wl%2C%20_p%20in%20zip(_wn%2C%20(WL1%2C%20WL2%2C%20WL3)%2C%20POLS)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_f%20%3D%20lambda%20q%2C%20_wl%3D_wl%2C%20_p%3D_p%3A%20float(getattr(x%2C%20q)(_wl%2C%20ANGLE_PM%2C%20T_C%2C%20pol%3D_p))%0A%20%20%20%20%20%20%20%20%20%20%20%20_waves.append((_lbl%2C%20_wl%2C%20_p%2C%20_f(%22n%22)%2C%20_f(%22ng%22)%2C%20_f(%22GV%22)%2C%20_f(%22GVD%22)%2C%20_f(%22TOD%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_f(%22woa_theta%22)%20*%201e3%2C%20_f(%22dndT%22)))%0A%20%20%20%20%20%20%20%20%23%20delays%20are%20quoted%20against%20the%20wave%20that%20comes%20in%3A%20wl1%20for%20SFG%2C%20the%20pump%0A%20%20%20%20%20%20%20%20%23%20for%20DFG.%201%2Fv_g(a)%20-%201%2Fv_g(ref)%20%3E%200%20means%20a%20falls%20behind%20the%20reference.%0A%20%20%20%20%20%20%20%20_ref%20%3D%202%20if%20IS_DFG%20else%200%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20pump%20for%20DFG%2C%20wl1%20for%20SFG%0A%20%20%20%20%20%20%20%20_oth%20%3D%20%5Bi%20for%20i%20in%20(0%2C%201%2C%202)%20if%20i%20!%3D%20_ref%5D%0A%20%20%20%20%20%20%20%20_gvm%20%3D%20lambda%20a%2C%20b%3A%201%20%2F%20_waves%5Ba%5D%5B5%5D%20-%201%20%2F%20_waves%5Bb%5D%5B5%5D%20%20%20%23%20fs%2Fum%0A%20%20%20%20%20%20%20%20_GVMS%20%3D%20%5B(_oth%5B0%5D%2C%20_ref)%2C%20(_oth%5B1%5D%2C%20_ref)%2C%20(_oth%5B1%5D%2C%20_oth%5B0%5D)%5D%0A%20%20%20%20%20%20%20%20_L_um%20%3D%20L.value%20*%201e3%0A%0A%20%20%20%20%20%20%20%20%23%20%22-%20Point%20group%20%3A%203m%20(C3v)%22%20in%20the%20docstring%3B%20the%20group%20mixin's%20name%20otherwise%0A%20%20%20%20%20%20%20%20_pg%20%3D%20next((l.split(%22%3A%22%2C%201)%5B1%5D.split(%22%3B%22)%5B0%5D.strip()%20for%20l%20in%20(x.__doc__%20or%20%22%22).splitlines()%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20%22Point%20group%22%20in%20l)%2C%20None)%0A%20%20%20%20%20%20%20%20if%20_pg%20is%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_pg%20%3D%20next((c.__name__.split(%22_%22%2C%201)%5B1%5D%20for%20c%20in%20type(x).__mro__%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20c.__name__.startswith((%22Uniax_%22%2C%20%22Biax_%22)))%2C%20%22n%2Fa%22)%0A%20%20%20%20%20%20%20%20_lines%20%3D%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22ndispers%20phase-matching%20report%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22generated%20%7B_dt.date.today().isoformat()%7D%20%20by%20ndispers%20%7Bnd.__version__%7D%20%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%22(akihiko-shimura.github.io%2Fndispers%2Fphasematching)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22CRYSTAL%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22%20%20%7Bxtal.value%7D%20%20%20point%20group%20%7B_pg%7D%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%2B%20(f%22%2C%20%7Bx.plane%7D%20plane%22%20if%20x.plane%20!%3D%20%22arb%22%20else%20%22%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20for%20l%20in%20_section(x.__doc__%2C%20%22Ref%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20l%3A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20headings%20flush%2C%20references%20indented%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_lines.append((%22%20%20%22%20if%20l.endswith(%22%3A%22)%20else%20%22%20%20%20%20%22)%20%2B%20l)%0A%20%20%20%20%20%20%20%20_vr%20%3D%20_section(x.__doc__%2C%20%22Validity%20range%22)%0A%20%20%20%20%20%20%20%20if%20_vr%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_lines.append(%22%20%20validity%3A%20%20%20%22%20%2B%20%22%20%22.join(_vr))%0A%20%20%20%20%20%20%20%20_lines%20%2B%3D%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22CONDITIONS%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20(%22%20%20process%20%20%20%20%20%20%20%20DFG%20%2F%20OPA%20%2F%20OPO%20%20%201%2Fwl_p%20%3D%201%2Fwl_s%20%2B%201%2Fwl_i%22%20if%20IS_DFG%20else%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%22%20%20process%20%20%20%20%20%20%20%20SFG%20%201%2Fwl3%20%3D%201%2Fwl1%20%2B%201%2Fwl2%22%20%2B%20(%22%20%20%20(SHG%3A%20wl1%20%3D%20wl2)%22%20if%20shg.value%20else%20%22%22))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20(f%22%20%20pump%2C%20signal%2C%20idler%20%20%7BWL3%20*%201e3%3A.2f%7D%2C%20%7BWL1%20*%201e3%3A.2f%7D%2C%20%7BWL2%20*%201e3%3A.2f%7D%20nm%20%20%20(vacuum)%22%20if%20IS_DFG%20else%0A%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%20%20wl1%2C%20wl2%2C%20wl3%20%20%7BWL1%20*%201e3%3A.2f%7D%2C%20%7BWL2%20*%201e3%3A.2f%7D%2C%20%7BWL3%20*%201e3%3A.2f%7D%20nm%20%20%20(vacuum)%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22%20%20temperature%20%20%20%20%7BT_C%3Ag%7D%20degC%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22%20%20length%20%20%20%20%20%20%20%20%20%7BL.value%3Ag%7D%20mm%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22%20%20polarizations%20%20%7B'%20'.join(POLS)%7D%20%20%20(Type%20%7B_type%7D%3B%20order%20%3D%20%7B'%2C%20'.join(_wn)%7D)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22%20%20cut%20%20%20%20%20%20%20%20%20%20%20%20%7B_cut%7D%7B_cut_note%7D%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22PHASE%20MATCHING%20%20%20(all%20solutions%20at%20this%20T%20and%20these%20wavelengths%3B%20%7BANGLE_KEY%7D%20varies)%22%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20for%20_k%20in%20TYPE_OF%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_lines.append(f%22%20%20%7B_k%7D%20%20Type%20%7BTYPE_OF%5B_k%5D%3A%3C14s%7D%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20(f%22%7BANGLE_KEY%7D%20%3D%20%7BSOLUTIONS%5B_k%5D%3A.3f%7D%20deg%22%20if%20_k%20in%20SOLUTIONS%20else%20%22none%22))%0A%20%20%20%20%20%20%20%20_lines%20%2B%3D%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22AT%20THE%20SELECTED%20SOLUTION%20%20(%7B_pols%7D%2C%20%7BANGLE_KEY%7D%20%3D%20%7B_deg%3A.3f%7D%20deg)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%20%20wave%20%20%20%20%20wl(nm)%20%20pol%20%20%20%20%20%20n%20%20%20%20%20%20%20n_g%20%20%20%20v_g(um%2Ffs)%20%20GVD(fs%5E2%2Fmm)%20%20TOD(fs%5E3%2Fmm)%20%20walk-off(mrad)%20%20dn%2FdT(1%2FK)%22%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20for%20_w%20in%20_waves%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_lines.append(f%22%20%20%7B_w%5B0%5D%3A%3C6s%7D%7B_w%5B1%5D%20*%201e3%3A9.2f%7D%20%20%20%7B_w%5B2%5D%7D%20%20%20%7B_w%5B3%5D%3A8.5f%7D%20%7B_w%5B4%5D%3A8.5f%7D%20%20%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%7B_w%5B5%5D%3A9.6f%7D%20%20%20%7B_w%5B6%5D%3A11.2f%7D%20%20%20%7B_w%5B7%5D%3A11.1f%7D%20%20%20%7B_w%5B8%5D%3A12.3f%7D%20%20%20%20%7B_w%5B9%5D%3A10.3e%7D%22)%0A%20%20%20%20%20%20%20%20_lines.append(%22%20%20group-velocity%20mismatch%22)%0A%20%20%20%20%20%20%20%20_w1%20%3D%20max(len(_wn%5Bi%5D)%20for%20i%20in%20(0%2C%201%2C%202))%0A%20%20%20%20%20%20%20%20for%20_a%2C%20_b%20in%20_GVMS%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_v%20%3D%20_gvm(_a%2C%20_b)%0A%20%20%20%20%20%20%20%20%20%20%20%20_expr%20%3D%20f%221%2Fv_g(%7B_wn%5B_a%5D%7D)%20-%201%2Fv_g(%7B_wn%5B_b%5D%7D)%22%0A%20%20%20%20%20%20%20%20%20%20%20%20_lines.append(f%22%20%20%20%20%7B_expr%3A%3C%7B2%20*%20_w1%20%2B%2016%7Ds%7D%20%3D%20%7B_v%3A%2B.4f%7D%20fs%2Fum%20%20%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22(%7B_v%20*%20_L_um%3A%2B.0f%7D%20fs%20over%20L)%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20(%22%20%20%20%3C-%20sets%20the%20parametric%20gain%20bandwidth%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20IS_DFG%20and%20%7B_a%2C%20_b%7D%20%3D%3D%20%7B0%2C%201%7D%20else%20%22%22))%0A%20%20%20%20%20%20%20%20_lines%20%2B%3D%20%5B%22%22%2C%20f%22ACCEPTANCE%20%20(full%20width%20at%20%7Bthreshold.value%3A.2f%7D%20of%20sinc%5E2(dk%20L%2F2)%3B%20%20L%20%3D%20%7BL.value%3Ag%7D%20mm)%22%5D%0A%20%20%20%20%20%20%20%20_tag%20%3D%20%7B%22%CE%BB%E2%82%81%20only%2C%20%CE%BB%E2%82%82%20fixed%22%3A%20%22wl1%20only%2C%20wl2%20fixed%20%20%3C-%20what%20SNLO%20and%20tables%20report%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22%CE%BB%E2%82%81%20and%20%CE%BB%E2%82%82%20tuned%20together%22%3A%20%22wl1%20%3D%20wl2%20tuned%20together%20%20%20(broadband%20SHG)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22signal%20(pump%20fixed%2C%20idler%20follows)%22%3A%20%22signal%2C%20pump%20fixed%20%20%20(idler%20follows%3B%20gain%20bandwidth)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22pump%20(signal%20fixed%2C%20idler%20follows)%22%3A%20%22pump%2C%20signal%20fixed%20%20%20(idler%20follows)%22%7D%0A%20%20%20%20%20%20%20%20for%20_lbl%2C%20_v%2C%20_u%20in%20ACCEPT%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_name%20%3D%20_tag.get(_lbl%2C%20_lbl)%0A%20%20%20%20%20%20%20%20%20%20%20%20_u2%20%3D%20%7B%22%C2%B0C%22%3A%20%22degC%22%2C%20%22%C2%B5rad%22%3A%20%22urad%22%2C%20%22%C2%B5m%22%3A%20%22um%22%7D.get(_u%2C%20_u)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20np.isfinite(_v)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_lines.append(f%22%20%20%7B_name.split('%20%20%3C-')%5B0%5D.split('%20%20%20(')%5B0%5D%3A%3C26s%7D%20%7B_v%3A10.4g%7D%20%7B_u2%3A%3C5s%7D%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%7B_v%20*%20L.value%20*%200.1%3A10.4g%7D%20%7B_u2%7D.cm%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20(%22%20%20%20%3C-%20%22%20%2B%20_name.split(%22%3C-%20%22)%5B1%5D%20if%20%22%3C-%22%20in%20_name%20else%20%22%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20(%22%20%20%20%22%20%2B%20_name%5B_name.index(%22(%22)%3A%5D%20if%20%22(%22%20in%20_name%20else%20%22%22))%0A%20%20%20%20%20%20%20%20_lines%20%2B%3D%20%5B%22%22%2C%20%22EFFECTIVE%20NONLINEARITY%22%5D%0A%20%20%20%20%20%20%20%20if%20HAS_DEFF%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_phi%20%3D%20np.radians(phi_cut.value)%20if%20PHI_FREE%20else%20None%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_de%20%3D%20float(x.deff_sfg(WL1%2C%20WL2%2C%20*deff_args(ANGLE_PM%2C%20_phi)%2C%20T_C%2C%20*POLS))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_lines.append(f%22%20%20d_eff%20%3D%20%7B_de%3A%2B.4g%7D%20pm%2FV%20%20%20at%20%7B_cut%7D%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20(f%22%20%20%20(%7Cd_eff%7C%20max%20%7BDEFF_MAX%3A.4g%7D%20pm%2FV%20at%20phi%20%3D%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20%22%2C%20%22.join(f%22%7Bv%3A.0f%7D%22%20for%20v%20in%20PHI_OPT)%20%2B%20%22%20deg)%22%20if%20PHI_OPT%20else%20%22%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_lines.append(%22%20%20%22%20%2B%20%22%2C%20%22.join(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%7Bil%7D%20%3D%20%7Bx.d_sfg(il%2C%20WL1%2C%20WL2%2C%20T_C)%3A.4g%7D%22%20for%20il%20in%20x._d_ref)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20f%22%20pm%2FV%20%20%20(Miller-scaled%20to%20%7BWL1%20*%201e3%3Ag%7D%20%2B%20%7BWL2%20*%201e3%3Ag%7D%20nm%20from%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20%22%3B%20%22.join(f%22%7Bil%7D%20%3D%20%7Bd0%3Ag%7D%20%40%20%7Bw1%20*%201e3%3Ag%7D%2B%7Bw2%20*%201e3%3Ag%7D%22%20for%20il%2C%20(d0%2C%20w1%2C%20w2)%20in%20x._d_ref.items())%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2B%20%22)%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_lines.append(f%22%20%20note%3A%20%7Bx._d_note%7D%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20except%20ValueError%20as%20_e%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_lines.append(f%22%20%20not%20available%20for%20%7B_pols%7D%3A%20%7B_e%7D%22)%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_lines.append(%22%20%20centrosymmetric%20crystal%3A%20no%20second-order%20nonlinearity%22)%0A%20%20%20%20%20%20%20%20_lines%20%2B%3D%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22NOTES%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%20%20Angles%20are%20internal%20and%20measured%20from%20the%20optic%2Fdielectric%20axes%3B%20wavelengths%20are%20in%20vacuum.%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%20%20The%20sign%20of%20d_eff%20is%20a%20convention%3B%20walk-off%20is%20included%20via%20theta%20%2B%20rho.%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20(f%22%20%20Reproduce%3A%20%20nd.media.crystals.%7Bxtal.value%7D().pmAngles_dfg(%7BWL3%3Ag%7D%2C%20%7BWL1%3Ag%7D%2C%20%7BT_C%3Ag%7D%2C%20deg%3DTrue)%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20if%20IS_DFG%20else%0A%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%20%20Reproduce%3A%20%20nd.media.crystals.%7Bxtal.value%7D().pmAngles_sfg(%7BWL1%3Ag%7D%2C%20%7BWL2%3Ag%7D%2C%20%7BT_C%3Ag%7D%2C%20deg%3DTrue)%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20REPORT%20%3D%20%22%5Cn%22.join(_lines)%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22%23%23%20Report%5Cn%5CnEverything%20above%20as%20plain%20text%20%E2%80%94%20copy%20it%20into%20a%20lab%20notebook%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22or%20alongside%20the%20data%20it%20belongs%20to.%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.ui.code_editor(REPORT%2C%20language%3D%22text%22%2C%20disabled%3DTrue%2C%20min_height%3D200)%2C%0A%20%20%20%20%20%20%20%20%5D)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20REPORT%20%3D%20%22%22%0A%20%20%20%20%20%20%20%20_out%20%3D%20mo.md(%22%22)%0A%20%20%20%20_out%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A